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The past few years have witnessed a rapid growth in the number and variety 
of applications of fuzzy logic. The applications range from consumer products 
such as cameras, camcorders, washing machines, and microwave ovens to 
industrial process control, medical instrumentation, decision-support systems, 
and portfolio select ion. 

To understand the reasons for the growing use of fuzzy logic it is necessary, 
first, to clarify what is meant by fuzzy logic. 

Fuzzy logic has two different meanings. I n a narrow sense, fuzzy logic is a 
logical system, which is an extension of multivalued logic. But in a wider 
sense— which is in predominant use today— fuzzy logic (FL) is almost 
synonymous with the theory of fuzzy sets, a theory which relates to cl asses of 
objects with unsharp boundaries in which membership is a matter of degree. 

I n this perspective, fuzzy logic in its narrow sense is a branch of FL. What is 
important to recognize is that, even in its narrow sense, the agenda of fuzzy 
logic is very different both in spirit and substance from the agendas of 
traditional multivalued logical systems. 

I n the Fuzzy Logic Tool box, fuzzy logic should be interpreted as FL, that is, 
fuzzy logic in its wide sense. The basic ideas underlying FL are explained very 
clearly and insightfully in the I ntroduction. What might be added is that the 
basic concept underlying FL is that of a linguistic variable, that is, a variable 
whose values are words rather than numbers. I n effect, much of FL may be 
viewed as a methodology for computing with words rather than numbers. 
Although words are inherently less precisethan numbers, their use is closer to 
human intuition. Furthermore, computing with words exploits the tolerance 
for imprecision and thereby lowers the cost of solution. 

Another basic concept in FL, which plays a central role in most of its 
applications, is that of a fuzzy if-then rule or, simply, fuzzy rule. Although 
rule-based systems have a long history of use in Al , what is missing in such 
systems is a machinery for dealing with fuzzy consequents and/or fuzzy 
antecedents. I n fuzzy logic, this machinery is provided by what is called the 
calculus of fuzzy rules. The calculus of fuzzy rules serves as a basis for what 
might be called the Fuzzy Dependency and Command Language (FDCL). 
Although FDCL is not used explicitly in Fuzzy Logic Tool box, it is effectively 
one of its principal constituents. I n this connection, what is important to 


recognize is that in most of the applications of fuzzy logic, a fuzzy logic solution 
is in reality a translation of a human solution into FDCL. 

What makes the Fuzzy Logic Tool box so powerful is the fact that most of 
human reasoning and concept formation is I inked to the use of fuzzy rules. By 
providing a systematic framework for computing with fuzzy rules, the F uzzy 
Logic Toolbox greatly amplifies the power of human reasoning. Further 
amplification results from the use of MATLAB and graphical user interfaces- 
areas in which The MathWorks has unparalleled expertise. 

A trend which is growing in visibility relates to the use of fuzzy logic in 
combination with neurocomputing and genetic algorithms. More generally, 
fuzzy logic, neurocomputing, and genetic algorithms may be viewed as the 
principal constituents of what might be cal led soft computing. Unlike the 
traditional, hard computing, soft computing is aimed at an accommodation 
with the pervasive imprecision of the real world. The guiding principle of soft 
computing is: Exploit the tolerance for imprecision, uncertainty, and partial 
truth to achieve tractabi I ity, robustness, and low solution cost. I n coming 
years, soft computing is likely to play an increasingly important role in the 
conception and design of systems whose M I Q (Machine IQ) is much higher than 
that of systems designed by conventional methods. 

Among various combi nations of methodologies in soft computing, theone which 
has highest visibility at this juncture is that of fuzzy logic and neurocomputing, 
leading to so-called neuro-fuzzy systems. Within fuzzy logic, such systems play 
a particularly important role in the induction of rules from observations. An 
effective method developed by Dr. Roger J ang for this purpose is called ANFIS 
(Adaptive Neuro-F uzzy I nference System). This method is an important 
component of the Fuzzy Logic Tool box. 

The F uzzy LogicToolbox is highly impressive in all respects. It makes fuzzy 
logic an effective tool for the conception and design of intelligent systems. The 
Fuzzy LogicToolbox is easy to master and convenient to use. And last, but not 
least important, it provides a reader -friendly and up-to-date introduction tothe 
methodology of fuzzy logic and its wide-ranging applications. 


Lotfi A. Zadeh 
Berkeley, CA 
J anuary 10, 1995 
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This section describes how to use the F uzzy Logic Tool box. It explains how to 
usethis guide and points you to additional books for toolbox installation 
information. 


What Is the Fuzzy Logic Toolbox? 

The F uzzy Logic Tool box is a collection of functions built on the MATLAB ® 
numeric computing environment. It provides tools for you to create and edit 
fuzzy inference systems within the framework of MATLAB, or if you prefer you 
can integrate your fuzzy systems into simulations with Simulink® or you can 
even build stand-alone C programs that call on fuzzy systems you build with 
MATLAB. This tool box relies heavily on graphical user interface (GUI ) tools to 
help you accomplish your work, although you can work entirely from the 
command line if you prefer. 

The tool box provides three categories of tools: 

• Command line functions 

• Graphical, interactive tools 

• Simulink blocks and examples 

The fi rst category of tools is made up of functions that you can cal I from the 
command line or from your own applications. Many of these functions are 
MATLAB M-files, series of MATLAB statements that implement specialized 
fuzzy logic algorithms. You can view the MATLAB code for these functions 
usi ng the statement 

type function_name 

You can change the way any toolbox function works by copying and renaming 
the M -file, then modifying your copy. You can also extend the tool box by adding 
your own M-files. 

Secondly, the tool box provides a number of interactive tools that let you access 
many of the functions through a GUI . Together, the GUI - based tools provide 
an environment for fuzzy i nference system design, analysis, and 
implementation. 

The third category of tools is a set of blocks for use with the Simulink 
simulation software. These are specifically designed for high speed fuzzy logic 
inference in the Simulink environment. 



How to Use This Guide 

If you are new to fuzzy logic, begin with Chapter 1, "I ntroduction."This 
chapter introduces the motivation behind fuzzy logic and leads you smoothly 
intothe tutorial. 

If you are an experienced fuzzy logic user, you may want to start at the 
beginning of Chapter 2, 'Tutorial, "to makesureyou are comfortable with the 
fuzzy logic terminology in the F uzzy LogicToolbox. If you just want an 
overview of each graphical tool and examples of specific fuzzy system tasks, 
turn directly to the section in Chapter 2 entitled, "Building Systems with the 
Fuzzy Logic Tool box." This section does not include information on the 
adaptivedata modeling application covered by the tool box function anfis.The 
basicfunctionality of thistool can be found in the section in Chapter 2 entitled, 
"anfis and the AN FIS Editor GUI." 

If you just want to start as soon as possibleand experiment, you can open an 
example system right away by typing 

fuzzy tipper 

This brings up the Fuzzy I nference System (FI S) editor for an example decision 
making problem that has to do with how to tip in a restaurant. 

All toolbox users should useChapter 3, "Reference, "for information on specific 
tools or functions. Reference descriptions includea synopsis of the function's 
syntax, as well as a complete explanation of options and operation. Many 
reference descriptions also include helpful examples, a description of the 
function’s algorithm, and references to additional reading material. For 
GUI-based tools, the descriptions include options for invoking the tool. 

Installation 

Toinstall thistool box on a workstation or a large machine, seeth el nstallation 
Guide for UNIX. To install the tool box on a PC or Macintosh, seethe 
Installation Guide for PC and Macintosh. 

To determine if the Fuzzy LogicToolbox is already installed on your system, 
check for a subdirectory named fuzzy within the main toolbox directory or 
folder. 



Typographical Conventions 


Before You Begin 


To Indicate 

This Guide Uses 

Example 

Example code 

Monospace type 
(Use Code tag . ) 

To assign the value 5 to A, 
enter 

A = 5 

Function 

Monospace type 

The cos function finds the 

names 

(Use Code tag . ) 

cosine of each array 
element. 

Function 

Monospace type for text 

The magic function uses 

syntax 

that must appear as 
shown. (Use Code tag.) 

Monospace italics for 
components you can 
replace with any variable. 

(UseCode-ital tag.) 

the syntax 

M =magic(n) 

Keys 

Boldface with an initial 
capital letter 

(Use Menu-Bodytexttag.) 

Press the Return key. 

Mathematical 

Variables in italics. 

This vector represents the 

expressions 

Functions, operators, and 
constants in standard 
type. (Use 

EquationVariables tag.) 

polynomial 

p =x 2 +2x + 3 

MATLAB 

output 

Monospace type 
(Use Code tag . ) 

MATLAB responds with 

A = 


5 


8 


To Indicate 

This Guide Uses 

Example 

Menu names, 

Boldface with an initial 

Choose the File menu. 

menu items, 

capital letter 


and controls 

(Use Menu-Bodytext tag.) 


New terms 

NCS italics 

An array is an ordered 


(Use Body text-ital tag.) 

collection of information. 


Before You Begin 


10 



I ntroduction 


What Is Fuzzy Logic? 1-2 

Why Use Fuzzy Logic? 1-5 

When Not to Use Fuzzy Logic 1-6 

What Can the Fuzzy Logic Tool box Do? 1-6 

An Introductory Example: Fuzzy vs. Non-Fuzzy .... 1-8 

The N on-Fuzzy Approach 1-9 

The F uzzy Approach 1-13 

Some Observations 1-14 


1 


Introduction 


What Is Fuzzy Logic? 

Fuzzy logic is all about the relative importance of precision: How important is 
it to be exactly right when a rough answer will do? All books on fuzzy logic 
begin with afewgood quotes on this very topic, and this is no exception. Here 
is what some clever people have said in the past: 

Precision is not truth. 

—Henri Matisse 

Someti mes the more measurabl e dri ves out the most i mportant. 

—Rene Du bos 

Vagueness is no moreto bedoneaway with in theworld of logic than friction in 
mechanics. 

—Charles Sanders Peirce 

I belie/e that nothing is unconditionally true, and hence I am opposed to B/ery 
statement of posi ti ve truth and B/ery man who makes i t. 

— H . L. Mencken 

So far as the laws of mathematics refer to reality, they are not certain. And so 
far as they are certai n, they do not refer to real i ty. 

—Albert Einstein 

As comp I exi ty ri ses, preci sestatemen ts I osemeani ng and meani ngful statemen ts 
lose precision. 

— Lotfi Zadeh 

Some pearls of folk wisdom also echo these thoughts: 

Don't lose sight of the forest for thetrees. 

Don't be penny wiseand pound foolish. 

The F uzzy Logic Tool box for use with MATLAB is a tool for solving problems 
with fuzzy logic. Fuzzy logic is a fascinating area of research because it does a 
good job of trading off between significance and precision— something that 
humans have been managing for a very long time. 

Fuzzy logic someti mes appears exotic or intimidating tothose unfamiliar with 
it, but once you become acquainted with it, it seems almost surprising that no 
one attempted it sooner. I n this sense fuzzy logic is both old and new because, 
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What Is Fuzzy Logic? 


although the modern and methodical science of fuzzy logic is still young, the 
concepts of fuzzy logic reach right down to our bones. 


Precision and Significance in the Real World 


A 1 500 kg mass 
is approaching 
your head at 
45.3 m/sec. 






Precision 


Significance 


Fuzzy logic is a convenient way to map an input space to an output space. This 
is the starting point for everything else, and the great emphasis here is on the 
word "convenient." 

What do I mean by mappi ng i nput space to output space? H ere are a few 
examples: You tell me how good your service was at a restaurant, and I ’ll tell 
you what the tip should be. You tell me how hot you want the water, and I 'll 
adjust the faucet valve to the right setting. You tell me how far away the 
subject of your photograph is, and I 'll focus the lens for you. You tell me how 
fast the car is going and how hard the motor is working, and I 'll shift the gears 
for you. 
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A graphical example of an input-output map is shown below. 


Input Space Output Space 

(all possible service ( a || possible tips) 

quality ratings) 



An input-output map for the tipping problem: 
"Given the quality of service, how much should I tip?" 


It's all just a matter of mapping inputstotheappropriateoutputs. Between the 
input and the output we’ll put a black box that does the work. What could go in 
the black box? Any number of things: fuzzy systems, linear systems, expert 
systems, neural networks, differential equations, interpolated 
multi-dimensional lookup tables, or even a spiritual advisor, just to name a few 
of the possible options. Clearly the list could go on and on. 

Of the dozens of ways to make the black box work, it turns out that fuzzy is 
often the very best way. Why should that be? As Lotfi Zadeh, who is considered 
to be the father of fuzzy logic, once remarked: "I n almost every case you can 
build the same product without fuzzy logic, but fuzzy is faster and cheaper." 
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Why Use Fuzzy Logic? 

Here is a list of general observations about fuzzy logic, 

• F uzzy logic is conceptually easy to understand. 

The mathematical concepts behind fuzzy reasoning are very simple. What 
makes fuzzy nice is the "natural ness" of its approach and not its far-reaching 
complexity. 

• Fuzzy logic is flexible. 

With any given system, it’s easy to massage it or layer more functionality on 
top of it without starting again from scratch. 

• F uzzy logic is tolerant of imprecise data. 

Everything is imprecise if you look closely enough, but more than that, most 
things are imprecise even on careful inspection. Fuzzy reasoning builds this 
understanding into the process rather than tacking it onto the end. 

• F uzzy logic can model nonlinear functions of arbitrary complexity. 

You can create a fuzzy system to match any set of input-output data. This 
process is made particularly easy by adaptive techniques like AN FIS 
(AdaptiveNeuro-Fuzzy I nference Systems), which are avail able in the Fuzzy 
Logic Tool box. 

• Fuzzy logic can be built on top of the experience of experts. 

I n direct contrast to neural networks, which take training data and generate 
opaque, impenetrable models, fuzzy logic lets you rely on the experience of 
people who already understand your system. 

• Fuzzy logic can be blended with conventional control techniques. 

Fuzzy systems don't necessarily replace conventional control methods. I n 
many cases fuzzy systems augment them and simplify their implementation. 

• F uzzy logic is based on natural language. 

The basis for fuzzy logic is the basis for human communication. This 
observation underpins many of the other statements about fuzzy logic. 

The last statement is perhaps the most important one and deserves more 
discussion. Natural language, that which is used by ordinary peopleon a daily 
basis, has been shaped by thousands of years of human history tobe convenient 
and efficient. Sentences written in ordinary language represent a triumph of 
efficient communication. We are generally unaware of this because ordinary 
language is, of course, something we use every day. Since fuzzy logic is built 


1-5 



atop the structures of qualitative description used in everyday language, fuzzy 
logic is easy to use. 


When Not to Use Fuzzy Logic 

Fuzzy logic is not a cure-all. When should you not use fuzzy logic? The safest 
statement is the first one made in this introduction: fuzzy logic is a convenient 
way to map an input space to an output space. If you find it's not convenient, 
try something else. If asimpler solution already exists, useit. Fuzzy logic is the 
codification of common sense— use common sense when you implement it and 
you will probably make the right decision. Many controllers, for example, do a 
fine job without using fuzzy logic. However, if you take the time to become 
familiar with fuzzy logic, you’ll see it can be a very powerful tool for dealing 
quickly and efficiently with imprecision and nonlinearity. 


What Can the Fuzzy Logic Toolbox Do? 

The F uzzy Logic Tool box allows you to do several things, but the most 
important thing it lets you do is create and edit fuzzy inference systems. You 
can create these systems using graphical tools or command-line functions, or 
you can generate them automatically using either clustering or adaptive 
neuro-fuzzy techniques. 

I f you have access to Si mul i nk, you can easi ly test your fuzzy system i n a block 
diagram simulation environment. 

The tool box also lets you run your own stand-alone C programs directly, 
without the need for Si mul ink. This is made possible by a stand-alone F uzzy 
I nference Engine that reads thefuzzy systems saved froma MATLAB session. 



What Is Fuzzy Logic? 


You can customize the stand-alone engine to build fuzzy inference into your 
own code. All provided code is ANSI compliant. 



Fuzzy 

Fuzz Y W | 1 Logic 

Inference 
System 

4 ? 

Stand-alone W ~ . 

Fuzzy Engine ■ User-written 

1 y S “ 1 M-files 


Other toolboxes 



Because of the integrated nature of MATLAB's environment, you can create 
your own tools to customize the F uzzy Logic Tool box or harness it with another 
toolbox, such astheControl System, Neural Network, or Opti mi zati on Tool box, 
to mention only a few of the possibilities. 
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Introduction 


An Introductory Example: Fuzzy vs. Non-Fuzzy 

A specific example would be helpful at this point. To illustrate the value of 
fuzzy logic, we'll show two different approaches to the same problem: linear 
and fuzzy. First we will work through this problem the conventional 
(non-fuzzy) way, writing MATLAB commands that spell out linear and 
piecewise-l inear relations. Then we’ll take a quick look at the same system 
using fuzzy logic. 

Consider the tipping problem: what is the "right" amount totip your 
waitperson? Here is a clear statement of the problem. 

The Basic Tipping Problem. Given a number between 0 and 10 that 
represents the quality of service at a restaurant (where 10 is excellent), what 
should the tip be? 


Cultural footnote: This problem is based on tipping as it is typically practiced 
in the United States. An average tip for a meal in theU.S. is 15%, though the 
actual amount may vary depending on the quality of the service provided. 
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An Introductory Example: Fuzzy vs. Non-Fuzzy 


The Non-Fuzzy Approach 

Let's start with the simplest possible relationship. Suppose that the tip always 
equals 15% of the total bill. 

tip = 0.15 


0.25 1 t t t t 

0 . 2 - 

0.15 

Q. 

0.1 - 

0.05 - 

0 I i ■ ■ ■ 

0 2 4 6 8 10 

service 

This doesn't really take into account the quality of the service, so we need to 
add a new term to the equation. Since service is rated on a scale of 0 to 10, we 
might have the tip go linearly from 5% if the service is bad to 25% if the service 
is excellent. Now our relation looks I ike this: 

tip=0 .20/1 0*service+0 . 05 
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So far so good. The formula does what we want it to do, and it's pretty 
straightforward. H owever, we may want theti p to reflect the qual ity of thefood 
as well. This extension of the problem is defined as follows: 

The Extended Tipping Problem. Given two sets of numbers between 0 and 10 
(where 10 is excel lent) that respectively represent the qual ity of the service and 
the quality of thefood at a restaurant, what should the tip be? 

Let's see how the formula will be affected now that we've added another 
variable. Suppose we try: 

tip = 0.20/20*(service+food)+0.05; 



I n this case, the results look pretty, but when you look at them closely, they 
don't seem quite right. Suppose you want the service to be a more i mportant 
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An Introductory Example: Fuzzy vs. Non-Fuzzy 


factor than the food quality. Let's say that the service will account for 80% of 
the overall tipping "grade" and the food will make up the other 20%. T ry: 

servRatio=0.8; 

tip=servRatio*(0.20/10*service+0.05) + ... 

( 1-servRatio ) * (0 . 20/ 1 0*f ood+0 . 05) ; 


0.25 

0.2 

. 9 - 0.15 

0.1 

0.05 

10 


The response is still somehow too uniformly linear. Suppose you want moreof 
a flat response in the middle, i.e., you want to give a 15% tip in general, and 
will depart from this plateau only if the service is exceptionally good or bad. 
This, in turn, means that those nice linear mappings no longer apply. We can 
still salvagethings by using a piecewise linear construction. Let's return tothe 
one-dimensional problem of just considering the service. You can string 
together a simple conditional statement using breakpoints I ike this: 

if service<3, 

tip=(0. 10/3) *service+0 . 05; 
elseif service<7, 
tip=0. 15; 

elseif service<=10, 

tip= (0. 1 0/3) * (service-7) +0. 15; 
end 
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Introduction 


The plot looks likethis. 



If we extend this to two dimensions, where we take food into account again, 
something likethis results: 

servRatio=0.8; 
if service<3, 

tip=( (0.10/3)*service+0.05)*servRatio + ... 

( 1-servRatio) * (0 .20/1 0*f ood+0 . 05) ; 
elseif service<7, 

tip=(0. 15)*servRatio + ... 

( 1-servRatio) *(0. 20/ 10*f ood+0. 05) ; 

else, 

tip=( (0 . 10/3) * (service-7)+0. 15) *servRatio + ... 

( 1-servRatio) *(0. 20/ 10*f ood+0. 05) ; 
end 
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An Introductory Example: Fuzzy vs. Non-Fuzzy 


Wow! The plot looks good, but the function is surprisingly complicated. It was 
a little tricky to code this correctly, and it's definitely not easy to modify this 
code in the future. Moreover, it's even less apparent how the algorithm works 
to someone who didn't witness the original design process. 


The Fuzzy Approach 

It would be nice if we could just capture the essentials of this problem, leaving 
aside all the factors that could be arbitrary. If we make a list of what really 
matters in this problem, we might end up with the following rule descriptions: 

1. If service is poor, then tip is cheap 

2. 1 f service is good, then tip is average 

3. 1 f service is excellent, then tip is generous 

The order in which the rules are presented here is arbitrary. It doesn’t matter 
which rules come first. If we wanted to include the food's effect on the tip, we 
might add the following two rules: 

4. If food is rancid, then tip is cheap 

5. 1 f food is delicious, then tip is generous 

I n fact, we can combine the two different lists of rules into one tight list of three 
rules like so: 

1. If service is poor or the food is rancid, then tipischeap 

2. If service is good, then tip is average 

3. Ifserviceis excellent or food is delicious, then tip is generous 

These three rules are the core of our solution. And coincidentally, we'vejust 
defined the rules for a fuzzy logic system. Now if we give mathematical 
meaning to the linguistic variables (what is an "average" tip, for example?) we 
would have a complete fuzzy inference system. Of course, there's a lot left to the 
methodology of fuzzy logicthat we're not mentioning right now, things like: 

• How are the rules all combined? 

• Howdol define mathematically what an "average" tip is? 

These are questions we provide detailed answers to in the next few chapters. 
The details of the method don’t really change much from problem to problem 
— the mechanics of fuzzy logic aren’t terribly complex. What matters is what 
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we've shown in this preliminary exposition: fuzzy is adaptable, simple, and 
easily applied. 



0.25 

0.2 

. 9 - 0.15 

0.1 

0.05 

10 


food 0 0 service 


Here is the picture associated with the fuzzy system that solves this problem. 
The picture above was generated by the three rules above. The mechanics of 
how fuzzy inference works is explained 'The Big Picture" on page 2-18, 
"Foundations of Fuzzy Logic" on page 2-20, and in "Fuzzy I nference Systems" 
on page 2-36. 1 n the "Building Systems with the Fuzzy LogicToolbox" on page 
2-45, the entire tipping problem is worked through using the graphical tools in 
the Fuzzy LogicToolbox. 

Some Observations 

Here are some observations about the example so far. We found a piecewise 
linear relation that solved the problem. It worked, but it was something of a 
nuisance to derive, and once we wrote it down as code, it wasn't very easy to 
interpret. On the other hand, the fuzzy system is based on some "common 
sense" statements. Also, we were able to add two more rules to the bottom of 
the list that influenced the shape of the overall output without needing to undo 
what had already been done. I n other words, the subsequent modification was 
pretty easy. 

Moreover, by using fuzzy logic rules, the maintenance of the structure of the 
algorithm decouples along fairly clean lines. The notion of an averagetip might 
change from day to day, city to city, cou ntry to cou ntry, but t he u nderl yi ng I ogi c 
isthesame: iftheserviceisgood,thetipshould be average. You can recalibrate 
the method quickly by simply shifting the fuzzy set that defines average 
without rewriting the fuzzy rules. 
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You can do this sort of thing with lists of piecewise linear functions, but there 
is a greater likelihood that recalibration will not be so quick and simple. 

For example, here is the piecewise linear tipping problem slightly rewritten to 
make it more generic. It performs the same function as before, only now the 
constants can be easily changed. 

% Establish constants 
lowTip=0.05; averTip=0. 15; highTip=0 . 25 ; 
tipRange=highTip-lowTip; 
badService=0; okayService=3; 
goodService=7; greatService=1 0 ; 
serviceRange=greatService-badService; 
badFood=0; greatFood=1 0 ; 
foodRange=greatFood-badFood; 

% If service is poor or food is rancid, tip is cheap 
if service<okayService , 

tip=( ( (averTip-lowTip) / (okayService-badService) ) ... 
*service+lowTip) *servRatio + ... 

(1-servRatio) * (tipRange/foodRange*food+lowTip) ; 

% If service is good, tip is average 
elseif service<goodService , 

tip=averTip*servRatio + (1-servRatio)* ... 
(tipRange/foodRange*food+lowTip) ; 

% If service is excellent or food is delicious, tip is generous 
else, 

tip=( ( (highTip-averTip) / ... 

(greatService-goodService) ) * ... 
(service-goodService)+averTip) *servRatio + ... 
(1-servRatio) * (tipRange/foodRange*food+lowTip) ; 
end 

Notice the tendency here, as with all code, for creeping generality to render the 
algorithm more and more opaque, threatening eventually to obscure it 
completely. What we're doing here isn't (shouldn't be!) that complicated. T rue, 
we can fight this tendency to be obscure by adding still more comments, or 
perhaps by trying to rewrite it in slightly more self-evident ways, but the 
medium is not on our side. 
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Thetrulyfascinatingthingtonoticeisthatifwe remove everyt h i n g except for 
three comments, what remain are exactly the fuzzy rules we wrote down 
before: 

% If service is poor or food is rancid, tip is cheap 

% If service is good, tip is average 

% If service is excellent or food is delicious, tip is generous 

If, as with a fuzzy system, the comment is identical with the code, think how 
much more I ikely your code is to have comments! F uzzy logic lets the language 
that's clearest to you, high level comments, also have meaning to the machine, 
which is why it's a very successful technique for bridging the gap between 
people and machines. 

Or think of it this way: by making the equations as simple as possible (linear) 
we make things simpler for the machine but more complicated for us. But 
really the limitation is no longer the computer— it's our mental model of what 
the computer is doing. We all know that computers have the ability to make 
things hopelessly complex; fuzzy logic is really about reclaiming the middle 
ground and letting the machine work with our preferences rather than the 
other way around. It's about time. 
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The Big Picture 

We'll start with a little motivation for where we are headed in this chapter. The 
poi nt of fuzzy logic is to map an i nput space to an output space, and the pri mary 
mechanism for doing this is a list of if-then statements called rules. All rules 
are evaluated in parallel, and the order of the rules is unimportant. The rules 
themselves are useful because they refer to variables and the adjectives that 
describethose variables. Beforewecan build asystemthat interprets rules, we 
have to define all the terms we plan on using and the adjectives that describe 
them. If we want to talk about how hot the water is, we need to define the range 
that the water's temperature can be expected to vary over as well as what we 
mean by the word hot. These are all things we’ll be discussing in the next 
several sections of the manual . The diagram below is something I ike a roadmap 
for the fuzzy inference process. It shows the general description of a fuzzy 
system on the left and a specific fuzzy system (the tipping example from the 
I ntroduction) on the right. 


The General Case... 


A Specific Example... 



To summarize the concept of fuzzy inference depicted in this figure, fuzzy 
inference is a method that interprets the values in the input vector and, based 
on some set of rules, assigns values to the output vector. 

This chapter is designed to guide you through the fuzzy logic process step by 
step by providing an introduction tothetheory and practice of fuzzy logic. The 
first three sections of this chapter are the most important— they move from 
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general to specific, first introducing underlying ideas and then discussing 
implementation details specific to the tool box. These three areas are 

• Foundations of fuzzy logic, which is an introduction to the general 
concepts. If you're already familiar with fuzzy logic, you may want to skip 
this section. 

• Fuzzy inference systems, which explains the specific methods of fuzzy 
inference used in the Fuzzy Logic Tool box. Si nee the field of fuzzy logic uses 
many terms that do not yet have standard interpretations, you should 
consi der readi ng this secti on j ust to become fami I i ar wi th the fuzzy i nference 
process as it is employed here. 

• Building systems with the Fuzzy Logic Toolbox, which goes into detail 
about how you build and edit a fuzzy system using this toolbox. This 
introduces the graphical user interface tools available in the Fuzzy Logic 
Toolbox and guides you through the construction of a complete fuzzy 

i nference system from start to fi nish. I f you just want to get up to speed as 
quickly as possible, start here. 

After this there are sections that touch on a variety of topics, such as Simulink 
use, automatic rulegeneration, and demonstrations. But from the point of view 
of getting to know the toolbox, these first three sections are the most crucial. 
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Foundations of Fuzzy Logic 

Everything is vague to a degreeyou do not realize til I you havetried tomakeit 
precise. —Bertrand Russell 

Fuzzy Sets 

F uzzy logic starts with the concept of a fuzzy set. A fuzzy set is a set without a 
crisp, clearly defined boundary. It can contain elements with only a partial 
degree of membershi p. 

To understand what a fuzzy set is, first consider what is meant by what we 
might call a classical set. A classical set is a container that wholly includes or 
wholly excludes any given element. For example, the set of days of the week 
unquestionably includes Monday, Thursday, and Saturday. It just as 
unquestionably excludes butter, liberty, and dorsal fins, and soon. 


Shoe 

Polish 


Butter 



Liberty 


Dorsal 

Fins 


Days of the week 


We call this set a classical set because it's been around for such a longtime. It 
was Aristotle who first formulated the Law of the Excluded M iddle, which says 
X must either be in set A or in set not -A. Another version runs I ike this: 

Of any subject, one thing must be either asserted or denied. 

Here is a restatement of the law with annotations: "Of any subject (say 
Monday), one thing (being a day of the week) must be either asserted or denied 
(I assert that Monday is a day of the week)." This law demands that opposites, 
the two categories A and not-A, should between them contain the entire 
universe. Everything falls into either one group or the other. There is nothing 
that is both a day of the week and not a day of the week. 
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Now consider the set of days comprising a weekend. Thediagram below is one 
attempt at classifying the weekend days. 


Shoe 

Polish 

Monday 

Butter 



Liberty 

Friday Thursday 

Dorsal 
Fins 


Days of the weekend 


Most would agree that Saturday and Sunday belong, but what about Friday? 
It "feels" like a part of the weekend, but somehow it seems like it should be 
technically excluded. So in thediagram above F riday tries its best tosit on the 
fence. Classical or "normal"sets wouldn't toleratethis kind of thing. Either 
you're in or you're out. Fluman experience suggests something different, 
though: fence sitting is a part of life. 

Of course we’re on tricky ground here, because we're starting to take individual 
perceptions and cultural background into account when we define what 
constitutes the weekend. But this is exactly the point. Even the dictionary is 
imprecise, defining the weekend as "the period from F riday night or Saturday 
to M onday morni ng." We're enteri ng the real m where sharp edged yes-no logic 
stops being helpful. Fuzzy reasoning becomes valuable exactly when we're 
talking about how people really perceive the concept "weekend" as opposed to 
a simple-minded classification useful for accounting purposes only. More than 
anything else, the following statement lays the foundations for fuzzy logic: 

/ n fuzzy logic, the truth of any statement becomes a matter of degree. 

Any statement can be fuzzy. The tool that fuzzy reasoning gives is the ability 
to reply to a yes-no question with a not-quite-yes-or-no answer. This is the kind 
of thing that humans do all the time (think how rarely you get a straight 
answer to a seemingly simple question) but it's a rather new trick for 
computers. 

Flow does it work? Reasoning in fuzzy logic is just a matter of generalizing the 
familiar yes-no (Boolean) logic. If we give "true" the numerical value of 1 and 
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"false" the numerical value of 0, we're saying that fuzzy logic also permits 
in-between values like 0.2 and 0.7453. For instance: 

Q: I s Saturday a weekend day? 

A: 1 (yes, or true) 

Q: Is Tuesday a weekend day? 

A: 0 (no, or false) 

Q: I s F riday a weekend day? 

A: 0.8 (for the most part yes, but not completely) 

Q: I s Sunday a weekend day? 

A: 0.95 (yes, but not quite as much as Saturday). 

Below on the left is a plot that shows the truth values for "weekend-ness" if we 
areforced to respond with an absolute yes or no response. On the right is a plot 
that shows the truth value for weekend-ness if we are allowed to respond with 
fuzzy in-between values. 



Days of the weekend two-valued membership Days of the weekend multivalued membership 


Technically, the representation on the right is from the domain of multivalued 
logic (or multivalent logic). If I ask the questi on "IsX a member of set A?" the 
answer might be yes, no, or any one of a thousand intermediate values in 
between. In other words, X might have partial membership in A. Multivalued 
logicstands in direct contrast to the more familiar concept of two-valued (or 
bivalent yes-no) logic. Two-valued logic has played a central role in the history 
of science since Aristotle first codified it, but the time has come for ittoshare 
the stage. 

To return to our example, now consider a continuous scale time plot of 
weekend-ness shown below. 
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Days of the weekend two-valued membership Days of the weekend multivalued membership 

By making the plot continuous, we're defining the degree to which any given 
instant belongs in the weekend rather than an entire day. I n the plot on the 
left, noticethat at midnight on Friday, just as thesecond hand sweeps past 12, 
the weekend-ness truth value jumps discontinuously from 0 to 1. This is one 
way to defi ne the weekend, and whi le it may be useful to an accountant, it 
doesn't really connect with our real-world experience of weekend-ness. 

The plot on the right shows a smoothly varying curve that accounts for the fact 
that al I of F riday, and, to a smal I degree, parts of Thursday, partake of the 
quality of weekend-ness and thus deserve partial membership in thefuzzy set 
of weekend moments. The curve that defines the weekend-ness of any instant 
in time is a function that maps the input space (time of the week) to the output 
space (weekend-ness). Specifically it is known as a membership function. We’ll 
discuss this in greater detail in the next section. 

As another example of fuzzy sets, consider the question of seasons. What 
season is it right now? I n the northern hemisphere, summer officially begins at 
the exact moment in the earth's orbit when the North Pole is pointed most 
directly toward the sun. It occurs exactly once a year, in latej une. Using the 
astronomical definitions for the season, we get sharp boundaries as shown on 
the left in the figure on the next page. But what we experience as the seasons 
varies more or less continuously as shown on the right below (in temperate 
northern hemisphere climates). 
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Membership Functions 

A membership function (MF) is a curvethat defines how each point in the input 
space is mapped to a membership value (or degree of membership) between 0 
and 1. The input space is sometimes referred to as the uni verse of discourse, a 
fancy name for a simple concept. 

One of the most commonly used examples of a fuzzy set is the set of tall people. 

I n this case the uni verse of discourse is all potential heights, say from 3 feet to 
9 feet, and the word "tall" would correspond to a curvethat defines the degree 
towhich any person is tall. If the set of tall people is given the well-defined 
(crisp) boundary of a classical set, we might say all people taller than six feet 
are officially considered tall. But such a distinction is clearly absurd. It may 
make sense to consider the set of al I real numbers greater than six because 
numbers belong on an abstract plane, but when we want to talk about real 
people, it is unreasonable to call one person short and another one tall when 
they differ in height by the width of a hair. 


excellent! 


IB 

1 

1 

1 

You must 
be taller 
than this 
line to be 
considered 
TALL 



But if the kind of distinction shown aboveis unworkable, then what isthe right 
way to define the set of tall people? M uch as with our plot of weekend days, the 
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figure below shows a smoothly varying curve that passes from not-tall to tall. 
The output-axis is a number known as the membership value between 0 and 1. 
The curve is known as a membership function and is often given the 
designation of p. This curve defines the transition from not tall to tall. Both 
people are tall to some degree, but one is significantly less tall than the other. 


1.0 

degree of 
membership, p 


0.0 


sharp-edged 
membership 
function for 
TALL 


height 



tall (pi = 1 .0) 


not tall (pi = 0.0) 



definitely a tall 
person (pi = 0.95) 


really not very 
tall at all (pi = 0.30) 


Subjective interpretations and appropriate units are built right into fuzzy sets. 
If I say "She's tall," the membership function "tall" should already take into 
account whether I 'm referring to a six-year-old or a grown woman. Similarly, 
the units are included in the curve. Certainly it makes nosensetosay "Is she 
tall in inches or in meters?" 

Membership Functions in the Fuzzy Logic Toolbox 

The only condition a membership function must really satisfy is that it must 
vary between 0 and 1. The function itself can bean arbitrary curve whose 
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shape we can define as a function that suits us from the point of view of 
simplicity, convenience, speed, and efficiency. 

A classical set might be expressed as 

A = jx| x >6} 

A fuzzy set is an extension of a classical set. If X is the universe of discourse 
and its elements are denoted by x, then a fuzzy set A in X is defined as a set of 
ordered pairs: 

A = fx, |i A (x) | xe X} 

|i A (x) is called the membership function (or MF) of x in A. The membership 
function maps each element of X to a membership value between 0 and 1. 

The Fuzzy LogicToolbox includes 11 built-in membership function types. 
These 11 functions are, in turn, built from several basic functions: piecewise 
linear functions, the Gaussian distribution function, the sigmoid curve, and 
quadratic and cubic polynomial curves. For detailed information on any of the 
membership functions mentioned below, turn to Chapter 3, "Reference" . By 
convention, all membership functions have the letters mf at the end of their 
names. 

The simplest membership functions are formed using straight lines. Of these, 
the simplest isth etriangular membership function, and it has the function 
name trimf . It's nothing more than a collection of three points forming a 
triangle. The trapezoidal membership function, trapmf, has a flat top and 
really is just a truncated triangle curve. These straight line membership 
functions have the advantage of simplicity. 



trimf, P = [3 6 81 trapmf, P = fl 5 7 81 

trimf trapmf 
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Two membership functions are built on the Gaussian distribution curve: a 
simple Gaussian curve and a two-sided composite of two different Gaussian 
curves. The two functions are gaussmf and gauss2mf . 

The general i zed bd I membership function is specified by three parameters and 
has the function name gbellmf. The bell membership function has one more 
parameter than the Gaussian membership function, so it can approach a 
non-fuzzy set if the free parameter is tuned. Because of their smoothness and 
concise notation, Gaussian and bell membership functions are popular 
methods for specifyi ng fuzzy sets. Both of these curves have the advantage of 
being smooth and nonzero at all points. 




gauss2mf 



gbellmf 


Although the Gaussian membership functions and bell membership functions 
achieve smoothness, they are unable to specify asymmetric membership 
functions, which are important in certain applications. Next we define the 
sigmoidal membership function, which is either open left or right. Asymmetric 
and closed (i.e. not open to the left or right) membership functions can be 
synthesized using two sigmoidal functions, so in addition to the basic sigmf , we 
also have the difference between two sigmoidal functions, dsigmf, and the 
product of two sigmoidal functions psigmf . 




psiqmf, P = f2 3 -5 81 


psigmf 


Polynomial based curves account for several of the membership functions in 
the tool box. Three related membership functions aretheZ, S, and Pi curves, all 
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named because of their shape. The function zmf is the asymmetrical 
polynomial curve open totheleft, smf is the mirror-image function that opens 
to the right, and pimf is zero on both extremes with a rise in the middle. 




zmf pimf smf 

There's a very wide selection to choose from when you're selecting your favorite 
membership function. And the Fuzzy Logic Tool box also allows you to create 
your own membership functions if you find this list too restrictive. On the other 
hand, if this list seems bewildering, just remember that you could probably get 
a I on g very wel I wi t h j u st on e or t wo types of membersh i p f u n ct i on s, f or exa m pi e 
the triangle and trapezoid functions. The selection is wide for those who want 
to explore the possibilities, but exotic membership functions are by no means 
required for perfectly good fuzzy i nference systems. Finally, remember that 
moredetails areavailableon all these functions in the reference section, which 
makes up the second half of this manual. 

Summary of Membership Functions 

• F uzzy sets describe vague concepts (fast runner, hot weather, weekend 
days). 

• A fuzzy set admits the possibility of partial membership in it. (F riday is sort 
of a weekend day, the weather is rather hot). 

• The degree an object belongs to a fuzzy set is denoted by a membership value 
between 0 and 1. (F riday is a weekend day to the degree 0.8). 

• A membership function associated with a given fuzzy set maps an input 
value to its appropriate membership value. 


Logical Operations 

We now know what's fuzzy about fuzzy logic, but what about the logic? 

The most important thing to realize about fuzzy logical reasoning is the fact 
that it is a superset of standard Boolean logic. I n other words, if we keep the 
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fuzzy values at their extremes of 1 (completely true), and 0 (completely false), 
standard logical operations will hold. As an example, consider the standard 
truth tables below: 


A BA and B 


0 

0 

0 

0 

1 

0 

1 

0 

0 

1 

1 

1 


AND 


A B 

A or B 

0 0 


0 

0 1 


1 

1 0 


1 

1 1 


1 


OR 


A 

not A 

0 

T 

1 

0 


NOT 


Now remembering that in fuzzy logicthetruth of any statement is a matter of 
degree, how will these truth tables be altered? The input values can be real 
numbers between 0 and 1. What function will preserve the results of the AND 
truth table (for example) and also extend to all real numbers between 0 and 1? 

One answer is the min operation. That is, resolve the statement A AND B, 
whereA and B are limited tothe range (0,1), by using the function min(A,B). 
Using the same reasoning, we can replace the OR operation with the max 
function, so that A OR 6 becomes equivalent to max(A,B). Finally, the 
operation NOT A becomes equivalent to the operation 1-A . Notice how the 
truth table above is completely unchanged by this substitution. 
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B 
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A B 

max(A,B) 
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T 

1 1 


T 


OR 
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0 


NOT 


Moreover, si nee there is a function behind the truth table rather than just the 
truth table itself, we can now consider values other than 1 and 0. 

The next figure uses a graph to show the same information. We've converted 
the truth table to a plot of two fuzzy sets applied together to create one fuzzy 
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set. The upper part of thefiguredisplays plots corresponding tothe two-valued 
truth tables above, whilethe lower part of thefiguredisplays how the 
operations work over a continuously varying range of truth values A and B 
according tothe fuzzy operations we've defined. 



AND OR NOT 

min(A,B) max(A,B) (1-A) 

Given these three functions, we can resolve any construction using fuzzy sets 
and the fuzzy logical operation AND, OR, and NOT. 

Additional Fuzzy Operators 

We've only defined here one particular correspondence between two-valued 
and multivalued logical operations for AND, OR, and NOT. This 
correspondence is by no means unique. 

I n more general terms, we're defining what are known as the fuzzy intersection 
or conjunction (AND), fuzzy union or disjunction (OR), and fuzzy complement 
(NOT). We have defined above what we'll call the classical operators for these 
functions: AND =min, OR - max, and NOT = additive complement. Typically 
most fuzzy logic applications make use of these operations and leave it at that. 
I n general, however, these functions are arbitrary to a surprising degree. The 
Fuzzy Logic Tool box uses the classical operator for the fuzzy complement as 
shown above, but also enables you to customize the AND and OR operators. 

The intersection of two fuzzy sets A and 6 is specified in general by a binary 
mapping T, which aggregates two membership functions as follows: 

PAnEiM =T(p a (x), p B (x)) 
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For example, the binary operator 7 may represent the multiplication of 
(o. A (x) and p B (x) . These fuzzy intersection operators, which are usually 
referred to as T-norm (T ri angular norm) operators, meet the following basic 
requirements. 

A T-norm operator is a binary mapping 7(- , ■) satisfying: 

boundary. 7(0, 0) = 0, T(a, 1) =T(1, a) -a 
monotonicity : T(a, b) <=T(c, d) if a <=cand b <=d 
commutativity: T(a, b) =T(b, a) 
associativity: T(a, T(b, c)) = T(T(a, b), c) 

The first requirement imposes the correct generalization tocrisp sets. The 
second requirement implies that a decrease in the membership values in A or 
B cannot produce an increase in the membership valuein A intersection 6. The 
third requirement indicates that the operator is indifferent to the order of the 
fuzzy sets to be combined. Finally, the fourth requirement allows ustotakethe 
intersection of any number of sets in any order of pairwise groupings. 

Like fuzzy intersection, the fuzzy union operator is specified in general by a 
binary mapping S: 

PaubM =S(u a (x), p b M) 

For example, the binary operator S can represent the addition of 

p A (x) and p B (x). These fuzzy union operators, which are often referred to as 

7-conorm (or S-norm) operators, must satisfy the foil owing basic requirements. 

A T-conorm (or S-norm ) operator is a binary mapping S(- , ■ ) satisfying: 

boundary: S(l, 1) =1, S(a, 0) =S( 0, a) =a 
monotonicity: S(a, b) <=S(c, d) if a <=cand b <=d 
commutativity: S(a, b) =S(b, a) 
associativity: S(a, S(b, c)) =S(S(a, b), c) 

Several parameterized 7-norms and dual 7-conorms have been proposed inthe 
past, such as those of Yager [Yag80], Dubois and Prade[Dub80], Schweizer and 
Sklar [Sch63], and Sugeno [Sug77], Each of these provides a way to vary the 
"gain"on the function so that it can be very restrictive or very permissive. 
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If-Then Rules 

Fuzzy sets and fuzzy operators are the subjects and verbs of fuzzy logic. These 
if-then rule statements are used to formulate the conditional statements that 
comprise fuzzy logic. 

A singlefuzzy if-then rule assumes the form 
if x is A then y is B 

where A and 6 are linguistic values defined by fuzzy sets on the ranges 
(universes of discourse) X and Y, respectively. The if-part of the rule "xisA" is 
called the an teceden tor premise, while the then-part of the rule "y is 6" is cal led 
the consequent or conclusion. An example of such a rule might be 

ifserviceis good then tip is average 

Note that good is represented as a number between 0 and 1, and sothe 
antecedent is an interpretation that returns a single number between 0 and 1. 
On theother hand, averageis represented as a fuzzy set, and so the consequent 
is an assignment that assigns the entire fuzzy set B tothe output variabley. In 
the if-then rule, the word "is" gets used in two entirely different ways 
depending on whether it appears in the antecedent or the consequent. I n 
MATLAB terms, this isthedistinction between a relational test using "="and 
a vari able assignment usi ng the "=" symbol . A less confusi ng way of wri ti ng the 
rule would be 

if service = good then tip =average 

I n general, the input to an if-then rule is the current value for the input 
variable (in this case, service) and the output is an entire fuzzy set (in this case, 
average). This set will later b eddfuzzified, assigning one value to the output. 
The concept of defuzzification is described in the next section, on page 2-41. 

I nterpreting an if-then rule involves distinct parts: first evaluating the 
antecedent (which involves fuzzi tying the input and applying any necessary 
fuzzy operators) and second applying that result to the consequent (known as 
implication). I n the case of two-valued or binary logic, if-then rules don't 
present much difficulty. If the premise is true, then the conclusion is true. If we 
relax the restrictions of two-valued logic and let the antecedent be a fuzzy 
statement, how does this reflect on the conclusion? The answer is a simple one: 
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if the antecedent is true to some degree of membership, then the consequent is 
also true to that same degree. In other words 

in binary logic: p -> q (p and q are either both true or both false) 

in fuzzy logic: 0.5 p -4 0.5 q (partial antecedents provide partial implication) 

The antecedent of a rule can have multiple parts: 

if sky is gray and wind is strong and baromdcer is falling, then ... 

in which case all parts of the antecedent are calculated simultaneously and 
resolved to a single number using the logical operators described in the 
preceding section. The consequent of a rule can also have multiple parts: 

if temperature is cold then hot water valveisopen and cold water valveisshut 

in which case all consequents are affected equally by the result of the 
antecedent. How is the consequent affected by the antecedent? The consequent 
specifies a fuzzy set be assigned to the output. The implication function then 
modifies that fuzzy set to the degree specified by the antecedent. The most 
common ways to modify the output fuzzy set are truncation using the min 
function (where the fuzzy set is "chopped off" as shown below) or scaling using 
the prod function (where the output fuzzy set is "squashed"). Both are 
supported by the Fuzzy Logic Tool box, but we use truncation for the examples 
in this section. 
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Antecedent 


Consequent 


e 


.A. 


> r 



1. Fuzzify 
inputs 


2. Apply 
OR operator 
(max) 


3. Apply 
implication 
operator (min) 


If service is excellent or food is delicious then tip = generous 



P(service==excellent) - 0 .0 P(food-=delicious) -0.7 


If ( 0.0 or 0.7 ) then tip = generous 



max(0.0, 0.7) =0.7 


If 


( 0.7 ) 


then tip = generous 


0.7 


A n 

generous I 1 


min(0. 7, generous) tipifuzzy) 


Summary of If-Then Rules 

Interpreting if-then rules is a three-part process. This process is explained in 
detail in the next section. 
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1 Fuzzify inputs: Resolve all fuzzy statements in the antecedent to a degree of 
membership between Oand 1. If there is only one part to the antecedent, this 
is the degree of support for the rule. 

2 Apply fuzzy operator to multiplepart antecedents : I f there are multi pie parts 
to the antecedent, apply fuzzy logic operators and resolve the antecedent to 
a single number between 0 and 1. This is the degree of support for the rule. 

3 Apply implication method: Use the degree of support for the entire rule to 
shape the output fuzzy set . The consequent of a fuzzy rule assi gns an enti re 
fuzzy set to the output. This fuzzy set is represented by a membership 
function that is chosen to indicate the qualities of the consequent. If the 
antecedent is only partially true, (i.e., is assigned a value less than 1), then 
the output fuzzy set is truncated according to the implication method. 

I n general, one rule by itself doesn't do much good. What's needed are two or 
more rules that can play off one another. The output of each rule is a fuzzy set. 
The output fuzzy sets for each rule are then aggregated into a single output 
fuzzy set. Finally the resulting set is dtfuzzified, or resolved to a single 
number. The next section shows how the whole process works from beginning 
to end for a particular type of fuzzy inference system called a M amdani type. 
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Fuzzy Inference Systems 

Fuzzy inference is the process of formulating the mapping from a given input 
to an output using fuzzy logic. The mapping then provides a basis from which 
decisions can be made, or patterns discerned. The process of fuzzy inference 
involves all of the pieces that are described in the previous sections: 
membership functions, fuzzy logic operators, and if-then rules. There are two 
types of fuzzy inference systems that can be implemented in the Fuzzy Logic 
Toolbox: M amdani-typeand Sugeno-type. These two types of i nference systems 
vary somewhat in the way outputs are determined. Descriptions of these two 
types of fuzzy i nference systems can be found i n the references, U an97, 
Mam75, Sug85], 

F uzzy inference systems have been successfully applied in fields such as 
automatic control, data classification, decision analysis, expert systems, and 
computer vision. Because of its multidisciplinary nature, fuzzy inference 
systems are associated with a number of names, such as fuzzy-rule-based 
systems, fuzzy expert systems, fuzzy modeling, fuzzy associative memory, 
fuzzy logic controllers, and simply (and ambiguously) fuzzy systems. Si nee the 
terms used to describe the various parts of the fuzzy inference process are far 
from standard, we will try to be as clear as possible about the different terms 
introduced in this section. 

Mamdani's fuzzy inference method is the most commonly seen fuzzy 
methodology. Mamdani's method was among the first control systems built 
using fuzzy set theory. It was proposed in 1975byEbrahimMamdani [Mam75] 
as an attempt to control a steam engine and boiler combination by synthesizing 
a set of linguistic control rules obtained from experienced human operators. 
Mamdani's effort was based on Lotfi Zadeh's 1973 paper on fuzzy algorithms 
for complex systems and decision processes [Zad73], Although the inference 
process we descri be i n the next few secti ons differs somewhat from the methods 
described in the original paper, the basic idea is much the same. 

Mamdani-type inference, as we have defined it for the Fuzzy Logic Tool box, 
expects the output membershi p functions to be fuzzy sets. After the 
aggregation process, there is a fuzzy set for each output variable that needs 
defuzzification. It's possible, and in many cases much more efficient, to use a 
si ngl e spi ke as the output membershi p function rather than a di stri buted fuzzy 
set. This is sometimes known as a singleton output membership function, and 
it can bethought of as a pre-defuzzified fuzzy set. It enhances the efficiency of 
the defuzzification process because it greatly simplifies the computation 
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required by the more general Mamdani method, which finds the centroid of a 
two-dimensional function. Ratherthan integrating acrossthe two-dimensional 
function to find the centroid, we use the weighted average of a few data points. 
Sugeno-type systems support this type of model . I n general, Sugeno-type 
systems can be used to model any inference system in which the output 
membership functions are either linear or constant. 


Dinner for Two, Reprise 

I n this section we provide the same two-input one-output three-rule tipping 
problem that you saw in the introduction, only in more detail. The basic 
structure of this example is shown in the diagram below. 

Dinner for two 

a 2 input, 1 output, 3 rule system 



The inputs are crisp 
(non-fuzzy) numbers 
limited to a specific 
range. 


All rules are 
evaluated in parallel 
using fuzzy 
reasoning. 


The results of the rules 
are combined and 
distilled (defuzzified). 


The result is a crisp 
(non-fuzzy) number. 


I nformation flows from left to right, from two inputs to a single output. The 
parallel nature of the rules is one of the more important aspects of fuzzy logic 
systems. I nstead of sharp switching between modes based on breakpoints, we 
will glide smoothly from regions where the system's behavior isdominated by 
either one rule or another. 

I n the Fuzzy Logic Tool box, there are five parts of the fuzzy inference process: 
fuzzification of theinput variables, application of the fuzzy operator (AND or 
OR) in the antecedent, implication from the antecedent to the consequent, 
aggregation of the consequents across the rules, and defuzzification. These 
sometimes cryptic and odd names have very specific meaning that we’ll define 
carefully as we step through each of them in more detail below. 
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Step 1 . Fuzzify Inputs 

The first step is to take the inputs and determine the degree to which they 
belong to each of the appropriate fuzzy sets via membership functions. I n the 
Fuzzy Logic Tool box, the input is always a crisp numerical value limited to the 
universe of discourse of the input variable (in this case the interval between 0 
and 10) and the output is a fuzzy degree of membership in the qualifying 
linguistic set (always the interval between 0 and 1). Fuzzification of the input 
amounts to either a table lookup or a function evaluation. 

The example we're using in this section is built on three rules, and each of the 
rules depends on resolving the in puts into a number of different fuzzy linguistic 
sets: service is poor, service is good, food is rancid, food is delicious, and soon. 
Before the rules can be evaluated, the inputs must befuzzified according to 
each of these linguistic sets. For example, towhat extent is the food really 
delicious? The figure below shows how well the food at our hypothetical 
restaurant (rated on a scale of Oto 10) qualifies, (via its membership function), 
as the linguistic variable "delicious." I n this case, we rated the food as an 8, 
which, given our graphical definition of delicious, corresponds top =0.7 for the 
"delicious" membership function. 


1. Fuzzify 
inputs. 


food is delicious 


food = 8 

input 

(The compliment to the chef would be "your food is delicious to the degree 0. 7.") 
In this manner, each input isfuzzified over all the qualifying membership 
functions required by the rules. 

Step 2. Apply Fuzzy Operator 

Once the inputs have been fuzzified, we know the degree to which each part of 
the antecedent has been satisfied for each rule. If the antecedent of a given rule 
has more than one part, thefuzzy operator is appl ied toobtai n one number that 
represents the result of the antecedent for that rule. This number will then be 
applied to the output function. The input to the fuzzy operator is two or more 


delicious 


/ 


0.7 


Result of 
fuzzification 
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membership values from fuzzified input variables. The output is a single truth 
value. 

As is described in the section on fuzzy logical operations, any number of 
well-defined methods can fill in for the AND operation or the OR operation. I n 
the Fuzzy LogicToolbox, two built-in AND methods are supported: min 
(minimum) and prod (product). Two built-in OR methods are also supported: 
max (maxi mum), and the probabilistic OR method probor. The probabilistic OR 
method (also known as the algebraic sum) is calculated according to the 
equation 

probor{a,b ) =a + b- ab 

I n addition to these built-in methods, you can create your own methods for 
AND and OR by writing any fundi on and setting that to be your method of 
choice. There will be more information on how to do this later. 


Shown below is an example of the OR operator max at work. We're evaluating 
the antecedent of the rule 3 for the tipping calculation. The two different pieces 
of the antecedent (service is excellent and food is delicious) yielded the fuzzy 
membership values 0.0 and 0.7 respedively. The fuzzy OR operator simply 
seleds the maximum of the two values, 0.7, and the fuzzy operation for rule 3 
is complete. I f we were usi ng the probabi I isti c OR method, the resu It wou Id sti 1 1 
be 0.7 in this case. 


l.Fuzzify 2. Apply 

inputs. OR operator (max). 




0.0 

delicious j 

r 






service is excellent or food is delicious 


0.7 

result of 
fuzzy operator 


service = 3 


food = 8 


input 1 


input 2 


Step 3. Apply Implication Method 

Before applying the implication method, we must take care of the rule's weight. 
Every rule has a wdght (a number bdween 0 and 1), which is applied to the 
number given by the antecedent. Generally this weight is 1 (as it is for this 
example) and so it has noeffed at all on the implication process. F rom timeto 
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time you may want to weight one rule relative to the others by changing its 
weight value to something other than 1. 

Once proper weighting has been assigned toeach rule, the implication method 
is implemented. A consequent is a fuzzy set represented by a membership 
function, which weights appropriately the linguistic characteristics that are 
attributed to it. The consequent is reshaped using a function associated with 
the antecedent (a single number). The input for the implication process is a 
single number given by the antecedent, and the output is a fuzzy set. 

I mpli cation is implemented for each rule. Two built-in methods are supported, 
and they are the same functions that are used by the AND method: min 
(minimum), which truncates the output fuzzy set, and prod (product), which 
scales the output fuzzy set. 


Antecedent Consequent 



service = 3 food = 8 


input 1 input 2 

Step 4. Aggregate All Outputs 

Since decisions are based on the testing of all of the rules in an FIS, the rules 
must be combined in some manner in order to make a decision. Aggregation is 
the process by which thefuzzy sets that represent the outputs of each rule are 
combined into a single fuzzy set. Aggregation only occurs once for each output 
variable, just prior to the fifth and final step, defuzzification. The input of the 
aggregation process is the list of truncated output functions returned by the 
implication process for each rule. The output of the aggregation process is one 
fuzzy set for each output variable. 
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Noticethat as long as the aggregation method is commutative (which it always 
should be), then the order in which the rules are executed is unimportant. 
Three built-in methods are supported: max (maximum), probor (probabilistic 
or), and sum (simply the sum of each rule's output set). 

I n the diagram below, all three rules have been placed together to show how 
the out put of each ru I e i s combi ned, or aggregated, i nto a si ngl e fuzzy set whose 
membership function assigns a weighting for every output (tip) value. 


1. Fuzzify inputs. 


2. Apply 
fuzzy 
operation 
(OR = max). 


3. Apply 
implication 
method (min). 


(OR=ma 

• \ poor \ rancid 

v — : ■ V J 



cheap 



1 



0 25% 0 

25% 

1 lf 

service is poor 

or 

food is rancid then 

tip = cheap 

L 


/ \ 


If service is good 


rule 2 has 
no dependency 
on input 2 


1 



JL 


25% 

then tip = average 


excellent/ L ►"l ft >• ••••: K 

* / delicious/ I generous / V / \ 


25% 

If service is excellent or food is delicious then tip = generous 


service = 3 

input 1 


r 

food = 8 

input 2 


125% 4. Apply 

aggregation 
method (max). 


0 „ , 259 

Result of 
aggregation 


Step 5. Defuzzify 

The input for the defuzzification process is a fuzzy set (the aggregate output 
fuzzy set) and the output is a single number. As much as fuzziness helps the 
rule evaluation duringtheintermediatesteps, thefinal desired output for each 
variable is generally a single number. However, the aggregate of a fuzzy set 
encompasses a range of output values, and so must be defuzzified in order to 
resolvea singleoutput value from the set. 
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Perhaps the most popular defuzzification method is the centroid calculation, 
which returns the center of area under the curve. There are five built-in 
methods supported: centroid, bisector, middleof maximum (the average of the 
maximum value of the output set), largest of maximum, and smallest of 
maximum. 



5. Defuzzify the 
aggregate output 
(centroid). 


tip = 1 6.7% 


Result of 
defuzzifica don 


The Fuzzy Inference Diagram 

The fuzzy i nference di agram is the composite of al I the smal ler di agrams we've 
been looking at so far in this section. It simultaneously displays all parts of the 
fuzzy inference process we've examined. I nformation flows through the fuzzy 
inference diagram as shown below. 



Notice how the flow proceeds up from the inputs in the lower left, then across 
each row, or rule, and then down the ruleoutputs to finish in thelower right. 
This is a very compact way of showing everything at once, from linguistic 
variable fuzzification all the way through defuzzification of the aggregate 
output. 
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Shown below is the real full-size fuzzy inference diagram. There's a lot to see 
in a fuzzy i nference diagram, but once you become accustomed to it, you can 
learn a lot about a system very quickly. For instance, from this diagram with 
these particular inputs, you can easily see that the implication method is 
truncation with the min function. The max function is being used for the fuzzy 
OR operation. Rule 3 (the bottom-most row in the diagram shown opposite) is 
having the strongest influence on the output. And so on. The Rule Viewer 
described in 'The Rule Viewer" on page 2-59 is a MATLAB implementation of 
the fuzzy inference diagram. 


1 . Fuzzify inputs. 


2. Apply 
fuzzy 
operation 
(OR = max). 


3. Apply 
implication 
method (min). 


1 . 

A 

poor 


\ rancid 

] 

^^^^cheap 


0 

10 

0 1C 


0% 25% 

If 

service is poor 

or 

food is rancid 

then 

tip = cheap 
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1 


A 


_A_ 

J 

\ 

good X. 


rule 2 has 
no dependency 
on input 2 

] 





25% 


service is good 


tip = average 



0 10 0 10 0% 25% 0% 25% 


If service is excellent or food is delicious then tip = generous 


4. Apply 
aggregation 
method (max). 


service = 3 food = 8 

input 1 input 2 


tip = 16.7% 



5. Defuzzify 
(centroid). 


output 


Customization 

One of the primary goals of the Fuzzy LogicToolbox istohavean open and 
easily modified fuzzy i nference system structure. Thus, the Fuzzy Logic 
Toolbox is designed to give you as much freedom as possible, within the basic 
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constrai nts of the process descri bed here, to customize the fuzzy i nference 
process for your application. For example, you can substitute your own 
MATLAB functions for any of the default functions used in the five steps 
detailed above: you make your own membership functions, AND methods, OR 
methods, implication methods, aggregation methods, and defuzzification 
methods. The next section describes exactly how to build and implement a 
fuzzy i nference system using the tools provided. 
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Building Systems with the Fuzzy Logic Toolbox 

Dinner for Two, from the Top 

Now we're going to work through a similar tipping example, only we’ll be 
building it using the graphical user interface (GUI ) tools provided by the Fuzzy 
LogicToolbox. Although it's possible to use the Fuzzy Logic Tool box by working 
strictly from the command line, in general it's much easier to build a system 
graphically. There are five primary GUI tools for building, editing, and 
observing fuzzy i nference systems in the Fuzzy LogicToolbox: the Fuzzy 
I nference System or FI S Editor, the Membership Function Editor, the Rule 
Editor, the Rule Viewer, and the Surface Viewer. These GUIs are dynamically 
linked, in that changes you make to the FI S using one of them, can affect what 
you see on any of the other open GUIs. You can have any or all of them open for 
any given system. 

I n addition to these five primary GUIs, the tool box includes the graphical 
ANFIS Editor GUI , which is used for building and analyzing Sugeno-type 
adaptive neural fuzzy inference systems. The ANFIS EditorGUI is discussed 
later in this chapter, in the section, "S u gen o-Type Fuzzy I nference" on page 
2 - 86 . 
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The FIS Editor handles the high level issues for the system: How many input 
and output variables? What are their names? The Fuzzy Logic Tool box doesn't 
limit the number of inputs. However, the number of inputs may be limited by 
the available memory of your machine. If the number of inputs istoolarge, or 
the number of membership functions is too big, then it may also be difficult to 
analyze the FIS using the other GUI tools. 

The Membership Function Editor is used to define the shapes of all the 
membership functions associated with each variable. 

The Rule Editor is for editing the list of rules that defines the behavior of the 
system. 
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The Rule Viewer and the Surface Viewer are used for looking at, as opposed to 
editing, the FI S. They are strictly read-only tools. The Rule Viewer is a 
M ATLAB-based display of the fuzzy inference diagram shown at theend of the 
last section. Used as a diagnostic, it can show (for example) which rules are 
active, or how individual membership function shapes are influencing the 
results. The Surface Viewer is used to display the dependency of one of the 
outputs on any one or two of the i nputs— that is, it generates and plots an 
output surface map for the system. 

Thischapter began with an illustration similar to the one below describing the 
main parts of a fuzzy inference system, only the one below shows how the three 
Editors fit together. The two Viewers examine the behavior of the entire 
system. 


The General Case... A Specific Example... 


Input 


Output 


Rules 


/ \ 

Input Output 

terms terms 

(interpret) (assign) 


service 


tip 


4 

if service is poor then tip is cheap 
if service is good then tip is average 
if service is excellent then tip is generous 


/ \ 

service = tip = 

(poor, (cheap, 

good, average, 

excellent} generous} 


The GUI Editors... 


The FIS Editor 



The Rule Editor 


/ \ 

The Membership 
Function Editor 


The five primary GUIs can all interact and exchange information. Any one of 
them can read and write both to the workspace and to the disk (the read-only 
viewers can still exchange plots with the workspace and/or the disk). For any 
fuzzy inference system, any or all of these five GUIs may be open. If more than 
one of these editors is open for a single system, the various GUI windows are 
aware of the existence of the others, and will, if necessary, update related 
windows. Thus if the names of the membership functions are changed using 
the Membership Function Editor, those changes are reflected in the rules 
shown in the Rule Editor. The editors for any number of different FIS systems 
may be open simultaneously. The FI S Editor, the Membership Function 
Editor, and the Rule Editor can all read and modify the FI S data, but the Rule 
Viewer and the Surface Viewer do not modify the FIS data in any way. 
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Getting Started 

We’ll start with a basic description of a two-input, one-output tipping problem 
(based on tipping practices in theU.S.). 

The Basic Tipping Problem. Given a number between 0 and 10 that represents the 
quality of service at a restaurant (where 10 is excellent), and another number 
between 0 and 10 that represents the quality of the food at that restaurant 
(again, 10 is excellent), what should the tip be? 

The starting point is to writedown the three golden rules of tipping, based on 
years of personal experience in restaurants. 

1. If the service is poor or the food is rancid, then tip ischeap. 

2. If the service is good, then tip is average 

3. If the service is excellent or the food is delicious, then tip is generous. 

We’ll assume that an average tip is 15%, a generous tip is 25%, and a cheap tip 
is 5%. It's also useful to have a vague idea of what the tipping function should 
look like. 



Obviously the numbers and the shape of the curve are subject to local 
traditions, cultural bias, and so on, but the three rules are pretty universal. 

Now we know the rules, and we have an idea of what the output should look 
like. Let's begin working with the GUI tools to construct a fuzzy inference 
system for this decision process. 
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The FIS Editor 


These menu items allow you to save, open, 
or edit a fuzzy system using any of the five 
basic GUI tools. 


Double-click on an input 
variable icon to open the 
Membership Function Editor. 


Double-click on the system 
diagram to open the 

Rule Editor. 


This status line describes the 
most recent operation. 


The name of the system is displayed 
here. It can be changed using one of the 
Save as... menu options. 


These pop-up menus are used to 
adjust the fuzzy inference functions, 
such as the defuzzification method. 



FIS Name: 


And method 
Or method 
Implication 
Aggregation 
Defuzzification 


tipper 


(mamdani] 



FIS Type: 


| centroid 


Current Variable 
Name 
Type 
Range 


Help 


System "tipper": 2 inputs, 1 output, and 3 rules 



Double-click on the icon for 
the output variable, tip, to 
open the Membership 
Function Editor. 


This edit field is used to name 
and edit the names of the 
input and output variable. 


The following discussion walks you through building a new fuzzy inference 
system from scratch. If you want to save time and follow along quickly, you can 
load the already built system by typing 

fuzzy tipper 

This will load the FI S associated with the file tipper .f is (the . f is is implied) 
and launch theFIS Editor. However, if you load the pre-bui It system, you will 
not be building rules and constructing membership functions. 

The FI S Editor displays general information about a fuzzy inference system. 
There's a simple diagram at the top that shows the names of each input 
variable on the left, and those of each output variable on the right. The sample 
membership functions shown in the boxes arejust icons and do not depict the 
actual shapes of the membership functions. 

Below the diagram is the name of the system and the type of inference used. 
The default, Mamdani-type inference, is what we've been describing so far and 
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what we’ll continue to use for this example. Another slightly different type of 
inference, called Sugeno-type inference, is also available. This method is 
explained in "Sugeno-Type Fuzzy Inference" on page 2-86. Below the name of 
the fuzzy inference system, on the left side of the figure, are the pop-up menus 
that allow you to modify the various pieces of the inference process. On the 
right side at the bottom of thefigure is the area that displays the name of either 
an input or output variable, its associated membership function type, and its 
range. The latter two fields are specified only after the membership functions 
have been. Below that region arethe Help and Close buttons that call up 
online help and close the window, respectively. At the bottom is a status line 
that relays information about the system. 

T o start this system from scratch, type 

fuzzy 

at the M ATLAB prompt. The generic untitled FIS Editor opens, with one input, 
labeled inputl, and one output, labeled outputl. For this example, we will 
construct a two-i nput, one output system, so go to the Edit menu and select 
Add input. A second yellow box labeled input2 will appear. The two inputs we 
will have in our example are service and food. Our one output is tip. We'd like 
tochange the variable names to reflect that, though. 

1 Click once on the left-hand (yellow) box marked inputl (the box will be 
highlighted in red). 

2 I n the white edit field on the right, change inputl to service and press 

Return. 

3 Click once on the left-hand (yellow) box marked input2 (the box will be 
highlighted in red). 

4 I n the white edit field on the right, change input2 to food and press Return. 

5 Click once on the right-hand (blue) box marked outputl. 

6 I n the white edit field on the right, change outputl to tip. 
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7 From the File menu select Save to workspace as.... 



8 Enter the variable name tipper and click on OK. 

You will see the diagram updated to reflect the new names of the input and 
output variables. There is now a new variable in the workspace cal led tipper 
that contains al I the i nformation about this system. By savi ng to the workspace 
with a new name, you also rename the entire system. Your window will look 
something I ike this. 



Leavethe inference options in the lower left in their default positions for now. 
You've entered all the information you need for this particular GUI. Next 
define the membership functions associated with each of the variables. To do 
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this, open the Membership Function Editor. You can open the Membership 
Function Editor in one of three ways: 

• Pull down the View menu item and select Edit Membership Functions.... 

• Double-click on the icon for the output variable, tip. 

• Type mf edit at the command line. 

The Membership Function Editor 


These menu items allow you 
to save, open, or edit a fuzzy 
system using any of the five 
basic GUI tools. 


These text fields display 
the name and type of , 
the current variable. 


This edit field lets 
you set the range of 
the current variable. 

This edit field lets you set 
the display range of the 
current plot. 


This is the "Variable Palette" 
area. Click on a variable here 
to make it current and edit its 
membership functions. 


This graph field displays all 
the membership functions 
of the current variable. 



input 


This status line describes! 
the most recent operation. \ 


Range 

\ 

Display Range 


[ 010 ] 


[ 010 ] 


Current Membership Function (click on MF to select) 
Name 
Type 
Params 


[15 5 ] 


Click on a line to select it and you 
can change any of its attributes, 
including name, type and 
numerical parameters. Drag your 
mouse to move or change the shape 
of a selected membership function. 


Help 


This edit field lets you 
change the name of the 
current membership 
function. 


This pop-up menu lets 
you change the type 
of the current 
membership function. 


This edit field lets 
you change the 
numerical 
parameters for the 
current membership 
function. 


s Ready 


The Membership Function Editor shares some features with the FIS Editor. In 
fact, all of the five basic GUI tools have similar menu options, status lines, and 
Help and Close buttons. The Membership Function Editor is the tool that lets 
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you display and edit all of the membership functions associated with all of the 
input and output variables for the entire fuzzy i nference system. 

When you open the Membership F unction Editor to work on a fuzzy inference 
system that does not already exist in the workspace, there are not yet any 
membership functions associated with the variables that you have just defined 
with the FIS Editor. 

On the upper left side of the graph area in the Membership F unction Editor is 
a "Variable Palette" that lets you set the membership functions for a given 
variable.Toset upyour membership functions associated with an input or an 
output variablefor the FI S, select an FIS variable in this region by clicking on 
it. 

Next select the Edit pull-down menu, and choose Add MFs.... A new window 
will appear, which allows you to select both the membership function type and 
the number of membership functions associated with the selected variable. I n 
thelower right corner of the window are the controls that let you changethe 
name, type, and parameters (shape), of the membership function, once it has 
been selected. 

The membershi p functions from the current vari able are displ ayed i n the mai n 
graph. These membership functions can be manipulated in two ways. You can 
first use the mousetoselect a particular membership function associated with 
a given variablequality, (such as poor, for the variable, service), and then drag 
the membership function from side to side. This will affect the mathematical 
description of the quality associated with that membership function for a given 
variable. The selected membership function can also be tagged for dilation or 
contraction by clicking on the small square drag points on the membership 
function, and then dragging the function with the mouse toward th eoutsidz 
for dilation, or toward th einside, for contraction. This will changethe 
parameters associated with that membership function. 

Below the Variable Palette is some information about the type and name of the 
current variable. There is a text field in this region that lets you changethe 
limits of the current variable's range (universe of discourse) and another that 
lets you set the limits of the current plot (which has no real effect on the 
system). 

The process of specifying the input membership functions for this two input 
tipper problem is as follows: 
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1 Select the input variable, service, by double-clicking on it. Set both the 
Range and the Display Range to the vector [0 10]. 

2 Select Add MFs... from the Edit menu. The window below pops open. 



3 Use the pul I -down tab to choose gaussmf for MF Typeand 3 for Number of 
MFs. This adds three Gaussian curves to the input variable service 

4 Click once on the curve with the leftmost hump. Change the name of the 
curve to poor. T o adjust the shape of the membership function, either use 
the mouse, as described above, or type in a desired parameter change, and 
then click on the membership function. The default parameter listing for 
this curve is [1.5 0]. 

5 Name the curve with the middle hump, good, and the curve with the 
rightmost hump, excellent. Reset the associated parameters if desired. 

6 Select the input variable, food, by clicking on it. Set both the Range and the 
Display Range to the vector [0 10]. 

7 Select Add MFs... from the Edit menu and add two trapmf curves to the 
input variable food. 

8 Click once directly on the curve with the leftmost trapazoid. Change the 
name of the curve to rancid. To adjust the shape of the membership 
function, either use the mouse, as described above, or type in a desired 
parameter change, and then cl ick on the membershi p function. The default 
parameter listing for this curve is [0 0 13]. 

9 Name the curve with the rightmost trapazoid, delicious, and reset the 
associated parameters if desired 
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Next you need tocreatethe membership functions for the output variable, tip. 
To create the output variable membership functions, usethe Variable Palette 
on the left, selecting the output variable, tip. The inputs ranged from 0 to 10, 
but the output scale is going to be a tip between 5 and 25 percent. 

Use triangular membership function types for the output. First, set the Range 
(and the Display Range) to [0 30], to cover the output range. I nitially, the 
cheap membership function will have the parameters [0 5 10], the average 
membership function will be [10 15 20], and the generous membership 
function will be [20 25 30]. Your system should look something like this. 



Now that the variables have been named, and the membership functions have 
appropri ate shapes and names, you 're ready to write down the rules. To cal I up 
the Rule Editor, gototheView menu and select Edit rules..., or type ruleedit 
at the command I i ne. 
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Tutorial 


The Rule Editor 


Input or output selection menus 


The menu items allow 
you to save, open, or 
edit a fuzzy system 
using any of the five 
basic GUI tools. 


Link input 
statements in rules. 


This status line 
describes the most 
recent operation. 


/ Rule Editor: tipper 


. File Edit View Options 


1 . If (service is poor) or (food is rancidl then (tip is cheap) (1 ) 


2. If (service is good) then (tip is averagef(1) 

3. If (service is excellent) or (food is dejicious) th^n (tip is generous) (1) 




poor 

3 

1 rancid 

good 


delicious 

excellent 

none 

d 

none 

not 

\ 


r not 

1 

r Connection 


weight: 




Negate input or output 
statements in rules. 


Create or edit rules with the GUI buttons and 
choices from the input or output selection menus. 


The rules are 
entered 
automatically 
using the GUI 
tools. 


The Help button 
gives some 
information about 
how the Rule 
Editor works, and 
the Close button 
closes the window. 


Constructing rules using the graphical Rule Editor interface is fairly 
self-evident. Based on the descriptions of the input and output variables 
defined with the FIS Editor, the Rule Editor allows you to construct the rule 
statements automatically, by clicking on and selecting one item in each input 
variable box, one item in each output box, and one connection item. Choosing 
none as one of the variable qualities will exclude that variablefrom a given 
rule. Choosing not under any variable name will negate the associated quality. 
Rules may be changed, deleted, or added, by clicking on the appropriate button. 

The Rule Editor also has some familiar landmarks, similar tothose in the FI S 
Editor and the Membership Function Editor, including the menu bar and the 
status line. The Format pop-up menu is avail able from the Options pull-down 
menu from the top menu bar— this is used to set the format for the display. 
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Similarly, Language can beset from under Options as well. The Help button 
will bring up a MATLAB Help window. 

To insert the first rule in the Rule Editor, select the following: 

• poor under the variable service 

• rancid under the variable food 

• the radio button, or, in the Connection block 

• cheap, under the output variable, tip. 

The resulting rule is: 

1. If (service is poor) or (food is rancid) then (tip is cheap) (1) 

The numbers in the parentheses represent weights that can be applied to each 
rule if desired. You can specify the weights by typing in a desired number 
between zero and one under the Weight: setting. If you do not specify them, the 
weights are assumed to be unity ( 1 ) . 

Follow a similar procedure to insert the second and third rules in the Rule 
Editor to get: 

2. If (service is poor) or (food is rancid) then (tip is cheap) (1) 

2. If (service is good) then (tip is average (1) 

3. If (service is excdlent) or (food is delicious) then (tip is generous) (1) 

Tochangea rule, first click on theruleto be changed. Next make the desired 
changes to that rule, and then click on Change rule. For example, to change 
the first ruleto 

2. If (service not poor) or (food not rancid) then (tip is not cheap) (1) 

click not under each variable, and then click Change rule. 

The Format pop-up menu from the Options menu indicates that you're looking 
at the verbose form of the rules. T ry changing it to symbolic. You will see 

2. (servi ce=poor) =>(tip=cheap) (1) 

2. (servi ce=good) => (ti p=average) (1) 

3. (servi ce=excd lent) => (ti p=generous) (1) 
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There is not much difference in thedisplay really, but it's slightly more 
language neutral, since it doesn’t depend on terms like "if" and "then.” If you 
change the format to indexed, you’ll see an extremely compressed version of the 
rules that has squeezed all the language out. 

1 , 1 ( 1 ) : 1 

2 , 2 ( 1 ) : 1 

3, 3 (1) : 1 

This is the version that the machine deals with. The first column in this 
structure corresponds totheinput variable, the second column corresponds to 
the output variable, thethird column displays the weight applied toeach rule, 
and the fourth column is shorthand that indicates whether this is an OR (2) 
rule or an AND (1) rule. The numbers in thefirst two columns refer to the index 
number of the membership function. A literal interpretation of rule 1 is: "if 
input 1 is M FI (thefirst membership function associated with input 1) then 
output lshould beMFl (thefirst membership function associated with output 
1) with the weight 1." Since there is only one input for this system, the AND 
connective implied by the 1 in the last column is of no consequence. 

The symbolic format doesn’t bother with the terms, if, then, and soon. The 
indexed format doesn’t even bother with the names of your variables. 
Obviously the functionality of your system doesn’t depend on how well you 
have named your variables and membership functions. The whole point of 
naming variables descriptively is, as always, making the system easier for you 
to interpret. Thus, unless you have some special purpose in mind, it will 
probably be easier for you to stick with the verbose format. 

At this point, the fuzzy i nference system has been completely defined, in that 
the variables, membership functions, and the rules necessary to calculate tips 
are in place. It would be nice, at this point, to look at a fuzzy inference diagram 
I ike the one presented at the end of the previous section and verify that 
everything is behaving the way wethink it should. This is exactly the purpose 
of theRuleViewer, the next oftheGUI toolswe’ll look at. From the View menu, 
select View rules.... 
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The Rule Viewer 


The menu items allow 
you to save, open, or edit 
a fuzzy system using any 
of the five basic GUI 
tools. 


Slide this line to change 
your input values, and 
generate a new output 
response. 


This edit field allows you 
to set the input 
explicitly. 


Each column of plots (yellow) shows how 
the input variable is used in the rules. The 
input values are shown here at the top. 


This column of plots (blue) 
shows how the output variable 
is used in the rules. 



/ Rule Viewer: tipper 


''File Edit View Optiory 

serve© =5 ' 




Each row of plots 
represents one rule (here 
there are 3). Click on a 
rule to display it in the 
status bar. 


tip = 1 5 


' Input: | [5 5] 


Plot points: I iQI 


Move: | e ft 


| right | down) 


Opened system tipper, 3 rules 


Help 


This line provides a 
defuzzified value. 


The bottom-right plot 
shows how the output of 
each rule is combined to 
make an aggregate 
output and then 
defuzzified. 


Shift the plots left, right, 
up, or down with these 
buttons. 


This status line describes the most recent operation. 


The Rule Viewer displays a roadmap of the whole fuzzy inference process. It's 
based on thefuzzy i nference diagram descri bed i n the previous section. Y ou see 
a single figure window with 10 small plots nested in it. The three small plots 
across the top of the f i gu re represent the antecedent and consequent of the f i rst 
rule. Each rule is a row of plots, and each column is a variable. The first two 
columns of plots (the six yellow plots) show the membership functions 
referenced by the antecedent, or the if-part of each rule. Thethird column of 
plots (the three blue plots) shows the membership functions referenced by the 
consequent, or the then-part of each rule. If you click once on a rule number, 
the corresponding rule wi 1 1 bedisplayed at the bottom of the figure. Noticethat 
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under food, there is a plot which is blank. This corresponds to the 
characterization of none for the variable food in the second rule. The fourth 
plot in the third column of plots represents the aggregate weighted decision for 
the given i nference system. This decision will depend on the input values for 
the system. 

There are also the now familiar items I ike the status line and the menu bar. I n 
the lower right there is a text field into which you can enter specific input 
values. For the two-in put system, you will enter an input vector, [9 8], for 
example, and then click on input. You can also adjust these input values by 
clicking anywhere on any of the three plots for each input. This will move the 
red index line horizontally, tothe point where you have clicked. You can also 
just click and drag this line in order to change the input values. When you 
release the line, (or after manually specifying the input), a new calculation is 
performed, and you can see the whole fuzzy i nference process take place. 
Where the index line representing service crosses the membership function 
line "service is poor" in the upper left plot will determine the degree to which 
rule one is activated. A yellow patch of color under the actual membership 
function curve is used to make the fuzzy membership value visually apparent. 
Each of the characterizations of each of the variables is specified with respect 
tothe input index line in this manner. If we follow rule 1 across the top of the 
diagram, we can seetheconsequent "tip is cheap" has been truncated toexactly 
the same degree as the (composite) antecedent— this is the implication process 
in action. The aggregation occurs down thethird column, and the resultant 
aggregate plot is shown in the single plot to be found in the lower right corner 
of the plot field. The defuzzified output value is shown by the thick line passing 
through the aggregate fuzzy set. 

The Rule Viewer allows you to interpret the entire fuzzy inference process at 
once. The Rule Viewer also shows how the shape of certain membership 
functions influences the overall result. Since it plots every part of every rule, it 
can become unwieldy for particularly large systems, but, for a relatively small 
number of inputs and outputs, it performs well (depending on how much screen 
space you devote to it) with up to 30 rules and as many as 6 or 7 variables. 

The Rule Viewer shows one calculation at a time and in great detail. In this 
sense, it presents a sort of micro view of the fuzzy inference system. Ifyou want 
to see the entire output surface of your system, that is, the entire span of the 
output set based on the enti re span of the i nput set, you need to open up the 
Surface Viewer. This is the last of our five basic GUI tools in the F uzzy Logic 
Toolbox, and you open it by selecting View surface... from the View menu. 
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The Surface Viewer 



Upon opening the Surface Viewer, we are presented with a two-dimensional 
curve that represents the mapping from service quality to tip amount. Since 
this is a one-input one-output case, we can see the entire mapping in one plot. 
Two-input one-output systems also work well, as they generate 
three-dimensional plots that MATLAB can adeptly manage. When we move 
beyond three dimensions overall, we start to encounter trouble displaying the 
results. Accordingly, the Surface Viewer isequipped with pop-up menus that 
let you select any two inputs and any one output for plotting. J ust below the 
pop-up menus are two text input fields that let you determine how many x-ax/'s 
and y-axis grid lines you want to include. This allows you to keep the 
calculation time reasonablefor complex problems. Pushing the Evaluate 
button initiates the calculation, and the plot comes up soon after the 
calculation is complete. To change the x-ax/'s or y-ax/'s grid after the surface is 
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in view, simply change the appropriate text field, and click on either X-grids: 
or Y-grids:, according to which text field you changed, to redraw the plot. 

The Surface Viewer has a special capability that is very helpful in cases with 
two (or more) inputs and one output: you can actually grab the axes and 
reposition them to get a different three-dimensional view on the data. The Ref. 
Input: field is used in situations when there are more inputs required by the 
system than the surface is mapping. Suppose you have a four-input one-output 
system and would liketoseethe output surface. The Surface Viewer can 
generate a three-di mensi onal output surface where any two of the i nputs vary, 
but two of the inputs must be held constant since computer monitors cannot 
display a five-dimensional shape. I n such a case the input would be a 
four-dimensional vector with NaNs holding the pi ace of the varying inputs while 
numerical values would indicate those values that remain fixed. An NaN is the 
IEEE symbol for "not a number." 

This concludes the quick walk-through of each of the main GUI tools. Notice 
that for the tipping problem, the output of the fuzzy system matches our 
origi nal idea of the shape of the fuzzy mappi ng from service to ti p fai rly wel I . 

I n hindsight, you might say, "Why bother? I could have just drawn a quick 
lookup table and been done an hour ago!" However, if you are interested in 
solving an entire class of similar decision-making problems, fuzzy logic may 
provide an appropriate tool for the solution, given its ease with which a system 
can be quickly modified. 


Importing and Exporting from the GUI Tools 

When you save a fuzzy system to disk, you're saving an ASCI I text F I S file 
representation of that system with the file suffix . f is. This text file can be 
edited and modified and is simple to understand. When you save your fuzzy 
system to the MATLAB workspace, you 're creating a variable (whose name you 
choose) that will act as a MATLAB structure for the FI S system. FI S files and 
F I S structures represent the same system. 


Note: If you do not save your FIS to your disk, but only save it to the 
MATLAB workspace, you wi II not be able to recover it for use in a new 
MATLAB session. 
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Customizing Your Fuzzy System 

If you want to include customized functions as part of your use of the F uzzy 
Logic T ool box, there are a few guidel i nes you need to fol I ow. The AN D method, 
OR method, aggregation method, and defuzzification method functions you 
provide need to work in a similar way to max, min, or prod in MATLAB. That is, 
they must be able to operate down the columns of a matrix. For example, the 
implication method does an element by element matrix operation, similar to 
the min function, as in 

a=[ 1 2 ; 34]; 
b=[2 2; 22]; 
min(a,b) 
ans = 

1 2 

2 2 

Custom Membership Functions 

You can create your own membership functions using an M-file. The values 
these functions can take must be between 0 and 1. There is a limitation on 
customized membership functions in that they cannot use morethan 16 
parameters. 

To define a custom membership function named custmf : 

1 Create an M-file for a function, custmf .m, that takes values between 0 and 
1, and depends on at most 16 parameters. 

2 Choose the Add Custom MF item in the Edit menu on the Membership 
Function Editor GUI. 

3 Enter your custom membership function M-filename, custmf, inthe M-file 
function name text box. 

4 E nter the vector of parameters you want to use to parameterize your 
customized membership function in the text box next to Parameter list. 

5 Give the custom membership function a name different from any other 
membership function name you will use in your FIS. 

6 Select OK. 
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Here is some sample code for a custom membership function, testmf 1 , that 
depends on eight parameters between 0 and 10. 

function out = testmf1(x, params) 
for i=1 : length (x) 
if x(i)<params(1 ) 
y(i)=params(1 ) ; 
elseif x(i)<params(2) 
y(i)=params(2) ; 
elseif x(i)<params(3) 
y(i)=params(3) ; 
elseif x(i)<params(4) 
y(i)=params(4) ; 
elseif x(i)<params(5) 
y(i)=params(5) ; 
elseif x(i)<params(6) 
y(i)=params(6) ; 
elseif x(i)<params(7) 
y(i)=params(7) ; 
elseif x(i)<params(8) 
y(i)=params(8) ; 
else 
y(i)=0; 
end 
end 

out=.1*y' ; 

You can try naming this file testmf 1 .m and loading it into the Membership 
Function Editor using the parameters of your choice. 
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Working from the Command Line 

The tipping example system is one of many example fuzzy i nference systems 
provided with the F uzzy Logic Tool box. The FI S is always cast as a MATLAB 
structure. To load this system (rather than bothering with creating it from 
scratch), type 

a = readfis( 'tipper. fis 1 ) 

MATLAB will respond with 

a = 


name : 

1 tipper 1 

type: 

1 mamdani 1 

andMethod : 

1 min 1 

orMethod : 

1 max 1 

def uzzMethod : 

1 centroid 1 

impMethod : 

1 min ' 

aggMethod : 

1 max 1 

input : 

[1x2 struct 

output : 

[1x1 struct 

rule : 

[1x3 struct 


The labels on the left of this listing represent the various components of the 
MATLAB structure associated with tipper. fis. You can access the various 
components of this structure by typing the component name after typing a. At 
the MATLAB command line, type 

a. type 

for example. MATLAB will respond with 

ans = 
mamdani 

The function 
getfis(a) 

returns almost the same structure information that typing a, alone does. 

getfis(a) returns 

Name = tipper 

Type = mamdani 
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Numlnputs = 2 
InLabels = 

service 

food 

NumOutputs = 1 
OutLabels = 
tip 

NumRules = 3 
AndMethod = min 
OrMethod = max 
ImpMethod = min 
AggMethod = max 
DefuzzMethod = centroid 

Noticethat some of these fields are not part of the structure, a. Thus, you 
cannot get information by typing a. inlabels, but you can get it by typing: 

getfis(a, 1 Inlabels 1 ) 

Similarly, you can obtain structure information using getfis in this manner. 

getfis(a, 1 input 1 , 1 ) 
getfis(a, 'output' ,1 ) 
getfis (a, 1 input 1 , 1 , 1 mf 1 , 1 ) 

The structure, field syntax also generates this information. For more 
information on thesyntaxfor MATLAB structures, seeChapter 13, "Structures 
and Cell Arrays," in Using MATLAB. 

For example, type 

a. input 


or 


a.input(1 ) .mf (1 ) 

Thefunction getfis is loosely modeled on theFlandleGraphics®function get. 
There is also a function called setfis that acts as the reciprocal to getfis. It 
al lows you to change any property of an FIS. F or example, if you wanted to 
change the name of this system, you could type 

a = setf is(a, 1 name 1 , 1 gratuity 1 ) ; 
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However, since a is already a MATLAB structure, you can set this information 
more simply by typing 

a. name = 'gratuity' ; 

N ow the F I S structure a has been changed to reflect the new name. I f you want 
a little more insight into this FIS structure, try 

showfis(a) 

This returns a printout listing all the information about a. This function is 
intended more for debugging than anything else, but it shows all the 
information recorded in the FI S structure 

Si nee the variable, a, designates the fuzzy tipping system, you can call up any 
of the GUIs for the tipping system directly from the command line. Any of the 
following will bring up the tipping system with the associated GUI . 

• f uzzy (a) : brings up the FIS Editor 

• mfedit(a): brings up the Membership Function Editor 

• ruleedit(a) : brings up the Rule Editor 

• ruleview(a) : brings up the Rule Viewer 

• surf view(a) : brings up the Surface Viewer 

If, in addition, a is a Sugeno-type FIS, then anfisedit(a) will bring upthe 
ANFIS Editor GUI. 

Once any of these GUIs has been opened, you can access any of the other GUIs 
using the pull-down menu rather than the command line. 


System Display Functions 

There are three functions designed to give you a high-level view of your fuzzy 
i nference system from the command line: plotf is, plotmf , and gensurf . The 
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first of these displays the whole system as a block diagram much as it would 
appear on the FIS Editor. 

plotfis(a) 



After closing any open MATLAB figures or GUI windows, the function plotmf 
plots all the membership functions associated with a given variable as follows: 

plotmf (a, 1 input 1 , 1 ) 


returns 
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service 


plotmf (a, 1 output 1 , 1 ) 



0 5 10 15 20 25 30 

tip 


These plots will appear in the Membership Function Editor GUI, or in an open 
MATLAB figure, if plotmf is called while either of these is open. 
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Finally, the function gensurf will plot any one or two inputs versus any one 
output of a given system. The result is either a two-dimensional curve, or a 
three-dimensional surface. Note that when there are three or more inputs, 
gensurf must be generated with all but two inputs fixed, as isdescribed in the 
description of genfis in Chapter 3, "Reference" . 

gensurf (a) 



Building a System from Scratch 

It is possible to use the F uzzy Logic Tool box without bothering with the GUI 
tools at all. For instance, to build the tipping system entirely from the 
command line, you would use the commands newfis, addvar, addmf, and 
addrule. 

Probably the trickiest part of this process is learning the shorthand that the 
fuzzy i nference systems use for building rules. This is accomplished using the 
command line function, addrule. 

Each variable, input, or output, has an index number, and each membership 
function has an index number. The rules are built from statements like this 

if inputl is M FI or input2 is MF3then outputl is MF2 (weight =0.5) 

This ruleisturned into a structure according to the foil owing logic: If there are 
m inputs to a system and n outputs, then the first m vector entries of the rule 
structure correspond to inputs 1 through m. The entry in column 1 is the index 
number for the membership function associated with input 1. The entry in 
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column 2 is the index number for the membership function associated with 
input 2. And so on. The next n columns work the same way for the outputs. 
Column m +n + 1 is the weight associated with that rule (typically 1) and 
column m +n +2 specifies the connective used (where AND =1 and OR =2) 
The structure associated with the rule shown above is 

1 3 2 0.5 2 

Here is one way you can build the entire tipping system from the command 
line, using the MATLAB structure syntax. 

a=newfis( 'tipper' ) ; 
a. input (1 ) . name= ' service ' ; 
a.input(1 ) .range=[0 10]; 
a. input (1 ) .mf (1 ) .name='poor' ; 
a. input (1 ) .mf (1 ) . type= ' gaussmf 1 ; 
a.input(1 ) .mf (1 ) .params=[1 .5 0]; 
a. input (1 ) .mf (2) . name= ' good ' ; 
a. input (1 ) .mf (2) . type= ' gaussmf 1 ; 
a.input(1 ) .mf (2) .params=[1 .5 5]; 
a. input (1 ) .mf (3) . name= ' excellent 1 ; 
a. input (1 ) .mf (3) . type= ' gaussmf 1 ; 
a. input (1 ) .mf (3) .params=[1 .5 10] ; 
a. input (2) .name='food ' ; 
a.input(2) .range=[0 10]; 
a . input (2) . mf ( 1 ) . name= ' rancid ' ; 
a. input (2) . mf (1 ) . type= ' trapmf ' ; 
a.input(2) .mf (1 ) .params=[ -2 0 1 3]; 
a.input(2) .mf (2) . name= ' delicious 1 ; 
a.input(2) .mf (2) . type= ' trapmf ' ; 
a.input(2) .mf (2) .params=[7 9 10 12]; 
a. out put (1 ) . name='tip' ; 
a.output(1 ) .range=[0 30]; 
a. out put (1 ) .mf (1 ) . name= 1 cheap ' 
a. out put (1 ) .mf (1 ) .type='trimf ' ; 
a.output(1 ) .mf (1 ) .params=[0 5 10]; 
a. out put (1 ) .mf (2) . name= 1 average ' ; 
a. out put (1 ) .mf (2) . type= 1 trimf ' ; 
a.output(1 ) .mf (2) .params=[10 15 20]; 
a. out put (1 ) .mf (3) . name= 1 generous 1 ; 
a. out put (1 ) .mf (3) . type= 1 trimf ' ; 
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a.output(1 ) .mf (3) .params=[20 25 30]; 

a.rule(1 ) . antecedent= [ 1 1]; 

a. rule(1 ) . consequent= [ 1 ] ; 

a. rule(1 ) .weight=1 ; 

a. rule(1 ) . connect ion=2 ; 

a. rule(2) .antecedent=[2 0]; 

a. rule (2) .consequent^ 2] ; 

a. rule(2) .weight=1 ; 

a. rule (2) . connect ion=1 ; 

a. rule(3) .antecedent=[3 2]; 

a. rule (3) .consequent^ 3] ; 

a. rule(3) .weight=1 ; 

a. rule (3) . connect ion=2 

Alternatively, here is how you can build the entire tipping system from the 

command line using F uzzy Logic Tool box commands. 

a=newfis( 'tipper' ) ; 

a=addmf (a, 1 input 1 , 1 , 1 service 1 , [0 10] ) ; 

a=addmf (a, 1 input 1 , 1 , 1 poor 1 , 1 gaussmf ',[1.5 0] ) ; 

a=addmf (a, 1 input 1 , 1 , 1 good 1 , 1 gaussmf ',[1.5 5] ) ; 

a=addmf (a, 1 input 1 , 1 , 1 excellent 1 , 1 gaussmf ',[1.5 10] ) ; 

a=addvar(a, 1 input 1 , 'food ' , [0 10] ) ; 

a=addmf (a, ' input ' ,2, ' rancid ' , ' trapmf ' , [ -2 0 1 3]); 

a=addmf (a, ' input ' ,2, ' delicious ',' trapmf 1 , [7 9 10 12]); 

a=addvar (a, 1 output ' , ' tip ' , [0 30] ) ; 

a=addmf (a, ' output ', 1 ,' cheap ',' trimf 1 , [0 5 10]); 

a=addmf (a, ' output ', 1 ,' average ',' trimf ', [10 15 20]); 

a=addmf (a, ' output ', 1 ,' generous ',' trimf ', [20 25 30]); 

ruleList=[ . . . 

11112 
2 0 2 1 1 
32312]; 
a=addrule (a, rule List) ; 
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FIS Evaluation 

To evaluate the output of a fuzzy system for a given input, use the function 
evalfis. For example, the following script evaluates tipper at the input, [1 2], 

a = readfis( 'tipper' ) ; 
evalfis([1 2], a) 
ans = 

5.5586 

This function can also be used for multi pie col lections of inputs, si nee different 
input vectors are represented in different parts of the input structure. By doing 
multiple evaluations at once, you get a tremendous boost in speed. 

evalfis( [35; 27], a) 
ans = 

12.2184 

7.7885 

The FIS Structure 

The FI S structure is the M ATLAB object that contains all the fuzzy inference 
system information. This structure is stored inside each GUI tool. Access 
functions such as getf is and sett is make it easy to exami ne this structure. 
You can also access the FIS structure information using the structure, field 
syntax (see the section, "Working from the Command Line" on page 2-65). 

All the information for a given fuzzy inference system is contained in the FI S 
structure, including variable names, membership function definitions, and so 
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on. This structure can itself bethought of as a hierarchy of structures, as shown 
in the diagram below: 



You can generate a listing of information on the FI S using the showf is 
command, as shown below. 


showfis(a) 


1 . 

Name 

tipper 

2. 

Type 

mamdani 

3. 

Inputs/Outputs 

[21] 

4. 

NumlnputMFs 

[32] 

5. 

NumOutputMFs 

3 

6. 

NumRules 

3 

7. 

AndMethod 

min 

8. 

OrMethod 

max 

9. 

ImpMethod 

min 

10. 

AggMethod 

max 

11 . 

Def uzzMethod 

centroid 

12. 

InLabels 

service 

13. 


food 

14. 

OutLabels 

tip 
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15. 

InRange 

[ 0 10 ] 

16. 


[ 0 10 ] 

17. 

OutRange 

[ 0 30 ] 

18. 

InMFLabels 

poor 

19. 


good 

20. 


excellent 

21 . 


rancid 

22. 


delicious 

23. 

OutMFLabels 

cheap 

24. 


average 

25. 


generous 

26. 

InMFTypes 

gaussmf 

27. 


gaussmf 

28. 


gaussmf 

29. 


trapmf 

30. 


trapmf 

31 . 

OutMFTypes 

trimf 

32. 


trimf 

33. 


trimf 

34. 

InMFParams 

[ 1.5000 

35. 


[ 1.5500 

36. 


[1.51000 

37. 


[0013] 

38. 


[ 7 9 10 10 

39. 

OutMFParams 

[ 0 5 10 0 ] 

40. 


[ 10 15 20 0 

41 . 


[ 20 25 30 0 

42. 

Rule Antecedent 

[11] 

43. 


[20] 

44. 


[32] 

42. 

Rule Consequent 

1 

43. 


2 

44. 


3 

42. 

Rule Weigth 

1 

43. 


1 

44. 


1 

42. 

Rule Connection 

2 

43. 


1 

44. 


2 
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The list of command line functions associated with FI S construction include 
getfis, setfis, showfis, addvar, addmf, addrule, rmvar, and rmmf. 

Saving FIS Files on Disk 

A specialized text file format is used for saving fuzzy inference systems to a 
disk. The functions readfis and writefis are used for reading and writing 
these files. 

I f you prefer, you can modify the F I S by editi ng its . f is text fi le rather than 
using any of theGUIs. You should be aware, however, that changing oneentry 
may oblige you to change another. For example, if you deletea membership 
function using this method, you also need to make certain that any rules 
requiring this membership function are also deleted. 

The rules appear in "indexed" format in a .fis text file. Here is the file 
tipper .fis. 

[System] 

Name= 1 tipper 1 
Type= 1 mamdani ' 

Numlnputs=2 
Num0utputs=1 
NumRules=3 
AndMethod= 1 min ' 

OrMethod= ' max ' 

ImpMethod= 1 min 1 
AggMethod= 1 max 1 
Def uzzMethod= ' centroid ' 

[ Inputl ] 

Name= ' service ' 

Range=[0 10] 

NumMFs=3 

MF1 = 1 poor 1 : 1 gaussmf 1 ,[ 1 .5 0] 

MF2= 1 good 1 : 1 gaussmf 1 , [ 1 .5 5] 

MF3= ' excellent 1 : 1 gaussmf ',[1.5 10] 

[ Input2] 

Name= 1 food 1 
Range=[0 10] 

NumMFs=2 
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MF1 = 1 rancid' : 'trapmf' , [0 0 1 3] 

MF2= ' delicious 1 : 1 trapmf [7 9 10 10] 

[Outputl ] 

Name= 1 tip 1 
Range=[0 30] 

NumMFs=3 

MF1 =' cheap 1 trimf 1 , [0 5 10] 

MF2= ' average trimf 1 ,[ 1 0 15 20] 

MF3= 1 generous 1 : 1 trimf 1 , [20 25 30] 

[Rules] 

1 1 , 1 ( 1 ) : 2 
2 0 , 2 ( 1 ) : 1 
3 2, 3 (1) : 2 
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Working with Simulink 

The Fuzzy LogicToolbox is designed to work seamlessly with Simulink, the 
simulation software availablefromTheMathWorks. Once you've created your 
fuzzy system using the GUI tools or some other method, you're ready to embed 
your system directly into a simulation. 


An Example: Water Level Control 

Picture a tank with a pipe flowing in and a pipe flowing out. You can change 
the valve controlling the water that flows in, but the outflow rate depends on 
the diameter of the outflow pipe (which is constant) and the pressure in the 
tank (which varies with the water level). The system has some very nonlinear 
characteristics. 



A controller for the water level in the tank needs to know the current water 
level and it needs to be able to set the valve. Our controller's input will be the 
water level error (desired water level minus actual water level) and its output 
will be the rate at which the valve is opening or closing. A first pass at writing 
a fuzzy controller for this system might be the following. 

1. If (level is okay) then (valve is no_change) (1) 

2. If (level is low) then (valve is openjast) (1) 

3. If (I e/el is high) then (valve i s dosejast) (1) 

One of the great advantages of the Fuzzy LogicToolbox is the ability to take 
fuzzy systems directly into Simulink and test them out in a simulation 
environment. A Simulink block diagram for this system is shown below. It 
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contains a Simulink block called the F uzzy Logic Controller block. The 
Simulink block diagram for this system is sltank. Typing 

sltank 

at the command line, causes the system to appear. 



At the same time, the file tank.f is is loaded into the FI S structure tank. 
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Some experimentation shows that three rules are not sufficient, sincethe 
water level tendstooscillatearoundthe desired level. This isseen from the plot 
below. 



0 l 1 1 1 1 1 1 1 1 1 1 

0 10 20 30 40 50 60 70 80 90 100 

Time (second) 


We need to add another input, the water level's rate of change, to slow down 
the valve movement when we get close to the right level . 

4. if (ievd is good) and (rate is negative) then (valveis close_slow) (1) 

5. If(l&/el is good) and (rate is positive) then (valveis open _slow) (1) 

The demo, sltank is built with these five rules. You can examine With all five 
rules in operations, the step response by simulating this system. This is done 
by clicking on Start from the pul I -down menu under Simulate, and clicking on 
the Comparison block. The result looks like this 



One interesting feature of the water tank system is that the tank empties much 
moreslowly than it fills up because of the specific value of the outflow diameter 
pipe. Wecan deal with this by setting theclose_slow valve membership 
function to be slightly different from theopen_slow setting. A PI D controller 
does not have this capability. The valve command versus the water level 
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change rate (depicted as water) and the relative water level change (depicted 
as la/el) surface looks I ike this. If you look closely, you can see a slight 
asymmetry to the plot. 



Because the MATLAB technical computing environment supports so many 
tools (I ike the Control System Toolbox, theNeural Network Toolbox, the 
Nonlinear Control Design Blockset, and soon), you can, for example, easily 
make a comparison of a fuzzy controller versus a linear controller or a neural 
network controller. 

For a demonstration of how the Rule Viewer can be used to interact with a 
Fuzzy Logic Controller block in a Simulink model, type 

sltankrule 

This demo contains a block cal led the Fuzzy Controller With Rule Viewer block. 

I n this demo, the Rule Viewer opens when you start the Simulink simulation. 
This Rule Viewer provides an animation of how the rules are fired during the 
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water tank simulation. The windows that open when you simulatethe 
sltankrule demo are depicted as follows: 


File Edit View Simulation Fotmat Tools 

"oj I I I JJ ±J 



The Rule Viewer that opens during the simulation can be used to access the 
Membership Function Editor, the Rule Editor, or any of the other GUIs, (see 
'The Membership Function Editor"on page2-52, or 'TheRuleEditor"on page 
2-56, for more information). 
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For example, you may want to open the Rule Editor to change one of your rules. 
To do so, select the Edit rules menu item under the View menu of the open 
Rule Viewer. Now you can view or edit the rules for this Simulink model: 



I 


I 


I 


r 





It's best if you stop the simulation prior to selecting any of these editors to 
change your FIS. Remember to save any changes you make to your FIS to the 
workspace before you restart the simulation. 

Building Your Own Fuzzy Simulink Models 

To build your own Simulink systems that use fuzzy logic, simply copy the Fuzzy 
Logic Controller block out of sltank (or any of the other Simulink demo 
systems avail able with the tool box) and pi ace it in your own block diagram. You 
can also open the Simulink library called fuzblock, which contains the Fuzzy 
Logic Controller block, the Fuzzy Controller With RuleViewer block, and 
several demo blocks. To access these blocks, type 

fuzblock 
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at the M ATLAB prompt. The following library appears: 



When you use these blocks, make sure that the fuzzy i nference system (FIS) 
structure corresponding to your fuzzy system is both in the M ATLAB 
workspace, and referred toby name in the dialog box associated with the Fuzzy 
Logic Controller block. 

Double-click on the Fuzzy Controller With Rule Viewer block, and the following 
appears: 



This block uses the zero-order hold method for sampling. 

The Fuzzy Logic Controller block is a masked Simulink block based on the 
S-function sf f is . mex. This function is itself based on the same algorithms as 
the function evalf is, but it has been tailored to work optimally within the 
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Simulink environment. For more descriptions of these, seefuzblock on page 
3-27, and sff is on page 3-65. 
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Sugeno-Type Fuzzy Inference 

The fuzzy i nference process we've been referri ng to so far is known as 
Mamdani's fuzzy inference method. It's the most commonly seen fuzzy 
methodology. In this section we discuss the so-called Sugeno, or 
Takagi-Sugeno-Kang method of fuzzy i nference fi rst introduced in 1985 
[Sug85], It issimilartotheMamdani method in many respects. I n fact thefi rst 
two parts of the fuzzy inference process, fuzzifying the inputs and applying the 
fuzzy operator, are exactly the same. The main difference between 
M amdani-type of fuzzy i nference and Sugeno-type is that the output 
membership functions are only linear or constant for Sugeno-type fuzzy 
i nference. 

A typical fuzzy rule in a zero-order Sugeno fuzzy model has the form 
if x is A and y is B then z-k 

whereA and B are fuzzy sets in the antecedent, while k is a crisply defined 
constant in the consequent. When the out put of each rule is a constant I ike this, 
the similarity with Mamdani's method is striking. The only distinctions are the 
fact that all output membership functions aresingleton spikes, and the 
implication and aggregation methods are fixed and can not be edited. The 
implication method is simply multiplication, and the aggregation operator just 
includes all of the singletons. 
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1. Fuzzify inputs 


2. Apply 
fuzzy 




3. Apply 
implication 
method (min). 


1 (OR = max) 

1 


cheap 


J 






If service is poor or food is rancid 


tip = cheap 


2 . 




good 



average 


rule 2 has 
no dependency 
on input 2 


If service is good 


then tip = average 



delic 


ft] 


generous 


If service is excellent or food is delicious then tip = generous 


service = 3 

input 1 


food = 8 

input 2 


output 

tip = 16.3% 


□ 



1 

j 1 

, 





□ 




4. Apply 
aggregation 
method (max). 


5. Defuzzify 

(weighted 

average) 


The figure above shows the fuzzy tipping model developed in previous sections 
of this manual adapted for use as a zero-order Sugeno system. Fortunately it is 
frequently the case that singleton output functions are completely sufficient for 
a given problem's needs. As an example, the system tippersg.fis is the 
Sugeno-type representation of the now-familiar tipping model. If you load the 
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system and plot its output surface, you will see it is almost the same as the 
Mamdani system we've been looking at. 

a = readfis( 'tippersg 1 ) ; 
gensurf (a) 



The more general first-order Sugeno fuzzy model has rules of the form 

if x is A and y is B then z - p*x + q*y + r 

where/4 and B are fuzzy sets in the antecedent, while p, q, and r are all 
constants. The easi est way to vi su al i ze the f i rst -order system i s to thi nk of each 
rule as defining the location of a "moving singleton." That is, the singleton 
output spikes can move around in a linear fashion in the output space, 
depending on what the input is. This also tends to makethe system notation 
very compact and efficient. H igher order Sugeno fuzzy models are possible, but 
they introduce significant complexity with little obvious merit. Sugeno fuzzy 
models whose output membershi p functi ons are greater than fi rst order are not 
supported by the Fuzzy Logic Tool box. 

Because of the linear dependence of each rule on the system's input variables, 
the Sugeno method is ideal for acting as an interpolating supervisor of multiple 
linear controllers that areto be applied, respectively, todifferent operating 
conditions of a dynamic nonlinear system. For example, the performance of an 
aircraft may change dramatically with altitude and Mach number. Linear 
controllers, though easy to compute and well-suited to any given flight 
condition, must be updated regularly and smoothly to keep up with the 
changing state of the flight vehicle. A Sugeno fuzzy i nference system is 
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extremely well suited tothetask of smoothly interpolating the linear gains 
that would be applied across the input space; it's a natural and efficient gain 
scheduler. Similarly, a Sugeno system is suited for modeling nonlinear systems 
by interpolating multiple linear models. 


An Example: Two Lines 

T o see a specific exam pi e of a system with I i near output membershi p functi ons, 
consider the one input one output system stored in sugenol .fis. 

fismat = readfis( 1 sugenol 1 ) ; 
getf is (f ismat , 1 output 1 , 1 ) 

Name = output 
NumMFs = 2 
MFLabels = 
linel 
line2 

Range =[01] 

The output variable has two membership functions: 

getf is (fismat , 1 output 1 , 1 , 1 mf 1 , 1 ) 

Name = linel 
Type = linear 
Params = 

-1 -1 

getf is (fismat , 1 output 1 , 1 , 1 mf 1 ,2) 

Name = line2 
Type = linear 
Params = 

1 -1 

F urther, these membershi p functi ons are I i near functi ons of the i nput vari able. 
The membership function linel is defined by the equation 

output = (-1 )*input + (-1) 

and the membership function line2 is defined by the equation 
output = (1 )*input + (-1) 
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The input membership functions and rules define which of these output 
functions will be expressed and when. 

showrule(fismat) 
ans = 

1. If (input is low) then (output is linel ) (1) 

2. If (input is high) then (output is line2) (1) 

The function plotmf shows us that the membership function low generally 
refers to input values less than zero, while high refers to values greater than 
zero. Thefunction gensurf shows how the overall fuzzy system output switches 
smoothly from the line called linel to the line called line2. 

subplot (2, 1 , 1 ) , plotmf (fismat, 1 input 1 , 1 ) 
subplot(2, 1 ,2) , gensurf (fismat) 



input 



input 


This is just one example of how a Sugeno-type system gives you the freedom to 
incorporate linear systems into your fuzzy systems. By extension, you could 
build a fuzzy system that switches between several optimal linear controllers 
as a highly nonlinear system moves around in its operating space. 

Conclusion 

Because it is a more compact and computati on al I y efficient representation than 
a M amdani system, the Sugeno system lends itself to the use of adaptive 
techniques for constructing fuzzy models. These adaptive techniques can be 
used to customize the membership functions so that the fuzzy system best 
models the data. 
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Note on FIS Conversion: TheMATLAB command line function mam2sug can 
be used to convert a Mamdani system into a Sugeno system (not necessarily 
with a single output) with constant output membership functions. It uses the 
centroid associated with all of the output membership functions of the 
Mamdani system. See Chapter 3 for details. 


Here are some final considerations about the two different methods: 

Advantages of the Sugeno method 

• It's computationally efficient. 

• It works well with linear techniques (e.g., PI D control). 

• It works well with optimization and adaptive techniques. 

• It has guaranteed continuity of the output surface. 

• It's well -suited to mathematical analysis. 

Advantages of the Mamdani method 

• It's intuitive. 

• It has widespread acceptance. 

• It's well-suited to human input. 
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anfis and the ANFIS Editor GUI 

The basic structure of the type of fuzzy i nference system that we’ve seen thus 
far is a model that maps input characteristics to input membership functions, 
input membership function to rules, rules to a set of output characteristics, 
output characteristics to output membership functions, and the output 
membership function to a single-valued output or a decision associated with 
the output. We have only considered membership functions that have been 
fixed, and somewhat arbitrarily chosen. Also, we’ve only applied fuzzy 
inference to modeling systems whose rulestructure is essentially 
predetermined by the user's interpretation of the characteristics of the 
variables in the model. 

I n this section we discuss the use of the function anfis and the AN FI S Editor 
GUI in the F uzzy Logic Tool box. These tools apply fuzzy inference techniques 
to data modeling. As you have seen from the other fuzzy inference GUIs, the 
shape of the membership functions depends on parameters, and changing 
these parameters wi 1 1 change the shape of the membershi p fu ncti on . I nstead of 
just looking atthedatatochoosethemembershipfunction parameters, we will 
see how membership function parameters can be chosen automatically using 
these Fuzzy Logic Tool box applications. 


A Modeling Scenario 

Suppose you want to apply fuzzy inference to a system for which you already 
have a collection of input/output data that you would like to use for modeling, 
model -following, or some similar scenario. You don't necessarily have a 
predetermined model structure based on characteristics of variables in your 
system. 

Therewill be some modeling situations in which you can't just look at the data 
and discern what the membership functions should look like. Rather than 
choosing the parameters associated with a given membership function 
arbitrarily, these parameters could be chosen so as to tailor the membership 
functions to the input/output data in order to account for these types of 
variations in the data values. This is where the so-called neuro-adaptive 
learning techniques incorporated into anfis in theFuzzy Logic Tool box can 
help. 
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Model Learning and Inference Through ANFIS 

The basic idea behind these neuro-adaptive learning techniques is very simple. 
These techniques provide a method for the fuzzy modeling procedure to learn 
information about a data set, in order to compute the membership function 
parameters that best allow the associated fuzzy i nference system to track the 
given input/output data. This learning method works similarly tothat of 
neural networks. The Fuzzy Logic Tool box function that accomplishes this 
membership function parameter adjustment is called anf is. anf is can be 
accessed either from the command line, or through the AN FI S Editor GUI . 

Si nee the functionality of the command line function anf is and the AN FI S 
Editor GUI is similar, they are used somewhat interchangeably in this 
discussion, until we distinguish them through the description of the GUI . 

What Is ANFIS? 

The acronym ANFIS derives its name from adaptive neuro-fuzzy inference 
system. Using a given input/output data set, the toolbox function anf is 
constructs a fuzzy i nference system (FIS) whose membership function 
parameters are tuned (adjusted) using either a back propagation algorithm 
alone, or in combination with a least squares type of method. This allows your 
fuzzy systems to learn from the data they are modeling. 

FIS Structure and Parameter Adjustment 

A network-type structure similar tothat of a neural network, which maps 
inputs through input membership functions and associated parameters, and 
then through output membership functions and associated parameters to 
outputs, can be used to interpret the input/output map. 

The parameters associated with the membership functions will change 
through the learning process. The computation of these parameters (or their 
adjustment) is facilitated by a gradient vector, which provides a measure of 
how well the fuzzy i nference system is modeling the input/output data for a 
given set of parameters. Once the gradient vector is obtained, any of several 
optimization routines could be applied in order to adjust the parameters so as 
to reduce some error measure (usually defined by the sum of the squared 
difference between actual and desired outputs), anf is uses either back 
propagation or a combination of least squares estimation and backpropagation 
for membership function parameter estimation. 
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Familiarity Breeds Validation: Know Your Data 

The modeling approach used by anf is is similar to many system identification 
techniques. First, you hypothesize a parameterized model structure (relating 
inputs to membership functions to rules to outputs to membership functions, 
and soon). Next, you collect input/output data in a form that will be usable by 
anf is for training. You can then use anf is to train the FI S model to emulate 
the training data presented to it by modifying the membership function 
parameters according to a chosen error criterion. 

I n general, this type of modeling works well if the training data presented to 
anfis for training (estimating) membership function parameters isfully 
representative of the features of the data that the trained FI S is intended to 
model. This is not always the case, however. In some cases, data is collected 
using noisy measurements, and thetraining data cannot be representative of 
all the features of the data that will be presented tothe model. This is where 
model validation comes into play. 

Model Validation Using Checking and Testing Data Sets 

Model validation is the process by which the input vectors from input/output 
data sets on which the FIS was not trained, are presented to the trained FIS 
model, to see how well the F I S model predicts the corresponding data set 
output values. This is accomplished with the AN FI S Editor GUI using the 
so-called testing data set, and its use is described in a subsection that follows. 
You can also use another type of data set for model validation in anfis. This 
other type of val idation data set is referred to as the checking data sd: and this 
set is used to control the potential for the model overfitting the data. When 
checking data is presented to anfis as well as training data, the FIS model is 
sel ected to have parameters associ ated with the mi ni mum checki ng data model 
error. 

One problem with model validation for models constructed using adaptive 
techniques is selecting a data set that is both representative of the data the 
trained model is intended to emulate, yet sufficiently distinct from the training 
data set so as not to render the validation process trivial. Ifyou have collected 
a large amount of data, hopefully this data contains all the necessary 
representative features, so the process of selecting a data set for checking or 
testing purposes is made easier. Flowever, ifyou expect to be presenting noisy 
measurements to your model, it's possible the training data set does not 
include all of the representative features you want to model. 
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The basic idea behind using a checking data set for model validation is that 
after a certai n poi nt i n the trai ni ng, the model begi ns overfitti ng the trai ni ng 
data set. I n principle, the model error for the checking data set tends to 
decrease as the training takes pi ace up to the point that overfitting begins, and 
then the model error for the checking data suddenly increases. In the first 
example in the following section, two similar data sets are used for checking 
and training, but the checking data set is corrupted by a small amount of noise. 
This example illustrates of the use of theANFIS Editor GUI with checking 
data to reduce the effect of model overfitting. I n the second example, a training 
data set that is presented to anf is is sufficiently different than the applied 
checki ng data set. By exami ni ng the checki ng error sequence over the trai ni ng 
period, it is clear that the checking data set is not good for model validation 
purposes. This example illustrates the use of the AN FI S Editor GUI to compare 
data sets. 

Some Constraints of anfis 

anf is is much morecomplexthan the fuzzy inference systems discussed so far, 
and is not availablefor all of the fuzzy i nference system options. Specifically, 
anfis only supports Sugeno-type systems, and these must be: 

• F irst or zeroth order Sugeno-type systems 

• Single output, obtained using weighted average defuzzification (linear or 
constant output membership functions) 

• Of unity weight for each rule 

An error occurs if your FIS structure does not comply with these constrai nts. 

Moreover, anfis cannot accept all the customization options that basicfuzzy 
inference allows. That is, you cannot make your own membership functions 
and defuzzification functions; you’ll have to use the ones provided. 

The ANFIS Editor GUI 

To get started with theANFIS Editor GUI, type 
anf isedit 

The following GUI will appear on your screen. 


2-95 



Tutorial 



From this GUI you can 


• Load data (training, testing, and checking) by selecting appropriate radio 
buttons in the Load data portion of the GUI and then selecting Load Data... 

The loaded data is plotted on the plot region. 

• Generate an initial FIS model or load an initial FIS model using the options 
in the Generate FIS portion of the GUI 

• ViewtheFIS model structure once an initial FIS has been generated or 
loaded by selecting the Structure button 
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• ChoosetheFIS model parameter optimization method: backpropagation or a 
mixture of backpropagation and least squares (hybrid method) 

• Choose the number of training epochs and thetraining error tolerance 

• Train the FI S model by selecting theTrain Now button 

This training adjusts the membership function parameters and plots the 
training (and/or checking data) error plot(s) in the plot region. 

• View the FIS model output versus thetraining, checking, or testing data 
output by selecting the Test Now button 

Thisfunction plots thetest data against the FIS output in the plot region. 

You can also use the AN FI S Editor GUI menu bar to load an FI S training 
initialization, save your trained FIS, open a new Sugeno system, or open any 
of the other G U I s to i nterpret thetrainedFIS model . 

Data Formalities and the ANFIS Editor GUI: Checking and Training 

To start training an FI S using either anf is or the ANFIS Editor GUI , first you 
need to have a training data set that contains desired input/output data pairs 
of the target system to be modeled. Sometimes you also want to have the 
optional testing data set that can check the generalization capability of the 
resulting fuzzy i nference system, and/or a checking data set that helps with 
model overfitting during thetraining. The use of a testing data set and a 
checking data set for model validation isdiscussed in "Model Validation Using 
Checking and Testing Data Sets" on page 2-94. As we mentioned previously, 
overfitting is accounted for by testing the FIS trained on thetraining data 
agai nst the checki ng data, and choosi ng the membershi p function parameters 
to be those associated with the mini mum checking error if these errors indicate 
model overfitting. You will have to examine your training error plots fairly 
closely in order to determine this. These issues are discussed later in an 
example. Usually these training and checking data sets are collected based on 
observations of the target system and are then stored in separate files. 


Note on Data Format: Any data set you load into the ANFIS Editor GUI, (or 
that is applied to the command line function anf is) must be a matrix with the 
input data arranged as vectors in all but the last column. The output data 
must be in the last column. 
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Tutorial 


ANFIS Editor GUI Example 1: Checking Data Helps 
Model Validation 

I n this section we look at an examplethat loads similar training and checking 
data sets, only the checking data set is corrupted by noise. 

Loading Data 

To work both of the following examples, you load the training data sets 
(fuzexltrnData and fuzex2trnData) and the checking data sets 
(fuzexl chkData and fuzex2chkData), into the AN FIS Editor GUI from the 
workspace. You may also substitute your own data sets. 

To load these data sets from the directory fuzzydemos into the MATLAB 
workspace, type 

load fuzex1trnData.dat 
load fuzex2trnData.dat 
load fuzex1chkData.dat 
load fuzex2chkData.dat 

from the command I i ne. 


Note on loading data: You may also want to load your data set from the 
fuzzydemos or any other directory on thedisk, using theANFIS Editor GUI, 
directly. 


Open the AN FI S Editor GUI by typing anf isedit. To load the training data 
set, click on Training, worksp. and then Load Data..., 

The small GUI window that pops up allows you to type in a variable name from 
the workspace. Type in fuzexltrnData, as shown below. 
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anfis and the ANFIS Editor GUI 


The training data appears in the plot in the center of the GUI as a set of circles. 



Notice the horizontal axis is marked data set index. This index indicates the 
row from which that input data value was obtained (whether or not the input 
is a vector or a scalar). Next click on Checking in theType column of the Load 
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data portion of the GUI to load f uzexl chkData from the workspace. This data 
appears in theGUI plot as p/usses superimposed on the training data. 



+++ Checking data 
ooo Training data 


This data set will be used to train a fuzzy system by adjusting the membership 
function parameters that best model this data. The next step is to specify an 
initial fuzzy i nference system for anfis to train. 


Initializing and Generating Your FIS 

You can either initialize the FIS parameters to your own preference, or if you 
do not have any preference for how you want the initial membership functions 
to be parameterized, you can let anfis do this for you. 


Automatic FIS Structure Generation with ANFIS 

To initialize your FIS using anfis : 

1 Choose Grid partition, the default partitioning method. (The two partition 
methods, grid partitioning and subtractive clustering, aredescribed later in 
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"Fuzzy C-Means Clustering" on page 2-120, and in "Subtractive Clustering" 
on page 2-123. 

2 Click on the Generate FIS button. This brings up a menu from which you 
can choose the number of membership functions, MFs, and the type of input 
and output membership functions. Notice there are only two choices for the 
output membershipfunction: constant and linear. This limitation of output 
membership function choices is because anfis only operates on Sugeno-type 
systems. 

3 Fill in the entries as we've done below, and click on OK. 



You can also implement this FIS generation from the command line using the 
command genfisl (for grid partitioning) or genfis2 (for subtractive 
clustering). A command line language example illustrating the use of genfisl 
and anfis is provided later. 

Specifying Your Own Membership Functions for ANFIS 

Although we don’t expect you to do this for this example, you can choose your 
own preferred membership functions with specific parameters to be used by 
anfis as an initial FI S for training. 
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To define your own FIS structure and parameters: 

1 Open the Edit membership functions menu item from the View menu. 

2 Add your desired membership functions (thecustom membership option will 
be disabled for anf is). The output membership functions must either be all 
constant or all linear. For carrying out this and the following step, see 'The 
FIS Editor" on page 2-49 and 'The Membership Function Editor" on page 
2-52. 

3 Select the Edit rules menu item in the View menu. Use the Rule Editor to 
generate the rules (see'The Rule Editor" on page 2-56). 

4 Select the Edit FIS properties menu item from the View menu. Name your 
FIS, and save it to either the workspace or the disk. 

5 Use the View menu to return totheANFIS Editor GUI totrain theFIS. 

To load an existing FI S for AN FI S initialization, in the Generate FIS portion 
of the GUI , click on Load from worksp. or Load from disk. You will load your 
FIS from the disk if you have saved an FIS previously that you would liketo 
use. Otherwise you will be loading your FIS from the workspace. Either of 
these radio buttons toggle the Generate FIS button to Load.... Load your FIS 
by clicking on this button. 
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Viewing Your FIS Structure 

After you generate the FI S, you can view the model structure by clicking on the 
Structure button i n the middle of the right-hand side of the GU I . A new GU I 
appears, as follows: 



Return to other open GUIs 
using the Window menu. 


Color coding of branches 
characterizes the rules. 


Node labels - for 
example, leftmost 
node is the 
input node^ 


Node representing a normalization 
factor for the rules. 


The branches in this nodal graph are color coded to indicate whether or not 
and, not, or or, are used in the rules. Clicking on the nodes indicates 
information about the structure. 

You can view the membership functions or the rules by opening either the 
Membership Function Editor, or the Rule Editor from the View menu. 
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ANFIS Training 

Thetwoanfis parameter optimization method options availablefor FIS 
training are hybrid (the default, mixed least squares and backpropagation) 
and backpropa (backpropagation). The Error Tolerance is used to create a 
training stopping criterion, which is related tothe error size. The training will 
stop after thetraining data error remains within this tolerance. This is best left 
set to 0 if you don't know how your training error is going to behave. 

To start thetraining: 

• Leave the optimization method at hybrid. 

• Set the n umber of trai n i ng epochs to 40, under the E pochs I i sti ng on the G U I 
(the default value is 3). 

• Select Train Now. 

The following should appear on your screen: 



*** Training error 
... Checking error 


Notice how the checking error decreases up to a certain point in thetraining 
and then it increases. This increase represents the point of model overfitting, 
anf is chooses the model parameters associated with the minimum checking 
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error (just prior tothis jump point). This is an examplefor which the checking 
data option of anfis is useful. 

Testing Your Data Against the Trained FIS 

To test your FIS against the checking data, click on Checking data in the Test 
FIS portion of the GUI, and click on Test Now. Now when you test the checking 
data agai nst the F I S it looks pretty good: 



Note on loading more data with anfis: If you are ever loading data into 
anfis after clearing previously loaded data, you must make sure that the 
newly loaded data sets havethesame number of inputs as the previously 
loaded ones did. Otherwise you will have to start a new anfisedit session 
from the command I i ne. 
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Note on the Checking Data option and Clearing Data: If you don't want to 
use the checking data option of anfis, don't load any checking data before you 
train the FI S. If you decide to retrain your FIS with no checking data, you can 
unload the checking data in one of two ways. One method is to click on the 
Checking radio button in the Load data portion of the GUI and then click on 
Clear Data to unload the checking data. The other method you can use is to 
close the GUI and go to the command line and retype anf isedit. I n this case 
you will have to reload the training data. After clearing the data, you will need 
to regenerate your FIS. Once the FIS is generated you can use your first 
training experienced decide on the number of training epochs you want for 
the second round of training. 


ANFIS Editor GUI Example 2: Checking Data Doesn't 
Validate Model 

I n this example, we examine what happens when the training and checking 
data sets are sufficiently different. We see how the AN FIS Editor GUI can be 
used to learn something about data sets and how they differ. 

1 Clear both the training and checking data. 

2 You can press the Clear Plot button on the right, although you don’t have to. 

3 Load fuzex2trnData and f uzex2chkData (respectively, the training data 
and checking data) from the MATLAB workspace just as you did in the 
previous example. 
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You should get something that looks likethis: 


i ANFIS Editor: Untitled 


File Edit View 
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train data loaded 

Help j Close 



+++ Checking data 
ooo Training data 
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Train the FI S for this system exactly as you did in the previous example, except 
now choose 60 Epochs before training. You should get the following: 



N otice the checki ng error i s quite I arge. 1 1 appears that the mi ni mum checki ng 
error occurs within the first epoch. Recall that using the checking data option 
with anfis automatically sets the FIS parameters to be those associated with 
the minimum checking error. Clearly this set of membership functions would 
not be the best choice for modeling the training data. 

What's wrong here? This example illustrates the problem discussed earlier 
wherein the checking data set presented to anfis for training was sufficiently 
different from the training data set. As a result, the trained F I S did not capture 
the features of this data set very well. This illustrates the i mportance of 
knowing the features of your data set well enough when you select your 
training and checking data. When this is not the case, you can analyze the 
checking error plots to see whether or not the checking data performed 
sufficiently well with the trained model. I n this example, the checking error is 
sufficiently large to indicate that either moredata needs to be selected for 
training, or you may want to modify your membership function choices (both 
the number of membership functions and the type). Otherwise the system can 
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be retrai ned without the checki ng data, if you think the trai ni ng data captures 
sufficiently the features you are trying to represent. 

Tocompletethis example, let'stest thetrained FI S model against the checking 
data. To do so, click on Checking data in the Test FIS portion of the GUI , and 
click on Test Now. The foil owing plot in the GUI indicates that there is quite a 
discrepancy between the checking data output and the FIS output. 



anfis from the Command Line 

As you can see, generating an FI S using the AN FI S Editor GUI is quite simple. 
Flowever, as you saw in the last example, you need to be cautious about 
implementing the checking data validation feature of anfis. You must check 
that the checking data error does what is supposed to. Otherwise you need to 
retrain the FI S. 

I n this section wedescri be how to carry out the command I i ne features of anfis 
on a chaotic times-series prediction example. 
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Using anfis for Chaotic Time Series Prediction 

The demo mgtsdemo uses anfis topredictati me series that isgenerated by the 
following Mackey-Glass (MG) time-delay differential equation: 


x(t) = 


0.2 x(t-x) 

1 + X 10 (t-T) 


0.1x(t) 


This time series is chaotic, and so there is no clearly defined period. The series 
will not converge or diverge, and the trajectory is highly sensitive to initial 
conditions. This is a benchmark problem in the neural network and fuzzy 
modeling research communities. 

To obtain the time series value at integer points, we applied the fourth-order 
Runge-Kutta method to find the numerical solution totheaboveMG equation; 
the result was saved in the file mgdata.dat. Here we assume x(0) =1.2, x =17, 
and x(t) =0for t <0. To plot the MG time series, type 

load mgdata.dat 

t = mgdata(:, 1); x = mgdata(:, 2); plot(t, x); 


Mac key -Glass chaotic time series 



I n time-series prediction we want to use known values of the time series up to 
the point in time, say, t, to predict the value at some point in the future, say, 
t+P. The standard method for this type of prediction is to create a mapping 
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from D sample data points, sampled every A units in time, (x(t-(D-l)A),...,x(t-A), 
x(t)), toa predicted future valuex(t-fP). Foil owing the conventional settings for 
predicting the MG time series, we set D =4 and A =P =6. For each t, the input 
training data for anfis is a four dimensional vector of the following form: 

w(t) = [x(t-18) x(t-12) x(t-6) x(t)] 

The output training data corresponds to the trajectory prediction: 

s(t) = x(t+6) 

For each t, ranging in values from 118 to 1117, the training input/output data 
will be a structure whose first component is the four-dimensional input w, and 
whose second component is the output s. There will be 1000 input/output data 
values. We use the first 500 data values for the anfis training (these become 
the trai ni ng data set), whi le the others are used as checki ng data for val idati ng 
the identified fuzzy model. This results in two 500-point data structures: 
trnData and chkData. 

FI ere is the code that generates this data: 

for t=1 18:1117, 

Data (t - 117, :) = [x(t-18) x(t-12) x(t-6) x(t) x(t+6)]; 
end 

trnData=Data(1 :500, :); 
chkData=Data(501 :end, :); 

To start the training, we need an FI S structure that specifies the structure and 
initial parameters of the FI S for learning. This is the task of genf isl : 

fismat = genf isl (trnData) ; 

Since we did not specify numbers and types of membership functions used in 
the FIS, default values are assumed. These defaults provide two generalized 
bell membership functions on each of the four inputs, eight altogether. The 
generated FIS structure contains 16 fuzzy rules with 104 parameters. I n order 
to achieve good generalization capability, it is important to have the number of 
training data points be several times larger than the number parameters being 
estimated. I n this case, the ratio between data and parameters is about five 
(500/104). 
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The function genfisl generates initial membership functions that are equally 
spaced and cover the whole input space. You can plot the input membership 
functions using the following commands. 

subplot(2,2, 1 ) 

plotmf (f ismat , 'input 1 , 1) 

subplot(2,2,2) 

plotmf (f ismat , 'input 1 , 2) 

subplot(2,2,3) 

plotmf (f ismat , 'input', 3) 

subplot(2,2,4) 

plotmf (f ismat , 'input', 4) 

These initial membership functions are shown below. 






To start the training, type 

[fismatl ,error1 ,ss,fismat2,error2] = ... 

anfis(trnData,fismat, [ ] , [ ] ,chkData) ; 

This takes about four minutes on a Sun SPARCstation 2 for 10 epochs of 
training. Because the checking data option of anfis was invoked, the final FIS 
you choose would ordinarily be the one associated with the minimum checking 
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error. This is stored in fismat2. The following code will plot these new 
membership functions: 

subplot(2,2, 1 ) 

plotmf (fismat2, 'input 1 , 1) 

subplot(2,2,2) 

plotmf (fismat2, 'input 1 , 2) 

subplot(2,2,3) 

plotmf (fismat2, 'input', 3) 

subplot(2,2,4) 

plotmf (fismat2, 'input', 4) 

Here is the result:. 






T o plot the error signals type 
plot( [errorl ; error2]); 

Here errorl and error2 arethe root mean squared error for thetraining and 
checking data, respectively. 

I n addition to these error plots, you may want to plot the FI S output versus the 
training or checking data. To compare the original MG time series and the 
fuzzy prediction side by side, try 

anfis_output = evalfis( [trnData; chkData] , fismat2); 
index = 125:1124; 
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subplot(21 1 ) , plot(t(index) , [x(index) anf is_output] ) ; 
subplot (21 2) , plot(t(index) , x(index) - anf is_output ) ; 


MG Time Serise and ANFIS Prediction 



Note that the difference between the original MG time series and the anf is 
estimated values is very small. This is why you can only see one curve in the 
first plot. The prediction error is shown in the second plot with a much finer 
scale. Notethat we have only trained for lOepochs. Better performance is 
expected if we apply more extensive training. 

More on anfis and the ANFIS Editor GUI 

The command anfis takes at least two and at most six input arguments. The 
general format is 

[fismatl ,trnError,ss,fismat2,chkError] = ... 

anf is (trnData,f ismat , trnOpt , dispOpt ,chkDat a, method) ; 

wheretrnOpt (training options), dispOpt (display options), chkData (checking 
data), and method (training method), areoptional. All of theoutput arguments 
are also optional . I n this section we discuss the arguments and range 
components of the command line function anfis, as well as the analogous 
functionality of the AN FIS Editor GUI. 

When theANFIS Editor GUI is invoked using anfisedit, only the training 
data set must exist prior to implementing anfis .in addition, the step-size will 
be fixed when the adaptive neuro-fuzzy system is trained using thisGUI tool. 
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Training Data 

The training data, trnData, is a required argument to antis, as well as to the 
AN FI S Editor GUI . Each row of trnData is a desired input/output pair of the 
target system to be modeled. Each row starts with an input vector and is 
followed by an output value. Therefore, the number of rows of trnData is equal 
tothe number of training data pairs, and, si nee there is only one output, the 
number of columns of trnData is equal tothe number of inputs plus one. 

Input FIS Structure 

The input FIS structure, f ismat, can be obtained either from any of the fuzzy 
editors: the FI S Editor, the Membership F unction Editor, and the Rule Editor 
from the ANFIS Editor GUI , (which allows an FIS structure to be loaded from 
the disk or the workspace), or from the command line function, genf isl (for 
which you only need to give numbers and types of membership functions). The 
FIS structurecontains both the model structure, (which specifies such items as 
the number of rules in the FI S, the number of membership functions for each 
input, etc.), and the parameters, (which specify the shapes of membership 
functions). There are two mdthods that anfis learning employs for updating 
membership function parameters: backpropagation for all parameters (a 
steepest descent method), and a hybrid method consisting of backpropagation 
for the parameters associated with the input membership functions, and least 
squares estimation for the parameters associated with the output membership 
functions. As a result, the training error decreases, at least locally, throughout 
the learning process. Therefore, the more the initial membership functions 
resemble the optimal ones, the easier it will be for the model parameter 
training to converge. FI uman expertise about the target system to be modeled 
may aid in setting up these initial membership function parameters in the FI S 
structure. 

Note that genf isl produces an FIS structure based on a fixed number of 
membership functions. This invokes the so-called curse of dimensionality, and 
causes an explosion of the number of rules when the number of inputs is 
moderately large, that is, more than four or five. The Fuzzy Logic Tool box offers 
a method that will provide for somedimension reduction in the fuzzy inference 
system: you can generate an FI S structure using the clustering algorithm 
discussed in "Subtractive Clustering" on page 2-123. FromtheANFIS Editor 
GUI , this algorithm is selected with a radio button beforetheFI S is generated. 
This subtractive clustering method partitions the data into groups called 
clusters, and generates an FIS with the minimum number rules required to 
distinguish the fuzzy qualities associated with each of the clusters. 
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Training Options 

The AN FIS Editor GUI tool allows you to choose your desired error tolerance 
and number of training epochs. 

Training option trnOpt for the command line anfis is a vector that specifies 
the stopping criteria and the step-size adaptation strategy: 

• trnOpt (1 ) : number of training epochs, default =10. 

• trnOpt ( 2 ) : error tolerance, default =0. 

• trnOpt (3): initial step-size, default=0.01. 

• trnOpt ( 4 ): step-size decrease rate, default =0.9. 

• trnOpt (5): step-size increase rate, default =1.1. 

If any element of trnOpt is an NaN or missing, then the default value is taken. 
The training process stops if the designated epoch number is reached or the 
error goal is achieved, whichever comes first. 

Usually we want the step-size profile to be a curve that increases initially, 
reaches some maximum, and then decreases for the remainder of the training. 
This ideal step-size profile is achieved by adjusting the initial step-size and the 
increase and decrease rates (trnOpt (3) - trnOpt (5)). The default values are 
set up to cover a wide range of learning tasks. For any specific application, you 
may want to modify these step-size options in order to optimize the training. 
However, as we mentioned previously, there are no user-specified step-size 
options for trai ni ng the adaptive neuro fuzzy i nference system generated usi ng 
theANFIS Editor GUI. 

Display Options 

Display options only apply to the command line function, anfis. 

Forthe command lineanf is, thedisplay options argument, dispOpt, isavector 
of either ones or zeros that specifies what information to display, (print in the 
MATLAB command line window), before, during, and after the training 
process. One is used to denote print this option, whereas zero denotes don't 
print this option. 

• disp0pt(l ): display ANFIS information, default =1. 

• dispOpt ( 2 ): display error (each epoch), default =1. 

• dispOpt (3) : display step-size (each epoch), default =1. 

• dispOpt ( 4 ): display final results, default =1. 
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The default mode is displays all available information. If any element of 
dispOpt is NaN or missing, the default value will betaken. 

Method 

Both the AN FI S Editor GUI and the command line anf is apply either a 
backpropagation form of the steepest descent method for membership function 
parameter estimation, or a combination of backpropagation and the 
least-squares method to estimate membership function parameters. The 
choices for this argument are hybrid or backpropagation. These method 
choices are designated in the command line function, anf is, by 1 and 0, 
respectively. 

Output FIS Structure for Training Data 

fismatl isthe output FIS structure corresponding to a minimal training error. 
This is the FI S structure that you will use to represent the fuzzy system when 
there is no checking data used for model crossvalidation. This data also 
represents the FI S structure that issaved bytheANFIS Editor GUI when the 
checking data option is not used. 

When the checking data option is used, the output saved isthat associated with 
the minimum checking error. 

Training Error 

The training error is the difference between the training data output value, 
and the output of the fuzzy inference system corresponding to the same 
training data input value, (the one associated with that training data output 
value). Thetraining error trnError records the root mean squared error 
(RMSE) of thetraining data set at each epoch, fismatl is the snapshot of the 
FIS structure when thetraining error measure is at its minimum. The ANF I S 
Editor GUI will plot thetraining error vs. epochs curve as the system is 
trained. 

Step-size 

You cannot control the step-size opt ions with theAN FIS Editor GUI. Using the 
command line anf is, the step-size array ss records the step-size during the 
training. Plotting ss gives the step-size profile, which serves as a reference for 
adjusting the initial step-size and the corresponding decrease and increase 


2-117 



rates. The step-size (ss) for the command line function antis is updated 
according to the following guidelines: 

• Ifthe error undergoes four consecutive reductions, increase the step-size by 
multiplying it by a constant (ssinc) greater than one. 

• If the error undergoes two consecutive combi nations of one increase and one 
reduction, decrease the step-size by multiplying it by a constant (ssdec) less 
than one. 

The default value for the initial step-size is 0.01; the default values for ssinc 
and ssdec are 1.1 and 0.9, respectively. All the default values can be changed 
via the training option for the command line antis. 

Checking Data 

The checking data, chkData, is used for testing the generalization capability of 
the fuzzy inference system at each epoch. The checking data has the same 
format as that of the training data, and its elements are generally distinct from 
those of the training data. 

The checking data is important for learning tasks for which the input number 
is large, and/or the data itself is noisy. I n general we want a fuzzy inference 
system to track a given input/output data set well. Si nee the model structure 
used for antis is fixed, there is a tendency for the model to overfit the data on 
which is it trained, especially for a large number of training epochs. If 
overfitting does occur, we cannot expect the fuzzy inference system to respond 
well toother independent data sets, especially if they are corrupted by noise. A 
validation or checking data set can be useful for these situations. This data set 
is used to crossvalidate the fuzzy inference model. This crossvalidation is 
accomplished by applying the checking data to the model, and seeing how well 
the model responds to this data. 

When the checking data option is used with antis, either via the command 
line, or using the AN FI S Editor GUI , the checking data is applied to the model 
at each training epoch. When the command line antis is invoked, the model 
parameters that correspond to the minimum checking error are returned via 
the output argument fismat2. The FIS membership function parameters 
computed using the AN FI S Editor GUI when both training and checking data 
are loaded are associ ated with the trai ni ng epoch that has a mi ni mum check i ng 
error. 
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The use of the minimum checking data error epoch to set the membership 
function parameters assumes 

• The checking data is similar enough to the training data that the checking 
data error wi 1 1 decrease as the trai ni ng begi ns 

• The checking data increases at some point in thetraining, after which data 
overfitting has occurred. 

As discussed in "ANFIS Editor GUI Example2: Checking Data Doesn’t 
Validate Model" on page 2-106, depending on the behavior of the checking data 
error, the resulting FIS may or may not be the one you should be using. 

Output FIS Structure for Checking Data 

The output of the command lineanfis, fismat2, is the output FIS structure 
with the minimum checking error. This is the FIS structure that should be 
used for further calculation if checking data is used for cross validation. 

Checking Error 

The checking error is the difference between the checking data output value, 
and the output of the fuzzy inference system corresponding to the same 
checking data input value, (the one associated with that checking data output 
value). The checking error chkError records the RMSE for the checking data 
at each epoch, f ismat2 is the snapshot of the FI S structure when the checking 
error is at its minimum. TheANFIS Editor GUI will plot the checking error vs. 
epochs curve as the system is trained. 
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2 


Tutorial 


Fuzzy Clustering 

Clustering of numerical data forms the basis of many classification and system 
modeling algorithms. The purpose of clustering is to identify natural groupings 
of data from a large data set to produce a concise representation of a system's 
behavior. The F uzzy LogicToolbox is equipped with some tools that allow you 
to find clusters in input-output training data. You can use the cl uster 
information to generate a Sugeno-type fuzzy inference system that best models 
the data behavior using a minimum number of rules. The rules partition 
themselves according to the fuzzy qualities associated with each of the data 
clusters. This type of FIS generation can be accomplished automatically using 
the command line function, genfis2. 

Fuzzy C-Means Clustering 

F uzzy c-means { FCM) is a data clustering technique wherein each data point 
belongs to a cluster to some degree that is specified by a membership grade. 
This technique was originally introduced byj im Bezdek in 1981 [Bez81] as an 
improvement on earlier clustering methods. It provides a method of how to 
group data points that populate some multidimensional space into a specific 
number of different clusters? 

The Fuzzy Logic Tool box command line function fern starts with an initial 
guess for the cluster centers, which are intended to mark the mean location of 
each cluster. The initial guess for these cl uster centers is most likely incorrect. 
Additionally, f cm assigns every data point a membership grade for each 
cluster. By iteratively updating the cl uster centers and the membership grades 
for each data point, fem iteratively moves the cl uster centers to the "right" 
location within a data set. This iteration is based on minimizing an objective 
function that represents the distance from any given data point to a cluster 
center weighted by that data point's membership grade. 

fern is a command line function whose output is a list of cluster centers and 
several membership grades for each data point. You can use the information 
returned by fem to help you build a fuzzy inference system by creating 
membership functions to represent the fuzzy qualities of each cluster. 


2-120 



Fuzzy Clustering 


An Example: 2-D Clusters 

Let's use some quasi-random two-dimensional data to illustrate how FCM 
clustering works. Load a data set and take a look at it. 

load fcmdata.dat 

plot (f cmdata( : ,1 ) ,fcmdata( : ,2) , ’o’ ) 
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Now we invoke the command line function, fcm, and ask it to find two clusters 
in this data set 

[center, U,objFcn] = fcm(fcmdata,2) ; 

Iteration count = 1, obj . fen = 8.941176 
Iteration count = 2, obj. fen = 7.277177 

and soon until the objective function is no longer decreasing much at all. 

The variable center contains the coordinates of the two cluster centers, U 
contains the membership grades for each of the data points, and obj Fen 
contains a history of the objective function across the iterations. 

The fcm function is an iteration loop built on top of several other routines, 
namely initfem, which initializes the problem, distf cm, which is used for 
distance calculations, and stepfem, which steps through one iteration. 
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Plotting the objective function shows the progress of the clustering, 
plot (obj Fen) 


objective function values 



Finally here is a plot displaying the two separate cl usters as classified by the 
fem routine. The foil owing figure is generated using: 

load fcmdata.dat 

[center, U, obj_fcn] = fcm(fcmdata, 2); 
maxU = max(U) ; 

indexl = find(U(1, :) == maxU); 
index2 = find(U(2, :) == maxU); 

line(fcmdata(index1 , 1), f cmdata(index1 , 2), 1 linestyle 1 , . . . 

1 none 1 , 1 marker 1 , 'o',' color ' , ' g ' ) ; 

line(fcmdata(index2, 1 ) ,fcmdata(index2,2) , 1 linestyle ' , . . . 


' none ' , ' marker ' , 

1 x 1 , ' color ' , 

1 r 1 ) ; 



hold on 

plot(center(1 ,1 ) 

,center(1 ,2) , 

' ko ' , 

' markersize ' 

, 15, ' LineWidth ' ,2 

plot(center(2,1 ) 

,center(2,2) , 

' kx ' , 

' markersize ' 

,15, 'LineWidth' ,2 
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Cluster centers are indicated in the figure below by the large characters. 



Subtractive Clustering 

Suppose we don't have a clear idea how many clusters there should be for a 
given set of data. Subtractive clustering, [Chi 94], is a fast, one-pass algorithm 
for estimating the number of clusters and the cluster centers in a set of data. 
The cluster estimates obtained from the subclust function can be used to 
initialize iterative optimization-based clustering methods (fcm) and model 
identification methods (likeanf is). The subclust function finds the clusters by 
using the subtractive clustering method. 

The gent is2 function builds upon the subclust function to provide a fast, 
one-pass method to take input-output training data and generate a 
Sugeno-type fuzzy i nference system that models the data behavior. 

An Example: Suburban Commuting 

I n this example we apply the gent is2 function to model the relationship 
between the number of automobi I e tri ps generated from an area and the area's 
demographics. Demographic and trip data are from 100 traffic analysis zones 
in New Castle County, Delaware. Five demographic factors are considered: 
population, number of dwelling units, vehicle ownership, median household 
income, and total employment. Hencethe model has five input variables and 
one output variable. 
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Load the data by typi ng 
tripdata 

subplot(2, 1 , 1 ) , plot(datin) 
subplot(2, 1 ,2) , plot(datout) 



tripdata creates several variables in the workspace. Of the original 100 data 
points, we will use 75 data points as training data (datin and datout) and 25 
data points as checking data, (as well as for test data to validate the model). 
The checking data input/output pairs aredenoted by chkdatin and chkdatout. 
The genf is2 function generates a model from data using clustering, and 
requires you to specify a cluster radius. The cluster radius indicates the range 
of influence of a cluster when you consider the data space as a unit hypercube. 
Specifying a small cluster radius will usually yield many small clusters in the 
data, (resulting in many rules). Specifying a large cluster radius will usually 
yield a few large clusters in the data, (resulting in fewer rules). The cluster 
radius is specified as the third argument of genf is2. Here we cal I the genf is2 
function using a cluster radius of 0.5. 

fismat=genfis2( datin, datout ,0.5) ; 

genf is2 is a fast, one-pass method that does not perform any iterative 
optimization. An FI S structure is returned; the model type for the FI S 
structure is a first order Sugeno model with three rules. We can useevalf is to 
verify the model. 

fuzout=evalfis(datin,fismat) ; 

trnRMSE=norm(f uzout-datout) /sqrt (length (fuzout) ) 



Fuzzy Clustering 


trnRMSE = 

0.5276 

The variable trnRMSE istheroot mean square error of the system generated by 
the training data. To validate the general izabi I ity of the model, we apply test 
data to the FIS. For this example, we use the checking data for both checking 
and testi ng the F I S parameters. 

chkf uzout=evalf is (chkdatin ,f ismat ) ; 

chkRMSE=norm(chkf uzout-chkdatout ) /sqrt (length (chkf uzout) ) 
chkRMSE = 

0.6170 

Not surprisingly, the model doesn't do quite as good a job on the testing data. 
A plot of the testing data reveals the difference. 

plot (chkdatout) 
hold on 

plot(chkfuzout, 'o' ) 
hold off 



At this point, we can usethe optimization capability of anf is to improve the 
model. First, we will try using a relatively short anfis training (50 epochs) 
without implementing the checking data option, but test the resulting FIS 
model against the test data. The command line version of this is as follows: 

fismat2=anfis( [datin datout] ,fismat, [50 0 0.1]); 

After the training is done, we type: 
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fuzout2=evalfis(datin,fismat2) ; 

trnRMSE2=norm(f uzout2-datout) /sqrt(length(f uzout2) ) 
trnRMSE2 = 

0.3407 

chkf uzout2=evalf is (chkdatin , f ismat2) ; 

chkRMSE2=norm(chkf uzout2-chkdatout) /sqrt(length(chkfuzout2) ) 
chkRMSE2 = 

0.5827 

The model has improved a lot with respect to the training data, but only a little 
with respect to the checking data. Here is a plot of the improved testing data. 

plot (chkdatout) 
hold on 

plot(chkfuzout2, 'o' ) 
hold off 



Here we see that genfis2 can be used as a stand-alone, fast method for 
generating a fuzzy model from data, or as a pre-processor to anf is for 
determining the initial rules. An important advantage of using a clustering 
method to find rules is that the resultant rules are more tailored to the input 
data than they are in an FI S generated without clustering. This reduces the 
problem of combinatorial explosion of rules when the input data has a high 
dimension (the dreaded curse of dimensionality). 

Overfitting 

Now let’s consider what happens if we carry out a longer (200 epoch) training 
of this system using anf is, including its checking data option. 
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[fismat3,trnErr,stepSize,fismat4,chkErr]= . . . 

anfis( [datin datout] ,fismat2, [200 0 


0 . !],[], 


[chkdatin chkdatout]); 


The long I ist of output arguments returns a history of the step-sizes, the RM SE 
versus the training data, andtheRMSE versus the checking data associated 
with each training epoch. 

ANFIS training completed at epoch 200. 

Minimal training RMSE = 0.326566 
Minimal checking RMSE = 0.582545 

This looks good. The error with the training data isthe lowest we've seen, and 
the error with the checking data is also lower than before, though not by much. 
This suggests that maybe we had gotten about as close as possible with this 
system al ready. M aybe we have even gone so far as to overfit the system to the 
training data. Overfitting occurs when we fit the fuzzy system to the training 
data so wel I that it no I onger does a very good job of fitti ng the checki ng data. 
The result is a loss of generality. A look at the error history against both the 
training data and the checking data reveals much. 




— Checking Error 

- - Training Error 


20 40 60 80 100 120 140 160 180 200 


Epochs 


Here we can see that the training error settles at about the 50 th epoch point. I n 
fact, the smallest value of the checking data error occurs at epoch 52, after 
which it increases slightly, even as anfis continues to minimize the error 
against the training data all the way to epoch 200. Depending on the specified 
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error tolerance, this plot al so i ndi cates the model 's abi I ity to general ize the test 
data. 

A Clustering GUI Tool 

There is also the Clustering GUI , which implements fcm and subclust, along 
with all of their options. Its use is fairly self-evident. 

The clustering GUI looks likethis, and is invoked using the command line 
function, f indcluster. 


File Edit Window Help 




Load a data set (*.dat) 
into your directory. 


Choose fcm 
or subtractive 
clustering method. 


Options change 
with method 


Start clustering 
the data. 


Save the value of 
the cluster center. 


You can invokefindcluster with a data set directly, in order toopen theGUI 
with a data set. The data set must have the extension .dat. For example, to 
load the data set, clusterdemo.dat, type f indcluster ( 1 clusterdemo.dat 1 ) . 

You use the pull-down tab under Method to change between fcm (fuzzy 
c-means) and subtractiv (subtractive clustering). More information on the 
options can be found in the entries for fcm on page 3-22, and subclust on page 
3-72, respectively. 
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TheClustering GUI works on multidimensional datasets, but only displays 
two of those dimensions. Use the pull-down tabs under X-axis and Y-axisto 
select which data dimension you want to view. 
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Stand-Alone C-Code Fuzzy Inference Engine 

I n the fuzzy /fuzzy directory of the tool box, you can find two C files, f ismain .c 
and f is . c, which are provided as the source codes for a stand-alone fuzzy 
inference engine. The stand-alone C-code fuzzy i nference engi ne can read an 
FI S file and an input data file to perform fuzzy inference directly, or it can be 
embedded in other external applications. 

To compile the stand-alone fuzzy inference engi neon a UNIX system, type 
% cc -0 -o f ismain f ismain. c -lm 

(Note that % is only symbolic of a UNIX prompt, and that you do not have to 
type f is. c explicitly, since it is included in f ismain. c.) Upon successful 
compilation, type the executable command to see how it works: 

% fismain 

This prompts the following message: 

% Usage: fismain data_file fis_file 

This means that fismain needs two files to do its work: a data file containing 
rows of input vectors, and an FI S file that specifies the fuzzy i nference system 
under consideration. 

For example, consider an FI S structure file named, mam2l . f is. We can prepare 
the input data file using MATLAB: 

[x, y] = meshgrid(-5:5, -5:5); 
input_data = [ x ( : ) y ( : ) ] ; 
save fis_in input_data -ascii 

This saves all the input data as a 121-by-2 matrix in the ASCI I file f is_in, 
where each row of the matrix represents an input vector. 

Now we can call the stand-alone code: 

% fismain fis_in mam21.fis 

This generates 121 outputs on your screen. You can direct the outputs to 
another file: 

% fismain fis in mam21.fis > fis out 
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Now the file f is_out contains a 121-by-l matrix. I n general, each row of the 
output matrix represents an output vector. The syntax of f ismain is similar to 
its M EX-file counterpart evalf is . m, except that all matrices are replaced with 
files. 

T o compare the results from the M ATL AB M EX-fi le and the stand-alone 
executable, type the following within MATLAB: 

fismat = readf is ( ' mam21 1 ) ; 

matlab_out = evalfis(input_data, fismat); 

load fis_out 

max(max(matlab_out - fis_out)) 
ans = 

4.9583e-13 

This ti ny difference comes from the I i mited length pri ntout i n the fi le f is_out . 
There are several things you should know about this stand-alone executable: 

• 1 1 is compati ble with both ANSI and K & R standards for C code, as long as 

STDC is defined in ANSI compilers. 

• Customized functions are not allowed in the stand-alone executable, so you 
are limited to the 11 membership functions that come with the toolbox, as 
well as other factory settings for AND, OR, IMP, and AGG functions. 

• f ismain. c contains only the main ( ) function and it is heavily documented 
for easy adaptation to other applications. 

• To add a new membership function or new reasoning mechanism into the 
stand-alone code, you need tochangethefilefis.c, which contains all the 
necessary functions to perform the fuzzy inference process. 

• For the Macintosh, the compiled command f ismain tries to find f ismain. in 
and f ismain .f is as input data and FIS description files, respectively. The 
output isstored in f ismain . out. These filenames aredefined within 
Macintosh-specific #define symbols in fismain.c and can be changed if 
necessary. 
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This section is designed to briefly explain some of the specialized terms that are 
derived from fuzzy logic. 

aggregation - the combination of the consequents of each rule in a Mamdani 
fuzzy i nference system in preparation for defuzzification. 

Adaptive Neuro-Fuzzy Inference System (AN FI S) - a technique for 
automatically tuning Sugeno-type inference systems based on training data. 

antecedent - the initial (or "if") part of a fuzzy rule. 

consequent - the final (or "then") part of a fuzzy rule. 

defuzzification - the process of transforming a fuzzy output of a fuzzy 
i nference system into a crisp output. 

degree of membership - the output of a membership function, this value is 
always limited to between 0 and 1. Also known as a membership value or 
membership grade. 

degree of fulfillment - see firing strength. 

firing strength - the degree to which the antecedent part of a fuzzy rule is 
satisfied. The firing strength may be the result of an AND or an OR operation, 
and it shapes the output function for the rule. Also known as degree of 
fulfillment. 

fuzzification - the process of generating membership values for a fuzzy 
variable using membership functions. 

fuzzy c-means clustering - a data clustering technique wherein each data 
point belongs to a cluster to a degree specified by a membership grade. 

fuzzy inference system (FIS) - the overal I name for a system that uses fuzzy 
reasoning to map an input space to an output space. 

fuzzy operators - AN D, OR, and N OT operators. These are also known as 
logical connectives. 

fuzzy set - a set which can contain elements with only a partial degree of 
membership. 

fuzzy singleton - a fuzzy set with a membership function that is unity at a one 
particular point and zero everywhere else. 



Glossary 


implication - the process of shaping the fuzzy set in the consequent based on 
the results of the antecedent in a Mamdani-type FI S. 

Mamdani-type inference - a type of fuzzy inference in which the fuzzy sets 
from the consequent of each rule are combined through the aggregation 
operator and the resulting fuzzy set is defuzzified to yield the output of the 
system. 

membership function (MF) - a function that specifies the degree to which a 
given input belongs to a set or is related to a concept. 

singleton output function - an output function that is given by a spike at a 
single number rather than a continuous curve. I n the Fuzzy Logic Tool box it is 
only supported as part of a zero-order Sugeno model. 

subtractive clustering - a technique for automatically generating fuzzy 
inference systems by detecting clusters in input-output training data. 

Sugeno-type inference - a type of fuzzy inference in which the consequent of 
each rule is a linear combination of the inputs. The output is a weighted linear 
combination of the consequents. 

T-conorm - (also known as S-norm) a two-input function that describes a 
superset of fuzzy union (OR) operators, including maximum, algebraic sum, 
and any of several parameterized T-conorms. 

T-norm - a two-input function that describes a superset of fuzzy intersection 
(AND) operators, including minimum, algebraic product, and any of several 
parameterized T-nor ms. 
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This section of the chapter contains brief descriptions of all the functions in the 
Fuzzy Logic Tool box. The foil owing tables contain the functions listed by topic. 

GUI Tools 


Function 

Purpose 

anf isedit 

ANFIS Editor GUI. 

fuzzy 

Basic FIS Editor. 

mf edit 

Membership Function Editor. 

ruleedit 

Rule Editor and parser. 

ruleview 

Rule Viewer and fuzzy inference diagram. 

surfview 

Output Surface Viewer. 


Membership Functions 


Function 

Purpose 

dsigmf 

Difference of two sigmoid membershi p functions. 

gauss2mf 

Two-sided Gaussian curve membership function. 

gaussmf 

Gaussian curve membership function. 

gbellmf 

Generalized bell curve membership function. 

pimf 

Pi-shaped curve membership function. 
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Function 

Purpose 

psigmf 

Product of two sigmoidal membership functions. 

smf 

S-shaped curve membership function. 

sigmf 

Sigmoid curve membership function. 

trapmf 

Trapezoidal membership function. 

trimf 

Triangular membership function. 

zmf 

Z-shaped curve membership function. 


FIS Data Structure Management 


Function 

Purpose 

addmf 

Add membership function to FI S. 

addrule 

Add ruletoFIS. 

addvar 

Add variable to FI S. 

defuzz 

Defuzzify membership function. 

evalf is 

Perform fuzzy inference calculation. 

evalmf 

Generic membership function evaluation. 

gensurf 

Generate FIS output surface. 

getf is 

Get fuzzy system properties. 

mf 2mf 

T ranslate parameters between functions. 

newt is 

Create new FIS. 

parsrule 

Parse fuzzy rules. 

plotf is 

Display FIS input-output structure. 

plotmf 

Plot all of the membership functions associated 
with a given variable. 
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Function 

Purpose 

readf is 

Load FIS from disk. 

rmmf 

Remove membership function from FIS. 

rmvar 

Remove variable from FIS. 

sett is 

Set fuzzy system properties. 

showf is 

Display annotated FIS. 

showrule 

Display FIS rules. 

writef is 

Save FIS to disk. 


Advanced Techniques 


Function 

Purpose 

anf is 

Training routine for a Sugeno-type FI S (MEX only). 

f cm 

Find clusters with FCM clustering. 

genf isl 

Generate FI S matrix using grid method. 

genf is2 

GenerateFIS matrix using subtractive clustering. 

subclust 

Find cluster centers with subtractive clustering. 


Simulink Blocks 


Function 

Purpose 

f uzblock 

Fuzzy logic controller blocks and demo blocks. 

sff is 

Fuzzy inference S-function. 
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Demos 


Function 

Purpose 

def uzzdm 

Defuzzification methods. 

f cmdemo 

FCM clustering demo (2-D). 

f uzdemos 

GUI for F uzzy Logic Tool box demos. 

gasdemo 

AN FI S demo for fuel efficiency using subclustering. 

juggler 

Ball-juggler with Rule Viewer. 

invkine 

1 nverse kinematics of a robot arm. 

irisf cm 

FCM clustering demo (4-D). 

noisedm 

Adaptive noise cancellation. 

slbb 

Ball and beam control (Simulink ). 

slcp 

Inverted pendulum control (Simulink ). 

sltank 

Water level control (Simulink). 

sltankrule 

Water level control with RuleViewer (Simulink). 

sltbu 

Truck backer-upper (Simulink only). 
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Purpose 

Synopsis 

Description 


Add a membership function to an FIS. 

a = addmf (a, ' varType ', varlndex ,' mf Name ',' mf Type mfParams) 

A membership function can only be added toa variable in an existing MATLAB 
workspace FI S. I ndices are assigned to membership functions in the order in 
which they are added, so the first membership function added toa variable will 
always be known as membership function number one for that variable. You 
cannot add a membership function to input variable number two of a system if 
only one input has been defined. 

The function requires six input arguments in this order: 

1 A MATLAB variable name of a FI S structure in the workspace 

2 A string representing the type of variable you want to add the membership 
function to ('input 1 or 'output ' ) 

3 The index of the variable you want to add the membership function to 

4 A string representing the name of the new membership function 

5 A stri ng representi ng the type of the new membershi p functi on 

6 The vector of parameters that specify the membership function 
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Example 


See Also 


a=newfis( 'tipper' ) ; 

a=addvar (a, ' input ' , ' service ' , [0 10] ) ; 

a=addmf (a, 1 input ' , 1 , ' poor 1 , ' gaussmf ',[1.5 0] ) ; 

a=addmf (a, 1 input ' , 1 , ' good ' , ' gaussmf ',[1.5 5] ) ; 

a=addmf (a, 1 input ' , 1 , ' excellent 1 , 1 gaussmf ',[1.5 10] ) ; 

plotmf (a, ' input ' , 1 ) 



addrule, addvar, plotmf, rmmf, rmvar 
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addrule 


Purpose 

Synopsis 

Description 


Example 


See Also 


Add a ruleto an FIS. 
a = addrule(a, ruleList) 

addrule has two arguments. The first argument istheMATLAB workspace 
variableFIS name. The second argument for addrule isa matrixof oneor more 
rows, each of which represents a given rule. Theformat that the rule list matrix 
must take is very specific. If there are m inputs to a system and?; outputs, there 
must be exactly m + ?? +2 columns to the rule list. 

The first ??? col umns refer to the inputs of the system. E ach column contains a 
number that refers tothe index of the membership function for that variable. 

The next ?? columns refer to the outputs of the system. Each column contains a 
number that refers tothe index of the membership function for that variable. 

The??? +?? +1 column contains the weight that is to be applied tothe rule. The 
weight must be a number between zero and one, and is generally left as one. 

The??? +?? +2 column contains a 1 if the fuzzy operator for the rule's antecedent 
is AND. It contains a 2 if the fuzzy operator is OR. 

ruleList=[ 

11111 

12211 ]; 

a = addrule(a, ruleList) ; 

If the above system a has two inputs and one output, the first rule can be 
interpreted as: "If input 1 is MF 1 and input 2 is MF 1, then output 1 is MF 1." 

addmf, addvar, rmmf, rmvar, parsrule, showrule 
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addvar 


Purpose 

Synopsis 

Description 


Example 


See Also 


Add a variable to an FIS. 

a = addvar (a, 1 varType varName 1 , varBounds) 

addvar has four arguments in this order: 

1 The name of a FI S structure in the MATLAB workspace 

2 A string representing the type of the variable you want to add ('input 1 or 
'output ' ) 

3 A string representing the name of the variable you want to add 

4 The vector describing the limiting range values for the variable you want to 
add 

I ndices are applied to variables in the order in which they are added, so the 
first input variable added to a system will always be known as input variable 
number one for that system. I nput and output variables are numbered 
independently. 

a=newfis( 'tipper' ) ; 

a=addvar(a, ' input ' , ' service ' , [0 10] ) ; 

getf is (a, ' input ' , 1 ) 

MATLAB replies 

Name = service 
NumMFs = 0 
MFLabels = 

Range = [0 10] 

addmf, addrule, rmmf, rmvar 
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anfis 


Purpose 

Synopsis 


Description 


Training routine for Sugeno-type FIS (ME X only). 

[fismat,error1 ,stepsize] = anfis (trnData) 

[fismat,error1 , stepsize] = anfis(trnData,fismat) 

[fismatl ,error1 , stepsize] = ... 

anfis (trnData , f ismat , trnOpt , dispOpt ) 

[fismatl ,error1 , stepsize, fismat2,error2] = ... 

anf is (trnData, trnOpt, dispOpt ,chkData) 

[fismatl ,error1 , stepsize, fismat2,error2] = ... 
anf is ( t rnData, t rnOpt, dispOpt ,chkData,optMethod) 

This is the major training routine for Sugeno-type fuzzy inference systems, 
anfis uses a hybrid learning algorithm to identify parameters of Sugeno-type 
fuzzy i nference systems. It applies a combination of the least -squares method 
and the backpropagation gradient descent method for training FIS 
membership function parameters to emulate a given training data set. anfis 
can also be invoked using an optional argument for model validation. The type 
of model validation that takes place with this option is a checking for model 
overfitting, and the argument is a data set called the checking data set. 

The arguments in the above description for anfis are as follows: 

• trnData: the name of a training data set. This is a matrix with all but the 
last column containing input data, whilethe last column contains a single 
vector of output data. 

• f ismat: the name of an FIS, (fuzzy inference system) used to provide anfis 
with an initial set of membership functions for training. Without this option, 
anfis will use genf isl to implement a default initial FI S for training. This 
default FIS will have two membership functions of the Gaussian type, when 
invoked with only one argument. If f ismat is provided as a single number (or 
a vector), it is taken as the number of membership functions (or the vector 
whose entries are the respective numbers of membership functions 
associated with each respective input when these numbers differ for each 
input). In this case, both arguments of anf is are passed to genf isl to 
generate a valid FI S structure before starting the training process. 


3-10 



anfis 


• trnOpt: vector of training options. When any training option is entered as 
NaN the default options will be in force. These options are as follows: 

trnOpt(i ) : training epoch number (default: 10) 
trnOpt (2) : training error goal (default: 0) 
trnOpt (3): initial step size (default: 0.01) 
trnOpt (4): step size decrease rate (default: 0.9) 
trnOpt (5) : step size increase rate (default: 1.1) 

• dispOpt: vector of display options that specify what message to display in the 
MATLAB command window during training. The default value for any 
display option is 1 , which means the corresponding information is displayed. 
A o means the corresponding information is not displayed on the screen. 
When any display option is entered as NaN, the default options will be in 
force. These opti ons are as fol lows: 

dispOpt (1 ) : AN FI S information, such as numbers of input and output 
membership functions, and soon (default: 1) 
dispOpt (2): error (default: 1) 

dispOpt (3) : step size at each parameter update (default: 1) 
dispOpt (4): final results (default: 1) 

• chkData: the name of an optional checking data set for overfitting model 
validation. This data set is a matrix in the same format as the training data 
set. 

• optMethod: optional optimization method used in membership function 
parameter training: either 1 for the hybrid method or o for the 
backpropagation method. The default method is the hybrid method, which is 
a combination of least squares estimation with backpropagation. The default 
method is invoked whenever the entry for this argument is anything but o. 

The training process stops whenever the designated epoch number is reached 
or the training error goal is achieved. 
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Note on antis arguments: When anf is is invoked with two or more 
arguments, any optional arguments will takeon their default values if they 
are entered asNaNsor empty matrices. Default values can be changed directly 
by modifying the file anf is .m. Either NaNs or empty matrices must be used as 
place-holders for variables if you don’t want to specify them, but do want to 
specify succeeding arguments, for example, when you implement the checking 
data option of anf is. 


The range variables in the above description for anf is are as follows: 

• fismatl is the FIS structure whose parameters are set according to a 
minimum training error criterion. 

• errorl or error2 is an array of root mean squared errors representing the 
training data error signal and the checking data error signal, respectively. 

• stepsize isan array of step sizes. Thestep size is decreased (by multiplying 
it with the component of the training option corresponding to the step size 
decrease rate) if the error measure undergoes two consecutive combi nations 
of an increase foil owed by a decrease. Thestep size is increased (by 
multiplying it with the increase rate) if the error measure undergoes four 
consecutive decreases. 

• f ismat2 is the F I S structure whose parameters are set according to a 
minimum checking error criterion. 



anfis 



Example x = ( 0 : 0 . 1 : 10 )'; 

y = sin(2*x) . /exp(x/5) ; 
trnData = [x y ] ; 
numMFs = 5; 
mfType = 1 gbellmf 1 ; 
epoch_n = 20; 

in_fismat = genfisl (trnData, numMFs, mfType) ; 
out_fismat = anfis(trnData,in_fismat,20) ; 
plot(x,y,x,evalfis(x,out_fismat) ) ; 
legend( 'Training Data 1 , 'ANFIS Output 1 ); 

See Also genfisl, anfis 

References J ang, J ,-S. R., "Fuzzy Modeling Using Generalized Neural Networks and 

Kalman Filter Algorithm,” Proc. oftheNinth National Conf. on Artificial 
I ntdligence(AAAI-91), pp. 762-767, J uly 1991. 

J ang, J ,-S. R., "ANFIS: Adaptive-Network-based Fuzzy I nference Systems," 
IEEE Transactions on Systems, Man, and Cybernetics, Vol. 23, No. 3, pp. 
665-685, May 1993. 
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anfisedit 


Purpose To open the AN FIS Editor GUI. 

Synopsis anfisedit ( 'a' ) 

anfisedit(a) 
anfisedit 

Description 



Using anfisedit, you bring up the AN FI S Editor GUI from which you can load 
a data set and train anf is. The AN FI S Editor GUI invoked using 
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anfisedit 


Menu Items 


See Also 


anf isedit ( 1 a 1 ), brings up the AN FI S Editor GUI from which you can 
implement anfis using a FIS structure stored as a fileon your disk called, 
a.fis. 

anf isedit (a) operates the same way for a FI S structure a, stored as a variable 
in the MATLAB workspace. 

Refer to "anfis and the AN FI S Editor GUI” on page 2-92 for more information 
about how to use anfisedit . 

On the AN FI S Editor GUI , there is a menu bar that allows you to open related 
GUI tools, open and save systems, and soon. The File menu is the same as the 
one found on the FI S Editor. Refer to fuzzy on page 3-29 for more information. 

• Use the following Edit menu item: 

Undo to undo the most recent change. 

• Use the foil owing View menu items: 

Edit FIS properties... to invoke the FI S Editor. 

Edit rules... to invokethe Rule Editor. 

Edit member shi p fu ncti ons... to i n voke t he M embersh i p F u net i on E d i tor . 
View rules... to invokethe Rule Viewer. 

View surface... to invokethe Surface Viewer. 

fuzzy, mfedit, ruleedit, ruleview, surfview 
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convertfis 


Purpose 

Synopsis 

Description 


Convert a Fuzzy Logic Tool box version 1.0 FI S matrix to a version 2.0 FI S 
structure. 

f is_new=convertf is (f is_old) 

convertfis takes a version 1.0 FIS matrix and converts it to a version 2.0 
structure. 
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defuzz 


Purpose 

Synopsis 

Description 


Examples 


Defuzzify membership function, 
out = defuzz(x,mf ,type) 

defuzz (x,mf ,type) returns a defuzzified value out, of a membership function 
mf positioned at associated variablevaluex, using one of several defuzzification 
strategies, according tothe argument, type. The variabletype can beoneofthe 
following. 

• centroid: centroid of area method 

• bisector: bisector of area method 

• mom: mean of maximum method 

• som: smallest of maximum method 

• lorn: largest of maximum method 

If type is not one of the above, it is assumed to be a user-defined function, x and 
mf are passed to this function to generate the defuzzified output. 

x = -10:0.1:10; 

mf = trapmf (x, [-10 -8 -4 7]); 

xx = def uzz (x, mf centroid ') ; 
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dsigmf 


Purpose 

Synopsis 

Description 


Example 


See Also 


Built-in membership function composed of the difference between two 
sigmoidal membership functions. 


y = dsigmf (x,[a1 cl a2 c2]) 


The sigmoidal membership function used here depends on the two parameters 
a and c and is given by 


f(x;a, c) 


1 + e 


a(x-c) 


The membership function dsigmf depends on four parameters, al, cl, a2, and 
c2, and is the difference between two of these sigmoidal functions: 

f 2 (x; al,C])- f 2 (x; a2, c 2 ) 

The parameters are listed in the order: [al c 2 a 2 c 2 ]. 


x=0:0. 1:10; 

y=dsigmf (x, [5257]); 
plot(x,y) 

xlabel( 'dsigmf , P=[5 2 5 7]') 



gaussmf , gauss2mf, gbellmf , evalmf, mf2mf , pimf , psigmf , sigmf , smf, trapmf , 
trimf, zmf 
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evalfis 


Purpose 

Synopsis 

Description 


Perform fuzzy inference calculations. 

output= evalfis(input,fismat) 
output= evalfis(input,fismat, numPts) 

[output, IRR, ORR, ARR] = evalfis ( input , fismat ) 

[output, IRR, ORR, ARR] = evalfis ( input , fismat , numPts) 

evalfis has the following arguments: 

• input: a number or a matrix specifying input values. If input is an M-by-N 
matrix, where N is number of input variables, then evalfis takes each row 
of input as an input vector and returns the M-by-L matrix to the variable, 
output, where each row is an output vector and L is the number of output 
variables. 

• fismat: an FIS structure to be evaluated. 

• numPts: an optional argument that represents the number of sample points 
on which to evaluate the membership functions over the input or output 
range. If this argument is not used, the default value of 101 point is used. 

The range labels for evalfis are as follows: 

• output: the output matrix of size M-by-L, where M represents the number of 
input values specified above, and L is the number of output variables for the 
FIS. 

The optional range variables for evalfis are only calculated when the input 
argument is a row vector, (only one set of inputs is applied). These optional 
range variables are: 

• IRR: the result of evaluating the input values through the membership 
functions. This is a matrix of size numRules-by-N , where numRules is the 
number of rules, and N is the number of input variables. 

• ORR: the result of evaluating the output values through the membership 
functions. This is a matrix of size numPts-by-numRu/es*L, where numRules 
is the number of rules, and L is the number of outputs. The first numRules 
columns of this matrix correspond to the first output, the next numRules 
columns of this matrix correspond to the second output, and so forth. 

• ARR:thenumPts-by-L matrixoftheaggregatevalues sampled at numPts along 
the output range for each output. 
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evalfis 


Example 


See Also 


When invoked with only one range variable, this function computes the output 
vector, output , of the fuzzy inference system specified by the structure, 
f ismat, for the input value specified by the number or matrix, input. 

fismat = readfis( 'tipper 1 ) ; 
out = evalfis([2 1; 4 9], fismat) 

This generates the response 

out = 

7.0169 

19.6810 

ruleview, gensurf 
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evalmf 


Purpose 

Synopsis 

Description 


Examples 


See Also 


Generic membership function evaluation, 
y = evalmf (x, mfParams, mfType) 

evalmf evaluates any membership function, where x is the variable range for 
the membership function evaluation, mfType is a membership function from 
the tool box, and mfParams are appropriate parameters for that function. 

If you want tocreateyour own custom membership function, evalmf will still 
work, because it evaluates any membership function whose name it doesn’t 
recognize. 


x=0:0. 1:10; 
mfparams = [2 4 6] ; 
mftype = 1 gbellmf 1 ; 
y=evalmf(x, mfparams, mftype) ; 
plot(x,y) 

xlabel( 'gbellmf , P=[2 46]') 



dsigmf , gaussmf , gauss2mf , gbellmf, mf2mf , pimf , psigmf , sigmf , smf, trapmf , 
trimf, zmf 
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fcm 


Purpose 

Synopsis 

Description 


F uzzy c-means cl usteri ng. 

[center, U,obj_fcn] = fcm(data, cluster_n) 

[center, U, obj_fcn] = fcm(data, cluster_n) applies the fuzzy c-means 
clustering method to a given data set. 

The input arguments of this function are: 

• data: data set to be clustered; each row is a sample data point 

• cluster_n: number of clusters (greater than one) 

The output arguments of this function are: 

• center: matrix of final cluster centers where each row provides the center 
coordinates 

• U: final fuzzy partition matrix (or membership function matrix) 

• obj_fcn: values of the objective function during iterations 

fcm(data,cluster_n, options) uses an additional argument variable, 
options, to control clustering parameters, introduce a stopping criteria, and/ 
or set the iteration information display: 

options (1 ) : exponent for the partition matrix u (default: 2.0) 
options(2): maximum number of iterations (default: 100) 
options(3): minimum amount of improvement (default: le-5) 
options (4): info display during iteration (default: 1) 

If any entry of options is NaN, the default value for that option is used instead. 
The clustering process stops when the maximum number of iterations is 
reached, or when the objective function improvement between two consecutive 
iterations is less than the minimum amount of improvement specified. 
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Example 


fcm 


data = rand(100, 2) ; 

[center, U,obj_fcn] = fern (data, 2); 
plot (data( : , 1 ) , data( : , 2) , 1 o 1 ) ; 
maxU = max (U) ; 

indexl = find ( U ( 1 , : ) == maxU); 
index2 = find(U(2, :) == maxU); 


line (data(index1 , 1 ) , 

' marker 1 , ' * ' , 
'color 1 , ' g ' ) ; 

data; 

[indexl , 

2), 

1 linestyle 1 , 

' none ' 

line (data(index2, 1 ) , 

' marker 1 , ' * 1 , 

1 color 1 , 1 r 1 ) ; 

data; 

[ index2 , 

2), 

1 linestyle 1 , 

' none ' 
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findcluster 


Purpose 

Synopsis 

Description 


I nteractive clustering GUI for fuzzy c-means and subclustering, 
findcluster 

f indcluster ( 1 file . dat 1 ) 

findcluster brings up a GUI to implement fuzzy c-means (fcm) and/or fuzzy 
subtractive clustering (subtractiv) using the pull-down tab under Method on 
the GUI . Data is entered using the Load Data button. The options for each of 
these methods are set to default values. These can be changed. A description of 
the options for fuzzy c-means is found in fcm on page 3-22. A description of the 
options for fuzzy subclustering is found in subclust on page 3-72. 

This tool works on multidimensional data sets, but only displays two of those 
dimensions. Use the pull-down tabs under X-axis and Y-axis to select which 
data dimension you want to view. For example, if you have data that is 
five-dimensional, this tool labels the data as data_l, data_2, data_3, data_4, 
data_5, in the order in which the data appears in the data set. Start will 
perform the clustering, and Save Center will save the cluster center. 

When operating on a data set, file.dat, findcluster (file.dat) loads the 
data set automatically, plotting up tothe first twodimensions of thedata only. 
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findcluster 


Example 
See Also 


You can still choose which two dimensions of the data you want to cluster after 
the GUI comes up. 



findcluster ( 'clusterdemo.dat 1 ) 
f cm, subclust 
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fuzarith 


Purpose 

Synopsis 

Description 

Remark 

Example 


To perform fuzzy arithmetic. 

C = fuzarith(X, A, B, operator) 

Using interval arithmetic, C = fuzarith (X, A, B, operator) returns a fuzzy set 
C as the result of applying the function represented by the string, operator, 
that performs a bi nary operation on the sampled convex fuzzy sets A and B . The 
elements of A and B are derived from convex functions of the sampled universe, 

x. 


• A, B, and x are vectors of the same dimension. 

• operator is one of the foil owing strings: 1 sum 1 , 1 sub 1 , 1 prod 1 , and 1 div 1 . 

• The returned fuzzy set C is a column vector with the same length as x. 

Fuzzy addition might generate the message: divideby zero, but this will not 
affect the correctness of this function. 

point_n = 101;% this determines MF's resolution 

min_x = -20; max_x = 20;% universe is [min_x, max_x] 
x = linspace (min_x , max_x, point_n)'; 

A = trapmf(x, [-10 -2 1 3]);% trapezoidal fuzzy set A 
B = gaussmf(x, [2 5]);% Gaussian fuzzy set B 
Cl = fuzarith(x, A, B, 'sum'); 
subplot (2, 1 , 1 ) ; 

plot(x, A, ' b - - ' , x, B, 'm:', x, Cl, ' c ' ) ; 
title( 'fuzzy addition A+B'); 

C2 = fuzarith(x, A, B, 'sub'); 
subplot(2,1 ,2) ; 

plot(x, A, ' b - - ' , x, B, 'm:', x, C2, ' c ' ) ; 
title( 'fuzzy subtraction A- B ' ) ; 

C3 = fuzarith(x, A, B, 'prod'); 
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fuzblock 


Purpose 

Synopsis 

Description 


See Also 


Simulink fuzzy logic controller block, 
fuzblock 

This command brings up a Simulink system that, in addition to some Si mu I ink 
demo blocks, contains two Simulink blocks you can use: 

• The Fuzzy Logic Controller 

• The Fuzzy Logic Controller With Rule Viewer, (see also ruleview on page 
3-61). This block forces the Rule Viewer to pop open during a Simulink 
simulation. 

The dialog box associated with either of these blocks is found by double-clicking 
on the Fuzzy Logic Controller block. This box contains the name of the FI S 
structure in the workspace that corresponds to the desired fuzzy system you 
want in your Simulink model. 

To open this dialog box for the Fuzzy Logic Controller With Rule Viewer block, 
you have to 

1 Double-click on this block, and a Simulink diagram with a Fuzzy Logic 
Controller block opens. 

2 Double-click on the second Fuzzy Logic Controller block that pops open. 

If the fuzzy i nference system has multiple inputs, theseinputs should be 
multiplexed together before feeding them intoeither the Fuzzy Logic 
Controller or the Fuzzy Logic Controller With RuleViewer block. Similarly, if 
the system has multiple outputs, these signals will be passed out of the block 
on one multiplexed line. 



sffis, ruleview 
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fuzdemos 


Purpose 

Synopsis 

Description 


List of all Fuzzy Logic Tool box demos, 
fuzdemos 

This function brings up a GUI that allows you to choose between any of several 
Fuzzy Logic Tool box demos listed under "Demos" on page 3-5. 
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Purpose 

Synopsis 


The Diagram 


fuzzy 


To invoke the basic FIS editor, 
fuzzy 

f uzzy (f ismat) 



This GUI tool allows you to edit the highest level features of the fuzzy inference 
system, such as the number of input and output variables, the defuzzification 
method used, and soon. Refer to'TheFIS Editor"on page 2-49 and ff., for more 
information about howto use the GUIs associated with fuzzy. 

TheFIS Editor is the high-level display for anyfuzzy logic inference system. It 
allows you to cal I the various other editors to operate on the FI S. This interface 
allows convenient access to all other editors with an emphasis on maximum 
flexibility for interaction with the fuzzy system. 

The diagram displayed at the top of the window shows the inputs, outputs, and 
a central fuzzy rule processor. Click on one of the variable boxes to make the 
selected box the current variable. You should seethe box highlighted in red. 
Double-click on one of the variables to bring up the Membership Function 
Editor. Double-click on the fuzzy rule processor to bring up the Rule Editor. If 
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fuzzy 


Menu Items 


a variable exists but is not mentioned in the rule base, it is connected tothe 
rule processor block with a dashed rather than a solid line. 

The FIS Editor displays a menu bar that allows you to open related GUI tools, 
open and save systems, and so on. 

• Under File select: 

New Mamdani FIS... to open a new Mamdani -style system with no 
variables and no rules called Untitled. 

New Sugeno FIS... to open a new Sugeno-style system with no variables 
and no rules called Untitled. 

Open from disk... to load a system from a specified . f is file on disk. 
Save to disk to save the current system to a . f is file on disk. 

Save to disk as... to save the current system to disk with the option to 
rename or relocate the file. 

Open from workspace... to load a system from a specified F I S structure 
variable in the workspace. 

Save to workspace... to save the system to the currently named F I S 
structure variable in the workspace. 

Save to workspace as... to save the system to a specified F I S structure 
variable in the workspace. 

Close window to close the GUI . 

• Under Edit select: 

Add input to add another input to the current system. 

Add output to add another output to the current system. 

Remove variable to delete a selected variable. 

Undo to undo the most recent change. 

• Under View select: 

Edit MFs... to invoke the Membership Function Editor. 

Edit rules... to invokethe Rule Editor. 

Edit anfis... toi nvoke the AN FI S Editor for single output Sugeno systems 
only. 

View rules... to invokethe Rule Viewer. 

View surface... to i nvoke the Surface Viewer. 
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Inference 
Method Pop-up 
Menus 


See Also 


fuzzy 


Five pop-up menus are provided to change the functionality of the five basic 

steps in the fuzzy implication process: 

• And method: Choose min, prod, or Custom, for a custom operation. 

• Or method: Choose max, probor (probabilistic or), or Custom, for a custom 
operation. 

• Implication method: Choose min, prod, or Custom, for a custom operation. 
This selection is not available for Sugeno-style fuzzy inference. 

• Aggregation method: Choose max, sum, probor, or Custom, for a custom 
operation. This selection is not available for Sugeno-style fuzzy inference. 

• Defuzzification method: For Mamdani -style inference, choose centroid, 
bisector, mom (middle of maximum), som (smallest of maximum), lorn 
(largest of maximum), or Custom, for a custom operation. For Sugeno-style 
inference, choose between wtaver (weighted average) or wtsum (weighted 
sum). 

mfedit, ruleedit, ruleview, surfview, anfisedit 
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gauss2mf 


Purpose 

Synopsis 

Description 


Gaussian combination membership function, 
y = gauss2mf (x, [ sig 1 cl sig2 c2]) 

The Gaussian function depends on two parameters sig and c as given by 

-(x-c) 2 

f 2 a 2 

f(x\a , c) = e 


The function gauss2mf is a combination of two of these. The first function, 
specified by sig I andc7, determines the shape of the leftmost curve. The second 
function specified by sig2 and c2 determines the shape of the right-most curve. 
Whenever cl <c2, thegauss2mf function reaches a maximum valueof 1. 
Otherwise, the maximum value is less than one. The parameters are listed in 
the order: 

[sigl, cl, sig2, c2\ 
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gauss2mf 


Examples 


X 


( 0 : 0 . 1 : 10 ) ' ; 

yl = gauss2mf(x, [2 4 1 8]) 
y2 = gauss2mf(x, [2517]) 
y3 = gauss2mf(x, [2616]) 
y4 = gauss2mf(x, [2715]) 
y5 = gauss2mf(x, [2814]) 
plot ( x , [yl y2 y3 y4 y5 ] ) ; 
set (get, 'name 1 , 'gauss2mf' 


1 numbertitle 1 , 1 off 1 ) ; 



See Also 


dsigmf, gauss2mf, gbellmf , evalmf, mf2mf, pimf , psigmf, sigmf , smf, trapmf, 
trimf, zmf 


3-33 



gaussmf 


Purpose 

Synopsis 

Description 


Example 


See Also 


Gaussian curve built-in membership function, 
y = gaussmf (x, [sig c] ) 

The symmetric Gaussian function depends on two parameters a and c as given 
by 

-( X-C ) 2 

,, , 2 o 2 

f(X;c>, C) = e 


The parameters for gaussmf represent the parameters o and c listed in order 
in the vector [sig c]. 

x=0:0. 1:10; 
y=gaussmf (x, [2 5] ) ; 
plot(x,y) 

xlabel( 'gaussmf , P=[2 5]') 



dsigmf , gaussmf, gbellmf, evalmf , mf2mf , pimf, psigmf , sigmf, smf, trapmf , 
trimf, zmf 
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gbellmf 


Purpose 

Synopsis 

Description 


Example 


See Also 


Generalized bell-shaped built-in membership function, 
y = gbellmf (x,params) 

The generalized bell function depends on three parameters a, b, and cas given 
by 


f(x;a, b, c) = Th 

1 + 

a 

where the parameter b is usually positive. The parameter c locates the center 
of the curve. Enter the parameter vector params, the second argument for 
gbellmf, as the vector whose entries are a, b, and c, respectively. 

x=0:0. 1:10; 
y=gbellmf (x, [2 4 6] ) ; 
plot (x,y) 

xlabel( 'gbellmf , P=[2 46]') 



dsigmf, gaussmf, gauss2mf , evalmf, mf2mf, pimf , psigmf, sigmf , smf, trapmf, 
trimf, zmf 
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genfisl 


Purpose 

Synopsis 

Description 


Generate an FIS structure from data without data clustering, 
fismat = genfisl (data) 

fismat = genfisl (data, numMFs , inmf type , outmftype) 

genfisl generates a Sugeno-type FI S structure used as initial conditions 
(initialization of the membership function parameters) for anfis training, 
genf isl(data, numMFs, inmf type, outmftype) generates a F I S structure from a 
training data set, data, using a grid partition on the data (no clustering). 

The arguments for genfisl are as follows: 

• data isthetrainingdata matrix, which must be entered with all but the last 
columns representing input data, and the last column representing the 
single output. 

• numMFs is a vector whose coordinates specify the number of membership 
functions associated with each input. If you want the same number of 
membership functions to be associated with each input, then it suffices to 
make numMFs a single number. 

• inmftype is a string array in which each row specifies the membership 
function type associated with each input. Again, this can be a 
one-dimensional singlestring if the type of membership functions associated 
with each input is the same. 

• outmftype is a stri ng that specifies the membershi p function type associated 
with the output. There can only be one output, since this is a Sugeno-type 
system. The output membership function type must be either linear or 
constant. 

The number of membership functions associated with the output is the same 
as the number of rules generated by genfisl . The default number of 
membership functions, numMFs, is 2; the default input or output membership 
function type is 'gbellmf '. These are used whenever genfisl is invoked 
without the last three arguments. 
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genfisl 


Example 


See Also 


data = [rand(10,1) 10*rand(10, 1 ) -5 rand (10,1)]; 
numMFs = [3 7] ; 

mfType = str2mat ( ' pimf 1 , 1 trimf ' ) ; 
fismat = genfisl (data, numMFs , mfType) ; 

[x,mf] = plotmf (fismat, 'input' ,1 ) ; 
subplot(2, 1 , 1 ) , plot(x,mf); 
xlabel( ' input 1 (pimf)'); 

[x,mf] = plotmf (fismat, 'input' ,2) ; 
subplot(2, 1 ,2) , plot(x,mf); 
xlabel( ' input 2 (trimf) 1 ); 



anf is 
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genfis2 


Purpose 

Synopsis 

Description 


Generate an FIS structure from data using subtractive clustering. 

fismat = genfis2(Xin, Xout, radii) 

fismat = genfis2(Xin , Xout , radii, xBounds) 

fismat = genfis2(Xin, Xout, radii, xBounds, options) 

Given separate sets of input and output data, genfis2 generates an FIS using 
fuzzy subtractive clustering. When there is only one output, genfis2 may be 
used togeneratean initial FIS for anfis training by first implementing 
subtractive clustering on the data, genf is2 accomplishes this by extracting a 
set of rules that models the data behavior. The rule extraction method first 
uses the subclust function to determine the number of rules and antecedent 
membership functions and then uses linear least squares estimation to 
determine each rule's consequent equations. This function returns an FIS 
structure that contai ns a set of fuzzy rules to cover the feature space. 

The arguments for genf is2 are as follows: 

• Xin is a matrix in which each row contains the input values of a data point. 

• Xout is a matrix in which each row contai ns the output values of a data point. 

• radii is a vector that specifies a cluster center's range of influence in each of 
the data dimensions, assuming the data falls within a unit hyperbox. For 
example, if the data dimension is 3 (e.g., Xin has two columns and Xout has 
one col umn), radii = [0.5 0.4 0.3] specifi es that the ranges of i nfl uence i n the 
first, second, and third data dimensions (i.e., the first column of Xin, the 
second column of Xin, and thecolumn of Xout) are 0.5, 0.4, and 0.3 times the 
width of the data space, respectively. If radii is a scalar, then the scalar 
value is applied to all data dimensions, i.e., each cluster center will have a 
spherical neighborhood of influence with the given radius. 

• xBounds is a 2-by -N optional matrix that specifies how to map the data in Xin 
and Xout into a unit hyperbox, where N is the data (row) dimension. The first 
row of xBounds contains the minimum axis range values and the second row 
contains the maximum axis range values for scaling the data in each 
dimension. For example, xBounds =[-10 0-1; 10 50 1] specifies that data 
values in the first data dimension are to be scaled from the range [-10 +10] 
into values in the range [0 1]; data values in the second data dimension are 
to be scaled from the range [0 50]; and data values in the third data 
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genfis2 


Examples 


See Also 


dimension are to be scaled from the range [-1 +1], If xBounds is an empty 
matrix or not provided, then xBounds defaults to the minimum and 
maximum data values found in each data dimension. 

• options is an optional vector for specifying algorithm parameterstooverride 
the default values. These parameters are explained in the help text for 
subclust on page 3-72. Default values are in placewhen this argument isnot 
specified. 

fismat = genfis2(Xin,Xout,0.5) 

This is the minimum number of arguments needed to use this function. Here a 
range of influence of 0.5 is specified for all data dimensions. 

fismat = genfis2(Xin,Xout, [0.5 0.25 0.3]) 

This assumes the combined data dimension is 3. Suppose Xin hastwocolumns 
and Xout hasonecolumn, then 0.5 and 0.25 are the ranges of influence for each 
of the Xin data dimensions, and 0.3 is the range of influence for the Xout data 
dimension. 

fismat = genfis2(Xin, Xout, 0.5, [ -10 -5 0; 10 5 20]) 

This specifies how to normalize the data in Xin and Xout into values in the 
range [0 l]for processing. Suppose Xin hastwocolumns and Xout has one 
column, then the data in the first column of Xin are scaled from [-10 +10], the 
data in the second column of Xin are scaled from [-5 +5], and the data in Xout 
are scaled from [0 20], 

subclust 
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gensurf 


Purpose 

Synopsis 

Description 


Generate an F I S output surface, 
gensurf (fis) 

gensurf (f is, inputs, out put) 

gensurf (fis, inputs, output, grids, ref input) 

gensurf (fis) generates a plot of the output surface of a given fuzzy inference 
system (fis) using the first two inputs and the first output. 

gensurf (fis, inputs, output) generates a plot using the inputs (one or two) 
and output (only one is allowed) given, respectively, by the vector, inputs, and 
the scalar, output. 

gensurf (fis, inputs, output .grids) allows you tospecify the number of grids 
in theX (first, horizontal) and Y (second, vertical) directions. If grids is a two 
element vector, the grids in theX and Y directions can be set independently. 

gensurf (fis, inputs, out put, grids, ref input) can be used if there are more 
than two outputs. The length of the vector ref input is the same as the number 
of inputs. 

• Enter NaNs for the entries of refinput corresponding tothe inputs whose 
surface is being displayed. 

• Enter real doublescalars to fix the values of other inputs. 

[x,y,z]=gensurf ( . . . ) returns the variables that define the output surface 
and suppresses automatic plotting. 
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gensurf 


Example 1 


Example 2 


See Also 


a = readfis( 'tipper' ) ; 
gensurf (a) 



a = gensurf (Temp, [1 2 ] , 1 , [ 20 20], [nan nan 0.2]); 

generates the surface of a three-i nput F I S named Temp from its fi rst two i nputs 
to its first output, while fixing a reference val ue for the third input at .2. 

evalfis, surfview 


3-41 


getfis 


Purpose 

Synopsis 


Description 


Get fuzzy system properties. 
getfis(a) 

getfis(a, 1 fisprop ' ) 

getfis(a, 1 vartype 1 , varindex, 1 varprop 1 ) 
getfis(a, 'vartype' , varindex, 'mf 1 ,mfindex) 
getfis (a, 1 vartype ' , varindex, 1 mf 1 , mf index, 'mfprop ' ) 

This is the fundamental access function for the FIS structure. With this one 
function you can learn about every part of the fuzzy i nference system. 

The arguments for getfis are as follows: 

• a: the name of a workspace variable FI S structure. 

• ' vartype 1 : a string indicating the type of variable you want (input or 
output). 

• varindex: an integer indicating the index of the variable you want (1, for 
input 1, for example). 

• 'mf 1 : a required string that indicates you are searching for membership 
function information. 

• mf index: the index of the membership function for which you are seeking 
information. 
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Examples 


getfis 


One input argument (output is the empty set) 

a = readfis( 'tipper' ) ; 
getfis (a) 

Name = tipper 
Type = mamdani 
Numlnputs = 2 
InLabels = 

service 

food 

NumOutputs = 1 
OutLabels = 
tip 

NumRules = 3 
AndMethod = min 
OrMethod = max 
ImpMethod = min 
AggMethod = max 
DefuzzMethod = centroid 

Two input arguments 

getfis (a, 1 type ' ) 

ans = 

mamdani 

Three input arguments (output is the empty set) 

getfis (a, 1 input ' , 1 ) 

Name = service 
NumMFs = 3 
MFLabels = 

poor 

good 

excellent 
Range = [0 10] 

Four input arguments 

getfis (a, 1 input ' , 1 , 1 name 1 ) 

ans = 

service 
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getfis 


See Also 


Five input arguments 

getfis (a, 1 input 1 , 1 , 1 mf 1 ,2) 

Name = good 
Type = gaussmf 
Params = 

1.5000 5.0000 

Six input arguments 

getfis (a, 1 input 1 , 1 , 1 mf 1 ,2, 1 name 1 ) 

ans = 

good 

setfis, showfis 
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mam2sug 


Purpose 

Synopsis 

Description 

Examples 


T ransform Mamdani FIS into a Sugeno FI S. 
sug_f is=mam2sug (mam_f is) 

mam2sug(mam_fis) transforms a (not necessarily single output) Mamdani FIS 
structure mam_f is intoa SugenoFIS structure sug_f is. The returned Sugeno 
system has constant output membership functions. These constants are 
determined by the centroids of the consequent membership functions of the 
original Mamdani system. The antecedent remains unchanged. 


mam_f ismat = readf is ( ' mam22 . f is 1 ) ; 
sug_fismat = mam2sug (mam_f ismat ) ; 




subplot(2,2, 1 ) ; gensurf (mam_f ismat , 
title ( 'Mamdani system (Output 1)'); 

[1 

2] . 

, i); 

subplot(2,2,2) ; gensurf (sug_f ismat , 
title( 'Sugeno system (Output 1)'); 

[1 

2] , 

, i); 

subplot(2,2,3) ; gensurf (mam_f ismat , 
title( 'Mamdani system (Output 2)'); 

[1 

2] , 

, 2); 

subplot(2,2,4) ; gensurf (sug_f ismat , 
title( 'Sugeno system (Output 2)'); 

[1 

2] , 

, 2); 
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mf2mf 


Purpose 

Synopsis 

Description 


Examples 


See Also 


Translates parameters between membership functions. 
outParams = mf 2mf (inParams , inType , outType) 

This function translates any built-in membership function type into another, 
in terms of its parameter set. I n principle, mf2mf mimics the symmetry points 
for both the new and old membership functions. Occasionally this translation 
results in lost information, sothat if theoutput parameters aretranslated back 
into the original membership function type, the transformed membership 
function will not look the same as it did originally. 

The input arguments for mf2mf are as follows: 

• inParams: the parameters of the membership function you are transforming 

• inType: a string name for the type of membership function you are 
transforming 

• outType: a string name for the new membership function you are 
transforming to 

x=0:0. 1:5; 
mf pi = [1 23]; 

mfp2 = mf2mf (mf pi , 1 gbellmf ' , ' trimf ' ) ; 
plot(x,gbellmf (x,mfp1 ) ,x, trimf (x,mfp2) ) 



dsigmf , gaussmf , gauss2mf, gbellmf, evalmf, pimf , psigmf, sigmf, smf, 
trapmf , trimf, zmf 
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mfedit 


Purpose 

Synopsis 

Description 


Membership function editor. 

mfedit ( 1 a 1 ) 
mfedit (a) 
mfedit 



mfedit ( 1 a 1 ) generates a membership function editor that allows you to 
inspect and modify all the membership functions for your FIS stored in thefile, 
a.f is. 

mfedit(a) operates on a MATLAB workspace variable, for a FIS structure, a. 

mfedit alone brings up the membership function editor with no FI S loaded. 

For each membership function you can change the name, the type, and the 
parameters. Eleven built-in membership functions are provided for you to 
choose from, although of course you can always create your own specialized 
versions. Refer to 'The Membership Function Editor" on page 2-52 for more 
information about howto use mfedit. 
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mfedit 


Menu Items 


Membership 
Function 
Pop-up Menu 

See Also 


Select the icon for thevariableon the upper left side of the diagram (under the 
heading "FIS Variables") to display its associated membership functions in the 
plot region. Select membership functions by clicking once on them or their 
labels. 

On the Membership F unction Editor, there is a menu bar that allows you to 
open related GUI tools, open and save systems, and so on. The File menu for 
the Membership Function Editor is the same as the one found on the FI S 
Editor. Refer to fuzzy on page 3-29 for more information. 

• Under Edit, select: 

Add MF... to add membership functions to the current variable. 

Add custom MF... to add a customized membership function to the current 
variable. 

Remove current MF to delete the current membership function. 

Remove all MFsto delete all membership functions of the current variable. 
Undo to undo the most recent change. 

• Under View, select: 

Edit FIS properties... to invokethe FIS Editor. 

Edit rules... to invokethe Rule Editor. 

View rules... to invokethe Rule Viewer. 

View surface... to invoke the Surface Viewer. 

There are 11 built-in membership functions to choose from, and you also have 
the option of installing a customized membership function. 


fuzzy, ruleedit, ruleview, surfview 
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newfis 


Purpose 

Synopsis 


Description 


Examples 


See Also 


Create new FIS. 

a=newf is (f isName , f isType , andMethod , orMethod , impMethod , . . . 

aggMethod , def uzzMethod ) 

Thisfunction creates new FIS structures, newfis has upto seven input 
arguments, and the output argument is an FIS structure. The seven input 
arguments are as follows: 

• f isName is the string name of the FIS structure, f isName. f is you create. 

• f isType is the type of FIS. 

• andMethod, orMethod, impMethod, aggMethod, and def uzzMethod, 
respectively providethe methods for AND, OR, implication, aggregation, and 
defuzzification. 

The following example shows what the defaults are for each of the methods: 

a=newf is ( 1 newsys 1 ) ; 
getf is (a) 

returns 

Name = newsys 
Type = mamdani 
Numlnputs = 0 
InLabels = 

NumOutputs = 0 
OutLabels = 

NumRules 0 
AndMethod min 
OrMethod max 
ImpMethod min 
AggMethod max 
DefuzzMethod centroid 

ans = 

[newsys] 
readf is, writef is 
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parsrule 


Purpose 

Synopsis 

Description 


Examples 

See Also 


Parse fuzzy rules. 

f is2 = parsrule(fis,txtRuleList) 

fis2 = parsrule (fis , txtRuleList , ruleFormat ) 

f is2 = parsrule(fis,txtRuleList,ruleFormat,lang) 

This function parses the text that defines the rules (txtRuleList) for a 
MATLAB workspaceFIS variable, fis, and returns a FIS structure with the 
appropriate rule list in place. If the original input FIS structure, fis, has any 
rules initially, they are replaced in the new structure, fis2. Three different 
ruleformats (indicated by ruleFormat) aresupported: 'verbose 1 , "symbolic," 
and "indexed. "The default format is "verbose". When the optional language 
argument, lang, is used, the rules are parsed in verbose mode, assuming the 
key words are in the language, lang. This language must be either 'english', 
'francais', or 'deutsch'.Thekey language words in English are: if, then, is, AND, 
OR, and NOT. 

a = readfis( 'tipper' ) ; 

ruleTxt = 'if service is poor then tip is generous 1 ; 
a2 = parsrule (a, ruleTxt ,' verbose 1 ) ; 
showrule (a2) 
ans = 

1. If (service is poor) then (tip is generous) (1) 
addrule, ruleedit, showrule 
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pimf 


Purpose 

Synopsis 

Description 


Examples 


n -shaped built-in membership function, 
y = pimf (x, [a b c d] ) 

This spline-based curve is so named because of its n -shape. This membership 
function isevaluated atthe points determined by the vector x. The parameters 
a and d locate the "feet" of the curve, while b and c locate its "shoulders." 

x=0:0. 1:10; 

y=pimf (x, [1 4 5 10] ) ; 

plot (x,y) 

xlabel( 'pimf , P=[1 4 5 10]') 



See Also 


dsigmf, gaussmf, gauss2mf , gbellmf , evalmf, mf2mf, psigmf , sigmf , smf, 
trapmf , trimf , zmf 
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plotfis 


Purpose 

Synopsis 

Description 


Examples 


See Also 


Plot an FIS. 

plotfis(fismat) 

This function displays a high level diagram of an FIS, f ismat. I nputs and their 
membership functions appear to the left of the FIS structural characteristics, 
while outputs and their membership functions appear on the right. 

a = readfis( 'tipper 1 ) ; 
plotfis(a) 



evalmf , plotmf 
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plotmf 


Purpose 

Synopsis 

Description 


Examples 


Plot all of the membership functions for a given variable, 
plotmf (fismat , varType , var Index) 

This function plots all of the membership functions in the FI S called fismat 
associated with a given variable whose type and index are respectively given 
by varType ('input' or 'output'), and varlndex. This function can also be 
used with the M ATLAB function, subplot. 

a = readfis( 'tipper' ) ; 
plotmf (a, ' input 1 , 1 ) 



service 


See Also 


evalmf , plotf is 
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psigmf 


Purpose 

Synopsis 

Description 


Examples 


See Also 


Built-in membership function composed of the product of two 
sigmoidally-shaped membership functions. 


y = psigmf (x,[a1 cl a2 c2]) 


The sigmoid curve plotted for the vector x depends on two parameters a and c 
as given by 


f(x;a, c) 


1 + e 


a(x-c) 


psigmf issimplythe product of two such curves plotted for the values of the 
vector x 

fj(x; al, Cj) * f 2 (x; a2, c 2 ) 

The parameters are listed in the order: [al c 2 a 2 c 2 \. 
x=0:0. 1:10; 

y=psigmf (x, [23 -58]); 
plot(x,y) 

xlabel( 1 psigmf , P=[2 3 -58]') 



dsigmf , gaussmf , gauss2mf , gbellmf, evalmf , mf2mf , pimf , sigmf , smf, trapmf , 
trimf, zmf 
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readfis 


Purpose 

Synopsis 

Description 

Examples 


See Also 


Load an FIS from the disk, 
fismat = readfis (' filename ' ) 

Read a fuzzy i nference system from a .fis file (named filename) on the disk 
and bri ng the resulti ng fi le i nto the workspace. 

fismat = readfis (no input arguments) brings up a uigetfile dialog box to 
assist with the name and directory location of the file. 

fismat = readfis( 'tipper' ) ; 
getfis(fismat) 

returns 

Name = tipper 
Type = mamdani 
Numlnputs = 2 
InLabels = 
service 
food 

NumOutputs = 1 
OutLabels = 
tip 

NumRules = 3 
AndMethod = min 
OrMethod = max 
ImpMethod = min 
AggMethod = max 
DefuzzMethod = centroid 
ans = 
tipper 


writef is 
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rmmf 


Purpose 

Synopsis 

Description 


Examples 


See Also 


To remove membership function from an FIS. 

fis = rmmf (f is , 1 varType varlndex, ' mf 1 , mf Index) 

fis = rmmf (fis, varType, varlndex, 'mf ' ,mf Index) removes the membership 
function, mf Index, of variable type, varType, of index varlndex, from the fuzzy 
inference system associated with the workspace FIS structure, fis: 

• The string varType must be 'input 1 or 'output 1 . 

• varlndex is an integer for the index of the variable. This index represents the 
order in which the variables are listed. 

•Theargument 'mf' is a string representing the membership function. 

• mf index is an integer for the index of the membership function. This index 
represents the order in which the membership functions are listed. 

a = newfis( 'mysys ' ) ; 

a = addvar (a, ' input ',' temperature ', [0 100]); 
a = addmf (a, 'input' ,1 , 'cold' , 'trimf' , [0 30 60]); 
getf is (a, ' input ' , 1 ) 

returns 


Name = temperature 
NumMFs = 1 
MFLabels = 

cold 

Range = [0 100] 
b = rmmf (a, ' input ', 1 ,' mf ', 1 ) ; 
getfis(b, ' input ' , 1 ) 

returns 


Name = temperature 
NumMFs = 0 
MFLabels = 

Range = [0 100] 

addmf, addrule, addvar, plotmf , rmvar 
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rmvar 


Purpose 

Synopsis 

Description 


Examples 


To remove variables from an FIS. 

[f is2, errorStr] = rmvar(fis, 1 varType ' , varlndex) 
fis2 = rmvar(fis, ' varType varlndex) 

f is2 = rmvar (f is , ' varType ' , varlndex) removes the variable 1 varType ' , of 
index varlndex, from the fuzzy i nference system associated with the 
workspace FIS structure, fis: 

• The string varType must be 'input 1 or 'output 1 . 

• varlndex is an integer for the index of the variable. This index represents the 
order in which the variables are listed. 

[fis2, errorStr] = rmvar(fis, ' varType varlndex) returns any error 
messages to the string, errorStr. 

This command automatically alters the rule list to keep its size consistent with 
the current number of variables. You must delete from the FIS any rule that 
contains a variableyou want to remove, before removing it. You cannot remove 
a fuzzy variable currently in use in the rule list. 

a = newfis( 'mysys 1 ) ; 

a = addvar (a, ' input 1 , 1 temperature ', [0 100]); 
getfis(a) 
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rmvar 


See Also 


returns 

Name = mysys 

Type = mamdani 

Numlnputs = 1 
InLabels = 

temperature 
NumOutputs = 0 
OutLabels = 

NumRules = 0 
AndMethod = min 
OrMethod = max 
ImpMethod = min 
AggMethod = max 
Def uzzMethod = centroid 

ans = 
mysys 

b = rmvar (a, ' input 1 , 1 ) ; 
getfis(b) 
returns 

Name = mysys 
Type = mamdani 
Numlnputs = 0 
InLabels = 

NumOutputs = 0 
OutLabels = 

NumRules = 0 
AndMethod = min 
OrMethod = max 
ImpMethod = min 
AggMethod = max 
DefuzzMethod = centroid 

ans = 
mysys 


addmf, addrule, addvar, rmmf 
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ruleedit 


Purpose 

Synopsis 

Description 


Rule editor and parser. 

ruleedit ( 1 a 1 ) 
ruleedit (a) 



Language and Format 
are options on this 
GUI. Languages are 
English, Deutsch, and 
Francois. Formats are 
verbose, symbolic, 
and indexed. 


The Rule Editor, when invoked using ruleedit ( 1 a 1 ) , is used to modify the 
rules of an FIS structure stored in a file, a.fis. It can also be used to inspect 
the rules being used by a fuzzy inference system. 

To use this editor to create rules, you must first define all of the input and 
output variables you want to use with the FIS editor. You can create the rules 
using the listbox and check box choices for input and output variables, 
connections, and weights. Refer to 'The Rule Editor" on page 2-56 for more 
information about how to use ruleedit . 

The syntax ruleedit (a) is used when you want to operate on a workspace 
variable for an FIS structure called a. 
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ruleedit 


Menu Items 


See Also 


On the Rule Editor, there is a menu bar that allows you to open related GUI 
tools, open and save systems, and soon. The File menu for the Rule Editor is 
the same as the one found on the FI S Editor. Refer to fuzzy on page 3-29 for 
more information. 

• Usethe following Edit menu item: 

Undo to undo the most recent change. 

• Use the following View menu items: 

Edit FIS properties... to invoke the FI S Editor. 

Edit membership functions... to invoke the Membership F unction 
Editor. 

View rules... to invoke the Rule Viewer. 

View surface... to i nvoke the Surface Viewer. 

• Usethe Options menu items: 

Language to select the language: English, Deutsch, and Francais 

Format to select the format: 

* verbose uses the words "if," "then," "AN D," "OR," and so on to create 
actual sentences. 

' symbolic substitutes some symbols for the words used in the verbose 
mode. For example, "if A AND 6 then C" becomes "A & B ->C." 

' indexed mirrors how the rule is stored in the FIS structure. 

addrule, fuzzy, mfedit, parsrule, ruleview, showrule, surfview 
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ruleview 


Purpose 

Synopsis 

Description 


Menu Items 


Rule viewer and fuzzy inference diagram. 
ruleview( 1 a 1 ) 



The Rule Viewer invoked using ruleview( 'a' ) depicts the fuzzy inference 
diagram for an FI S stored in a file, a. f is. It is used to view the entire 
implication process from beginning to end. You can move around the line 
indices that correspond to the inputs and then watch the system readjust and 
compute the new output. Refer to ‘The Rule Viewer" on page 2-59 for more 
information about how to use ruleview. 

On the Rule Viewer, there is a menu bar that allows you to open related GUI 
tools, open and save systems, and soon. The File menu for the Rule Viewer is 
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See Also 


the same as the one found on the FI S Editor. Refer to fuzzy on page 3-29 for 
more information. 

• Use the View menu items: 

Edit FIS properties... to invoke the FI S Editor 

Edit membership functions... to invoke the Membership F unction Editor 
Edit rules... to invoke the Rule Editor 
View surface... to invoke the Surface Viewer 

• Use the Options menu item: 

Rule display format to set the format in which the rule appears. If you 
click on the rule numbers on the left side of the fuzzy inference diagram, 
the rule associated with that number will appear in the Status Bar at the 
bottom of the Rule Viewer. 

fuzzy, mfedit, ruleedit, surfview 
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Purpose 

Synopsis 

Description 


Set fuzzy system properties, 
a = setfis (a ,' fispropname' ,' newfisprop' ) 

a = setf is( a, 1 vartype ' ,varindex, 1 varpropname ' , 'newvarprop ' ) 
a = setf is(a, 1 vartype ' ,varindex, 1 mf 1 ,mf index, ... 

1 mfpropname 1 , 1 newmf prop 1 ) ; 

The command setfis can be called with three, five, or seven input arguments, 
depending on whether you want to set a property of the entire FI S structure, 
for a particular variable belonging to that FIS structure, or for a particular 
membership function belonging to one of those variables. The arguments are: 

• a - a variable name of an FI S from the workspace 

• 'vartype' - a string indicating the variabletype: input or output 

• varindex - the index of the input or output variable 

• 1 mf 1 - a requi red stri ng for the fourth argument of a 7-argument cal I for 
setfis, indicating this variable is a membership function 

• mf index - the index of the membership function belonging to the chosen 
variable 

• 'fispropname' - a string indicating the property of the F I S field you want 
to set: name, type, andmethod, ormethod, impmethod, aggmethod, 
defuzzmethod 

• 'newfisprop' - a string describing the name of the FI S property or method 
you want to set 

• 'varpropname' - a string indicatingthenameof thevariablefield you want 
to set: name or range 

• 'newvarprop' - astringdescribingthenameofthevariableyou wanttoset 
(for name), or an array describing the range of that variable (for range) 

• 'mfpropname' - a string indicating the name of the membership function 
field you want to set: name, type, or params . 

• 'newmfprop' - a string describing the name or type of the membership 
function field want to set (for name or type), or an array describing the range 
of the parameters (for params) 
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Examples Called with three arguments, 

a = readfis( 'tipper' ) ; 

a2 = setfis(a, 'name 1 , 'eating'); 

getfis(a2, 'name'); 

Results in 

out = 
eating 

If used with five arguments, setfis will update two variable properties. 

a2 = setfis(a, 1 input 1 , 1 , 1 name 1 , 1 help 1 ) ; 
getfis(a2, ' input 1 , 1 , ' name ' ) 
ans = 
help 

If used with seven arguments, setfis will update any of several membership 
function properties. 

a2 = setfis(a, 1 input 1 , 1 , 1 mf 1 ,2, 1 name 1 , 'wretched 1 ) ; 
getf is(a2, ' input ' , 1 , ' mf ' , 2, ' name ' ) 
ans = 

wretched 

See Also getfis 
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Purpose 

Synopsis 

Description 


See Also 


Fuzzy inference S-functi on for Simulink. 
output = sffis(t,x,u,flag,fismat) 

This MEX-file is used by Simulink to undertake the calculation normally 
performed by evalfis. It has been optimized to work in theSimulink 
environment. This means, among other things, that sffis builds a data 
structure in memory during the initialization phase of a Simulink simulation, 
which it then continues to use until the simulation is complete. 

The arguments t, x, and flag are standard Simulink S-function arguments 
(see Chapter 8, "S-F unctions" in the Using Simulink documentation). The 
argument u is the input to the M ATLAB workspace FI S structure, f ismat. If, 
for example, there are two inputs to f ismat, then u will be a two-element 
vector. 

evalfis, fuzblock 
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Purpose 

Synopsis 

Description 


Examples 


Display annotated FIS. 
showfis(fismat) 

showfis(fismat) prints a version of the MATLAB workspace variable FIS, 
f ismat , allowing you to see the significance and contents of each field of the 
structure. 


a = readfis( 'tipper' ) ; 
showfis(a) 


returns 


1 . 

Name 

tipper 

2. 

Type 

mamdani 

3. 

Inputs/Outputs 

[2 1] 

4. 

NumlnputMFs 

[3 2] 

5. 

NumOutputMFs 

3 

6. 

NumRules 

3 

7. 

AndMethod 

min 

8. 

OrMethod 

max 

9. 

ImpMethod 

min 

10. 

AggMethod 

max 

11 . 

Def uzzMethod 

centroid 

12. 

InLabels 

service 

13. 


food 

14. 

OutLabels 

tip 

15. 

InRange 

[0 10] 

16. 


[0 10] 

17. 

OutRange 

[0 30] 

18. 

InMFLabels 

poor 

19. 


good 

20. 


excellent 

21 . 


rancid 

22. 


delicious 

23. 

OutMFLabels 

cheap 

24. 


average 

25. 


generous 

26. 

InMFTypes 

gaussmf 

27. 


gaussmf 
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See Also 


28. 


gaussmf 

29. 


trapmf 

30. 


trapmf 

31 . 

OutMFTypes 

trimf 

32. 


trimf 

33. 


trimf 

34. 

InMFParams 

[1.5 0 0 0] 

35. 


[1.5 5 0 0] 

36. 


[1.5 10 0 0] 

37. 


[0013] 

38. 


[7 9 10 10] 

39. 

OutMFParams 

[0 5 10 0] 

40. 


[10 15 20 0] 

41 . 


[20 25 30 0] 

42. 

Rule Antecedent 

[1 1] 

43. 


[2 0] 

44. 


[3 2] 

42. 

Rule Consequent 

1 

43. 


2 

44. 


3 

42. 

Rule Weigth 

1 

43. 


1 

44. 


1 

42. 

Rule Connection 

2 

43. 


1 

44. 


2 


getf is 
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Purpose 

Synopsis 


Description 


Display FIS rules. 

showrule (f is) 
showrule (f is , indexList ) 
showrule (f is , indexList , format) 
showrule (f is , indexList , format , Lang) 

This command is used to display the rules associated with a given system. It 
can be invoked with one to four arguments. The first argument, fis, is 
required. This is the M ATLAB workspace variable name for a FI S structure. 
The second (optional) argument indexList is the vector of rules you want to 
display. Thethird argument (optional) isthestring representing the format in 
which the rules are returned, showrule can return the rule in any of three 
different formats: 'verbose 1 (the default mode, for which English is the default 
language), 'symbolic ' , and 'indexed' , for membership function index 
refer end ng. 

When used with four arguments, the forth argument must be verbose, and 
showrule(fis, indexList .format, lang) displays the rules in the language 
given by lang, which must be either 'english', 'francais', or 'deutsch'. 
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Examples 


See Also 


a = readfis( 'tipper' ) ; 
showrule(a, 1 ) 
ans = 

1. If (service is poor) or (food is rancid) then (tip is cheap) (1) 

showrule(a,2) 
ans = 

2. If (service is good) then (tip is average) (1) 

showrule(a, [3 1 ],' symbolic ' ) 
ans = 

3. (service==excellent) | (food==delicious) => (tip=generous) (1) 
1. (service==poor) | (food==rancid) => (tip=cheap) (1) 

showrule (a, 1 :3, ' indexed ' ) 
ans = 

1 1, 1 (1) : 2 
2 0, 2 (1) : 1 
3 2, 3 (1) : 2 

parsrule, ruleedit, addrule 
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Purpose 

Synopsis 

Description 


Examples 


See Also 


Sigmoidally-shaped built-in membership function. 


y = sigmf (x, [a c] ) 


The sigmoidal function, sigmf (x, [a c] ), as given below by f(x, a, c), is a 
mapping on a vector x, and depends on two parameters a and c: 


f(x , a, c) = 


1 + e 


a(x-c ) 


Depending on the sign of the parameter a, the sigmoidal membership function 
is inherently open to the right or to the left, and thus is appropriate for 
representing concepts such as "very large" or "very negative." M ore 
conventional-looking membership functions can be built by taking either the 
product or difference of two different sigmoidal membership functions. You can 
find out more about this in this chapter's entries for dsigmf and psigmf . 

x=0:0. 1:10; 
y=sigmf (x, [2 4] ) ; 
plot (x,y) 

xlabel( 'sigmf , P=[2 4]') 



dsigmf, gaussmf , gauss2mf, gbellmf , evalmf, mf2mf, pimf , psigmf, smf, 
trapmf , trimf , zmf 
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Purpose 

Synopsis 

Description 


Examples 


S-shaped built-in membership function, 
y = smf (x, [a b] ) 

This spline-based curve is a mapping on the vector x, and is named because of 
its S-shape. The parameters a and b locate the extremes of the sloped portion 
of the curve. 

x=0:0. 1:10; 
y=smf (x, [ 1 8]); 
plot (x,y) 

xlabel( 'smf, P=[1 8] 1 ) 



See Also 


dsigmf, gaussmf, gauss2mf , gbellmf , evalmf, mf2mf, pimf , psigmf , sigmf, 
trapmf , trimf , zmf 
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subclust 


Purpose 

Synopsis 

Description 


Find cluster centers with subtractive clustering. 

[C,S] = subclust(X, radii, xBounds, options) 

This function estimates the cluster centers in a set of data by using the 
subtractive clustering method. The subtractive clustering method assumes 
each data point is a potential cluster center and calculates a measure of the 
likelihood that each data point would define the cluster center, based on the 
density of surrounding data points. The algorithm 

• Selects the data point with the highest potential tobe the first cluster center 

• Removes all data points in the vicinity of the first cluster center (as 
determined by radii), in order to determine the next data cluster and its 
center location 

• Iterates on this process until all of the data is within radii of a cluster center 

The subtractive clustering method is an extension of the mountain clustering 
method proposed by R. Yager. 

The matrix x contains the data to be clustered; each row of x is a data point. 
The variable radii is a vector of entries between 0 and lthat specifies a cluster 
center's range of influence in each of the data dimensions, assuming the data 
falls within a unit hyperbox. Small radii values generally result in finding a 
few large clusters. Good values for radii are usually between 0.2 and 0.5. 

For example, if the data dimension is two (x has two columns), 
radii =[0.5 0.25] specifies that the range of influence in the first data 
dimension is half the width of the data space and the range of influence in the 
second data dimension is one quarter the width of the data space. If radii is a 
scalar, then the scalar value is applied toall data dimensions, i.e., each cluster 
center will have a spherical neighborhood of influence with the given radius. 

xBounds is a 2-by-N matrix that specifies how to map the data in x into a unit 
hyperbox, where N is the data dimension. This argument is optional if x is 
already normalized. Thefirst row contains the mini mum axis range values and 
the second row contains the maxi mum axis range values for scaling the data in 
each dimension. For example, xBounds =[-10-5; 10 5] specifiesthat data values 
in thefirst data dimension are to be scaled from the range [-10 +10] into values 
in the range [0 1]; data values in the second data dimension are to be scaled 
from the range [-5 +5] into values in the range [0 1], If xBounds is an empty 
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Examples 


matrix or not provided, then xBounds defaults to the minimum and maximum 
data values found in each data dimension. 

The options vector can be used for specifyi ng cl usteri ng algorithm parameters 
to override the default values. These components of the vector options are 
specified as follows: 

•options(l) = quashFactor: This is the factor used to multiply the radii 
values that determine the neighborhood of a cluster center, so as to quash 
the potential for outlying points to be considered as part of that cluster, 
(default: 1.25) 

• options(2) = acceptRatio: This sets the potential, as a fraction of the 
potential of the first cluster center, above which another data point will be 
accepted as a cluster center, (default: 0.5) 

• options(3) = rejectRatio: This sets the potential, as a fraction of the 
potential of thefirst cluster center, below which a data point will berejected 
as a cluster center, (default: 0.15) 

• options(4) = verbose: If this term is not zero, then progress information 
will be printed as the clustering process proceeds, (default: 0) 

The function returns the cluster centers in the matrix C; each row of C contains 
the position of a cluster center. The returned S vector contains the sigma values 
that specify the range of influence of a cluster center in each of the data 
dimensions. All cluster centers share the same set of sigma values. 

[C,S] = subclust(X,0.5) 

This is the mini mum number of arguments needed to usethisfunction. A range 
of influence of 0.5 has been specified for all data dimensions. 

[C,S] = Subclust(X, [0.5 0.25 0.3] , [] , [2.0 0.8 0.7]) 

This assumes the data dimension is 3 (x has 3 columns) and uses a range of 
influence of 0.5, 0.25, and 0.3 for thefirst, second and third data dimension, 
respectively. The scaling factors for mapping the data into a unit hyperbox will 
be obtained from the minimum and maximum data values. The squashFactor 
is set to 2.0, indicating that we only want to find clusters that are far from each 
other. The acceptRatio is set to 0.8, indicating that we will only accept data 
points that have a very strong potential for being cluster centers. The 
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See Also 
References 


re ject Ratio is set to 0.7, indicating that we want to reject all data points 
without a strong potential. 

gent is2 

Chiu, S., "Fuzzy Model Identification Based on Cluster Estimation, '7 ournal of 
I intelligent & Fuzzy Systems, Vol. 2, No. 3, Sept. 1994. 

Yager, R. and D. Filev, "Generation of Fuzzy Rules by Mountain Clustering," 
J ournal of Intelligent & Fuzzy Systems, Vol. 2, No. 3, pp. 209-219, 1994. 
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Purpose 

Synopsis 

Description 


Output surface viewer, 
surf view( 1 a 1 ) 



The Surface Viewer invoked using surf view ( 1 a 1 ) is a GUI tool that lets you 
examine the output surface of an FIS stored in a file, a. f is, for any one or two 
inputs. Since it does not alter the fuzzy system or its associated FIS structure 
in anyway, it is a read-only editor. Using the pop-up menus, you select the two 
input variables you want assigned tothe two input axes (X and Y), as well the 
output variable you want assigned to the output (or Z) axis. Select the 
Evaluate button to perform the calculation and plot the output surface. 

By clicking on the plot axes and dragging the mouse, you can manipulate the 
surface so that you can view it from different angles. 

I f there are more than two i nputs to your system, you must suppl y the constant 
values associated with any unspecified inputs in the reference input section. 

Refer to 'The Surface Viewer" on page 2-61 for more information about howto 
use surfview. 
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Menu Items 


See Also 


On the Surface Vi ewer, there is a menu bar that allows you to open related GUI 
tools, open and save systems, and soon. The File menu for the Surface Viewer 
is the same as the one found on theFIS Editor. Refer tofuzzy on page 3-29 for 
more information. 

• Use the View menu items: 

Edit FIS properties... to invokethe FIS Editor. 

Edit membership functions... to invokethe Membership F unction Editor. 
Edit rules... to invokethe Rule Editor. 

View rules... to invokethe Rule Viewer. 

• Use the Options menu items: 

Plot to choose among eight different kinds of plot styles. 

Color Map to choose among several different color schemes. 

Always evaluate to automatically evaluate and plot a new surface every 
time you make a change that affects the plot (I ike changing the number of 
grid points). This is the default option. To deselect this option, select it once 
more. 

anfisedit, fuzzy, gensurf, mfedit, ruleedit, ruleview 
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Purpose 

Synopsis 

Description 


Examples 


Trapezoidal -shaped built-in membership function, 
y = trapmf (x, [a b c d] ) 

The trapezoidal curve is a function of a vector, x, and depends on four scalar 
parameters a, b, c, and d, as given by 


f(x;a, b , c, d) 


0 , 

x-a 

b~a' 

1 , 


x< a 
a<x<b 
b<x<c ■ 


d-x 

d~c' 


c<x<d 


0, d<x ) 


or, more compactly, by 

f(x-,a, b , c, d) = max^min^^— 1 , j oj 

The parameters a and d locate the "feet" of the trapezoid and the parameters/? 
and c locate the "shoulders." 

x=0:0. 1:10; 

y=trapmf (x, [1 5 7 8] ) ; 
plot(x,y) 

xlabel( 'trapmf , P=[1 5 7 8]') 
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See Also 


dsigmf , gaussmf , gauss2mf , gbellmf, evalmf , mf2mf , pimf , psigmf , sigmf , smf, 
trimf, zmf 
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Purpose 

Synopsis 

Description 


Examples 


Triangular-shaped built-in membership function. 


y = trimf (x,params) 
y = trimf (x, [a b c] ) 


The triangular curve is a function of a vector, x, and depends on three scalar 
parameters a, b, and c, as given by 


f(x-a , b , c) = 


0, x< a 
x-a 


, a<x<b 


b-a 
c-x , 

r - B , b<x<c 
0, c<x 


or, more compactly, by 

f(x-,a,b,c ) = max^min^^— °) 


The parameters a and c locate the "feet" of the triangle and the parameter c 
locates the peak. 

x=0:0. 1:10; 
y=trimf (x, [3 6 8] ) ; 
plot(x,y) 

xlabel( 'trimf , P=[3 68]') 
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See Also 


dsigmf , gaussmf , gauss2mf , gbellmf, evalmf , mf2mf , pimf , psigmf , sigmf , smf, 
trapmf 
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Purpose 

Synopsis 


Description 


Examples 


See Also 


Save an F I S to the disk. 

writefis (fismat ) 

writefis (fismat , ' filename ' ) 

writefis (fismat , 1 filename ' , ' dialog ' 

writefis saves a MATLAB workspace FI S structure, fismat, as a . f is file on 
disk. 

writefis (fismat ) brings up a dialog box to assist with the naming and 
directory location of the file. 

writefis (fismat , 1 filename ' ) writes a .f is file corresponding to the FI S 
structure, fismat , to a disk file cal led f ilename .f is. No dialog box is used and 
the file is saved to the current directory. 

writefis (fismat , 1 filename ' , 1 dialog 1 ) brings up a dialog box with the 
default name f ilename. f is supplied. 

The extension . f is is only added to filename if it is not already included in the 
name. 

a = newfis( 'tipper' ) ; 

a = addvar (a, ' input 1 , 1 service ', [0 10]); 

a = addmf (a, 'input' ,1 , 'poor' , 'gaussmf ' , [1 .5 0]); 

a = addmf (a, 'input' ,1 , 'good' , 'gaussmf ', [1 .5 5]); 

a = addmf (a, ' input ', 1 ,' excellent 1 , 1 gaussmf ',[ 1 .5 10]); 

writefis(a, 'my_file' ) 

readf is 
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Purpose 

Synopsis 

Description 

Examples 


Z-shaped built-in membership function, 
y = zmf (x, [a b] ) 

This spline-based function of x is so named because of its Z-shape. The 
parameters a and b locate the extremes of the sloped portion of the curve. 

x=0:0. 1:10; 
y=zmf (x, [3 7] ) ; 
plot(x,y) 

xlabel( 1 zmf, P=[3 7] 1 ) 



See Also 


dsigmf , gaussmf , gauss2mf , gbellmf, evalmf , mf2mf , pimf , psigmf , sigmf , smf, 
trapmf , trimf 
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addmf 72, 76, 3-6 
addrule 72, 76, 3-8 
addvar 76, 3-9 

aggregation 36, 40, 44, 60, 63, 86, 132, 133 
AND 30, 37, 63, 103 
ANFIS 93 

antis 7, 92, 95, 104, 109, 112, 114, 125, 3-10 
options 116 

ANFIS Editor GUI 7, 45, 92, 95, 106, 109, 114, 
3-14 

anf isedit 67, 93, 95, 3-14 
antecedent 35, 37, 59, 74, 86, 88, 132 


B 

backpropagation 104 


c 

chaotic time series 120 
checking data 94, 118 
checking error 119 
clustering 120, 121, 128, 132 
clustering algorithms 133 
clustering GUI 128 
consequent 32, 35, 40, 59, 74, 86, 132 
antecedent 37 
convertfis 3-16 

D 

defuzz 3-17 

defuzzification 32, 36, 40, 49, 63, 95, 132 
defuzzify 32, 41 

degree of membership 20, 24, 33, 35, 38 
distfcm 121 


dsigmf 27, 3-18 


E 

error tolerance 104 
evalf is 73, 126, 3-19, 3-26 
evalmf 3-21 


F 

f cm (fuzzy c-means) 120, 128, 3-22 
findcluster 128, 3-24 

F I S 43, 47, 49, 53, 58, 65, 67, 73, 84, 92, 95, 109, 
125, 132, 3-10 
C-code 130 
Editor 45, 49, 68, 115 
files 76 

generating 100 
Mamdani-type 49, 86, 91 
matrix 73 
saving a FIS 62 
structure 93, 115, 117 
Sugeno-type 86, 88, 95 
Sugeno-type See also Sugeno-type inference 
97 

fuzblock 83, 3-27 
f uzdemos 3-28 

fuzzification 32, 37, 42, 49, 132 
fuzzy 3-29 

fuzzy clustering 115, 120 
fuzzy c-means clustering 3-22 
F uzzy I nference System (F I S) 19, 45, 132 
fuzzy operators 29, 32, 37, 39, 44, 58, 63, 71, 131 
fuzzy set 20, 23, 26, 35, 38, 40, 60, 86, 88, 132 
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gauss2mf 27, 3-32 
gaussian 27 
gaussmf 27, 3-34 
gbellmf 27, 3-35 
genfis 111 
genfisl 101, 3-36 
genf is2 101, 126, 3-38 
gensurf 70, 88, 90, 3-40 
getf is 66, 76, 89, 3-42 
glossary 132 
grid partition 100 


H 

hybrid method 104 


I 

if-then rules 32 
antecedent 32 
consequent 32 

implication 32, 35, 37, 39, 43, 60, 63, 86 
implication Seealso if-then rules 32, 133 
initfcm 121 


L 

logical operations 28 


M 

mam2sug 3-45 
Mamdani's method 36 
M amdani -style inference 133 
M amdani -type inference 35, 49, 86, 90 
max 41, 43 


membership function 23, 28, 35, 53, 55, 118, 133 
mf editor 103 

Membership Function Editor 45, 47, 52 
membership functions 
bell 27 
custom 63 
Gaussian 27 
Pi 27 
S 27 

sigmoidal 27 
Z 27 

MF Seealso membership function 
mf2mf 3-46 
mfedit 3-47 
min 43 

model validation 94, 98 
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neuro-fuzzy inference 93 
newt is 72, 3-49 
NOT 30, 103 


o 

OR 30, 37, 43, 103 


P 

parsrule 3-50 
pimf 28, 3-51 
plotf is 68, 3-52 
plotmf 68, 90, 3-53 
probabilistic OR 39 
probor 41 
psigmf 27, 3-54 
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readf is 65, 73, 88, 3-55 
rmmf 76, 3-56 
rmvar 76, 3-57 
Rule Editor 45, 56 
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Rule Viewer 45, 47, 59 
ruleedit 3-59 
ruleview 3-61 
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showrule 3-68 
sigmf 27, 3-70 
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101 
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